Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,360,157 C→T 29.8% I138I (ATC→ATT CJ019_03160 → Tn3 family transposase ISPa42

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,360,1570CT29.8% 98.7 / 54.2 67I138I (ATC→ATTCJ019_03160Tn3 family transposase ISPa42
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (13/7);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-02

CCCTTTACCGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAA  >  ACB143/3360024‑3360290
                                                                                                                                     |                                                                                                                                     
cccTTTACCGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTaa                                                                                                                                     <  4:26316/136‑1 (MQ=255)
cccTTTACCGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCaa                                                                                                                                     <  6:56839/136‑1 (MQ=255)
      aCCGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCAt                                                                                                                               >  8:125813/1‑136 (MQ=255)
        cGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGc                                                                                                                             >  2:216478/1‑136 (MQ=255)
             ccAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCTGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTtgat                                                                                                                        <  6:194537/136‑1 (MQ=255)
             ccAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTtgat                                                                                                                        <  1:150829/136‑1 (MQ=255)
              cAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTtgatg                                                                                                                       >  4:81151/1‑136 (MQ=255)
                   agcGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTg                                                                                                                  <  1:216478/136‑1 (MQ=255)
                        tACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCa                                                                                                             <  3:288897/136‑1 (MQ=255)
                        tACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCa                                                                                                             <  5:83142/136‑1 (MQ=255)
                        tACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCa                                                                                                             <  3:17827/136‑1 (MQ=255)
                        tACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCa                                                                                                             <  2:257923/136‑1 (MQ=255)
                           ggCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCg                                                                                                          >  7:304484/1‑136 (MQ=255)
                             cAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTg                                                                                                        >  4:72774/1‑136 (MQ=255)
                                     gCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAg                                                                                                >  1:93911/1‑136 (MQ=255)
                                      cAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAg                                                                                                <  5:51943/135‑1 (MQ=255)
                                      cAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAg                                                                                                <  7:7140/135‑1 (MQ=255)
                                        ggTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCtt                                                                                             >  5:289032/1‑136 (MQ=255)
                                        ggTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCtt                                                                                             >  6:143188/1‑136 (MQ=255)
                                             cccccTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGc                                                                                        >  1:30885/1‑136 (MQ=255)
                                               cccTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCt                                                                                       >  1:148175/1‑135 (MQ=255)
                                                    cTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTAtt                                                                                   <  8:304484/134‑1 (MQ=255)
                                                    cTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTAtt                                                                                   <  5:54251/134‑1 (MQ=255)
                                                    cTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTAtt                                                                                   <  6:188886/134‑1 (MQ=255)
                                                    cTACAGCAGCGACGGGCTAGAGGTCACTGAGCAGATCGCACTTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTAtt                                                                                   <  7:141577/134‑1 (MQ=255)
                                                      aCAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGAc                                                                               >  1:134737/1‑136 (MQ=255)
                                                      aCAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGAc                                                                               >  3:86888/1‑136 (MQ=255)
                                                       cagcagCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCg                                                                                 <  1:259271/133‑1 (MQ=255)
                                                       cagcagCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCg                                                                                 <  5:177345/133‑1 (MQ=255)
                                                             cGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTcaa                                                                         >  5:8138/1‑135 (MQ=255)
                                                             cGACGGGCTA‑‑GGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTcaacaa                                                                      >  7:54007/1‑136 (MQ=255)
                                                               aCGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTcaacaa                                                                      >  5:245929/1‑136 (MQ=255)
                                                               aCGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTcaacaa                                                                      >  7:138023/1‑136 (MQ=255)
                                                               aCGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTcaaca                                                                       >  6:167868/1‑135 (MQ=255)
                                                                cGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCt                                                                            >  4:80461/1‑129 (MQ=255)
                                                                cGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCt                                                                            <  3:80461/129‑1 (MQ=255)
                                                                  ggCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCt                                                                   >  4:69809/1‑136 (MQ=255)
                                                                       gCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATGTCATTAATTCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCCACAACCACg                                                                 >  3:222810/1‑133 (MQ=255)
                                                                                gAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCt                                                     >  6:171643/1‑136 (MQ=255)
                                                                                gAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGCTTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCt                                                     >  8:212809/1‑136 (MQ=255)
                                                                                 aGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCtt                                                    <  3:72774/136‑1 (MQ=255)
                                                                                    aGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGc                                                 <  1:49822/136‑1 (MQ=255)
                                                                                    aGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGc                                                 <  6:216159/136‑1 (MQ=255)
                                                                                        cGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTAAGCGTTGCTATTAGCTTACGGCTTATTCGACGATCAACAACCTAGCAGAAACCGCCTTGGCGAGt                                             >  4:151834/1‑136 (MQ=255)
                                                                                         gCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCg                                                                                 >  8:237588/1‑99 (MQ=255)
                                                                                         gCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCg                                                                                 <  7:237588/99‑1 (MQ=255)
                                                                                           aCGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGcc                                          >  1:10080/1‑136 (MQ=255)
                                                                                              tGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAg                                       >  8:221887/1‑136 (MQ=255)
                                                                                                aCTTCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGa                                     <  2:319206/136‑1 (MQ=255)
                                                                                                   gcgcTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGt                                  >  5:227798/1‑136 (MQ=255)
                                                                                                    cgcTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGtt                                 <  8:54007/136‑1 (MQ=255)
                                                                                                    cgcTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTAGGCGAGTGCCCAGGAAGt                                  <  6:292767/135‑1 (MQ=255)
                                                                                                     gcTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTa                                >  2:83582/1‑136 (MQ=255)
                                                                                                     gcTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTa                                <  8:240522/136‑1 (MQ=255)
                                                                                                      cTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGTTGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTAc                               <  5:147919/136‑1 (MQ=255)
                                                                                                      cTGTAGTTCTTGAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTAc                               <  2:42271/136‑1 (MQ=255)
                                                                                                      cTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCTGCTTATTCGACGCTCAACAACCTCGAAGAAACCGACTTGGCGAGTGACCAGTAAGTTAc                               <  4:90426/136‑1 (MQ=255)
                                                                                                      cTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTAc                               <  5:265878/136‑1 (MQ=255)
                                                                                                          aGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCtt                            >  1:72184/1‑135 (MQ=255)
                                                                                                           gTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCa                          <  5:167868/136‑1 (MQ=255)
                                                                                                           gTTCTTGAGCAGCGAGGCGATCTGATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAACCCGCCTTGGCGAGTGCCCATGAAGTTACCTTCa                          <  6:129142/136‑1 (MQ=255)
                                                                                                              cTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAAcc                       >  2:109210/1‑136 (MQ=255)
                                                                                                                 gAGCAGCGAGTCGATCTCATTAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAg                    >  5:76576/1‑136 (MQ=255)
                                                                                                                   gcagcGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTc                  >  1:231227/1‑136 (MQ=255)
                                                                                                                       cGAGTCGATCGCATTAATGGCATGCTTGATGAGCTGATTCAGCGTTGCTATGAGCTTCCCGCTTTTTCGACGCTCAACAACCTAGCAGAAACCGCCTTGGCGAGTGCCCAGGAACGTACCTTCAACCTAGTCGTGa              >  6:21353/1‑136 (MQ=255)
                                                                                                                         aGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACt            >  1:173810/1‑136 (MQ=255)
                                                                                                                           tCGATCTGATTAATGCCATTCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCg          <  6:8138/136‑1 (MQ=255)
                                                                                                                                 tcatcaATGCCATGCTTGAAGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTAATCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCGacccc     >  2:183575/1‑131 (MQ=255)
                                                                                                                                   atcaATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCGAgcgc     <  3:290484/133‑1 (MQ=255)
                                                                                                                                   atcaATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCaa  <  4:86888/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CCCTTTACCGGCACCAAGCAGCGGTACGGCAGTACCTGCAGGTTACCCCCTTCTACAGCAGCGACGGGCTAGCGGTCACTGAGCAGATCGCACGTGACTGCGCTGTAGTTCTTGAGCAGCGAGTCGATCTCATCAATGCCATGCTTGATGAGCTGATTCAGCGTGGCTATGAGCTTCCGGCTTATTCGACGCTCAACAACCTCGCAGAAACCGCCTTGGCGAGTGCCCAGGAAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAA  >  ACB143/3360024‑3360290

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: