Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F50 I1 R3 779 308.9 10220750 96.7% 9883465 138.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP025268 3,743,290 +C intergenic (‑66/‑116) yiaJ ← / → galK IclR family transcriptional regulator YiaJ/galactokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0252683,743,2891.C99.1% 372.3 / ‑3.3 115intergenic (‑65/‑117)yiaJ/galKIclR family transcriptional regulator YiaJ/galactokinase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (64/50);  total (64/51)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 121 total aligned reads displayed.

ACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCATCAATCCTGTTGACAATTAATCATCGGCATAGTATATCGGCATAGTATAATACGACAAGGTGAGGAACTAAACCCAGGAGGCAGATCATGAGTCTGAAAGAAAAA  >  NZ_CP025268/3743150‑3743423
                                                                                                                                             |                                                                                                                                      
aCTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATaa                                                                                                                                         >  1:3888623/1‑140 (MQ=255)
aCTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATaa                                                                                                                                         >  1:3803514/1‑140 (MQ=255)
 ctctGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAc                                                                                                                                       <  1:1412387/140‑1 (MQ=255)
  tctGGCTTCCGGCTGGACGCTCTTGTTCCTGCGCCATCTCTTTCTCTTTTTTTCCCTTCACTTCTTTCCCCATTGTGTGGCGGCCTGATGGGCGCGCAAAGTTGGCCGTAGTTCACGATCTCGATGGATGATTTAGAAcc                                                                                                                                      >  2:2816221/1‑139 (MQ=255)
     ggCTTCCGGCTTGACGGTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCAACCCTTATTGCCCCATTTTGGCGCGTCCTGATGGTGGAGCAAAGTGTGCAGTAGGTCACGAACTAGAAAGATAAGTTATAActaat                                                                                                                                   >  2:60688/1‑136 (MQ=255)
          ccGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGAtt                                                                                                                               >  1:729003/1‑140 (MQ=255)
          ccGGCTGGACGCTCTGTTTCCTGCGCCGTCTCGTTCTGTTTTTCTCCCGTGACTTCTTTACCCATTTTGTCGCGTCCTGATGGGAGCGCAAAGGGTGCCGTAGGTCACGCTCCTGGCAGATAATTGAGAACCAATTGagt                                                                                                                               >  2:3174960/1‑138 (MQ=255)
           cGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGAttt                                                                                                                              >  1:1053564/1‑140 (MQ=255)
               tGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTAt                                                                                                                          <  1:4693082/140‑1 (MQ=255)
                ggACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATg                                                                                                                         >  2:1492761/1‑140 (MQ=255)
                         ttttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGaa                                                                                                                <  1:3381819/140‑1 (MQ=255)
                         ttttCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGaa                                                                                                                <  2:4604877/140‑1 (MQ=255)
                            tCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAAtt                                                                                                             <  1:2994157/140‑1 (MQ=255)
                               tGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcat                                                                                                          <  1:458177/140‑1 (MQ=255)
                                gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcatc                                                                                                         >  1:1673632/1‑140 (MQ=255)
                                gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcatc                                                                                                         <  2:258294/140‑1 (MQ=255)
                                gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcatc                                                                                                         >  1:4277869/1‑140 (MQ=255)
                                gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcatc                                                                                                         >  2:4477226/1‑140 (MQ=255)
                                gcgcCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATtcatc                                                                                                         <  1:679228/140‑1 (MQ=255)
                                      tctcGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCag                                                                                                   >  1:2955368‑M1/1‑138 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  2:1263929‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  1:3474799‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  1:4010759‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  1:2376907‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  2:3786616‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  1:4278471‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  2:1325138‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  2:1897379‑M1/140‑9 (MQ=255)
                                            tctcTTTTTTTCCAATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcc                                                                                             <  2:2004775‑M1/140‑9 (MQ=255)
                                             ctctTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcct                                                                                            <  1:1155665‑M1/140‑10 (MQ=255)
                                             ctctTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgcct                                                                                            <  1:2441791‑M1/140‑10 (MQ=255)
                                                 ttttttCCCATCACTTCTTTCACCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgtt                                                                                        <  2:2027612‑M1/140‑14 (MQ=255)
                                                       cccATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaag                                                                                  <  2:4626246‑M1/140‑20 (MQ=255)
                                                        ccATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagt                                                                                 <  1:1973658‑M1/140‑21 (MQ=255)
                                                          aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTTATTTTTATGTCTTTTGAAATTCATCAATGagattgccttgttaaaaagtga                                                                               >  2:481916‑M1/1‑118 (MQ=255)
                                                          aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtga                                                                               <  2:243805‑M1/140‑23 (MQ=255)
                                                          aTCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtga                                                                               >  1:4998032‑M1/1‑118 (MQ=255)
                                                           tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATACCCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgat                                                                              <  1:1137701‑M1/140‑24 (MQ=255)
                                                           tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgat                                                                              <  2:2391483‑M1/140‑24 (MQ=255)
                                                           tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgat                                                                              <  1:1723183‑M1/140‑24 (MQ=255)
                                                           tCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgat                                                                              <  1:2955731‑M1/140‑24 (MQ=255)
                                                                  tttCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatata                                                                       >  2:3652474‑M1/1‑110 (MQ=255)
                                                                      cccATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                                                   >  2:596007‑M1/1‑106 (MQ=255)
                                                                      cccATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                                                   <  1:4941367‑M1/140‑35 (MQ=255)
                                                                      cccATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                                                   >  2:595504‑M1/1‑106 (MQ=255)
                                                                      cccATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGGTCACGATCTCGACAGATAATTTATAACCAATTGATTTCTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                                                   <  1:4150417‑M1/140‑35 (MQ=255)
                                                                         aTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaa                                                                >  1:573562‑M1/1‑103 (MQ=255)
                                                                          ttttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaat                                                               <  1:1114045‑M1/140‑39 (MQ=255)
                                                                           tttGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatc                                                              >  1:3656283‑M1/1‑101 (MQ=255)
                                                                                      cTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcat                                                   >  2:1393131‑M1/1‑90 (MQ=255)
                                                                                          tGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattg                                               <  1:4810570‑M1/140‑55 (MQ=255)
                                                                                             tAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaa                                            >  1:191651‑M1/1‑83 (MQ=255)
                                                                                              aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaat                                           >  2:2758963‑M1/1‑82 (MQ=255)
                                                                                              aGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATCTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaat                                           >  1:1228904‑M1/1‑82 (MQ=255)
                                                                                               gcgcAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaatt                                          >  1:4263155‑M1/1‑81 (MQ=255)
                                                                                                  cAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaatttta                                       <  1:1216273‑M1/140‑63 (MQ=255)
                                                                                                  aaaagTGTGCCGTAGTTCACGATCTCGACAGATAATTTAGAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaatttta                                       <  1:3830795‑M1/139‑63 (MQ=255)
                                                                                                   aaaGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaa                                      >  1:1963141‑M1/1‑77 (MQ=255)
                                                                                                     aGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagc                                    >  1:1524699‑M1/1‑75 (MQ=255)
                                                                                                         tgCCGTAGTTCACGATCTCGACAGATAATTTATAACCACTTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaa                                >  2:1830733‑M1/1‑71 (MQ=255)
                                                                                                         tgCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGAGTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaa                                >  2:1515215‑M1/1‑71 (MQ=255)
                                                                                                               aGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcg                          <  2:1853022‑M1/140‑76 (MQ=255)
                                                                                                                gTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcga                         >  2:641213‑M1/1‑64 (MQ=255)
                                                                                                                 ttCACGATCTCGCCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaatgttaagccaaaaaagcgat                        >  2:1446434‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  1:4485326‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  2:2877732‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  2:4481505‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  2:2763296‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  1:3588671‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  1:4026836‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  1:2387545‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  1:512802‑M1/1‑63 (MQ=255)
                                                                                                                 ttCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgat                        >  2:1446433‑M1/1‑63 (MQ=255)
                                                                                                                   cACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatca                      >  1:3416601‑M1/1‑61 (MQ=255)
                                                                                                                    aCGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaa                     <  1:1501898‑M1/140‑81 (MQ=255)
                                                                                                                    aCGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaa                     <  2:1301333‑M1/140‑81 (MQ=255)
                                                                                                                    aCGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaa                     <  2:1301755‑M1/140‑81 (MQ=255)
                                                                                                                     cGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaa                    >  2:4761011‑M1/1‑59 (MQ=255)
                                                                                                                        tctcGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaa                 >  1:2241057‑M1/1‑56 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttagaaagtgatcgatatatttgaaatcaagtttcgcatatggaaattttaagccaaaaaagcgatcaaaaaaaca              >  2:4532541‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  2:4104634‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  1:1087328‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  1:214759‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  2:4906653‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  1:4626807‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTAGGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagctatcaaaaaaaca              >  2:3882901‑M1/1‑53 (MQ=255)
                                                                                                                           cGACAGATAATTTATAACCAATTGATTTTTAAGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaaca              >  2:3882913‑M1/1‑53 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  1:3349513‑M1/140‑90 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  1:190319‑M1/140‑90 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  2:1769017‑M1/140‑90 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  1:1485261‑M1/140‑90 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgacatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  2:1622259‑M1/140‑90 (MQ=255)
                                                                                                                             aCAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccctgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaag            <  1:4076022‑M1/140‑90 (MQ=255)
                                                                                                                               aGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaagga          <  1:232404‑M1/140‑92 (MQ=255)
                                                                                                                                gATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgataaaaaaaacaaggaa         >  1:4613418‑M1/1‑48 (MQ=255)
                                                                                                                                 aTAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaaggaag        >  1:3338567‑M1/1‑47 (MQ=255)
                                                                                                                                   aaTTTATAACCAATTGATTTGTATGTCTTTTGACATTCATCAATCagattgccttgttaaaaagtgatcgatatattagaaatcaaggttcgcatattgaaattgtaagccaaaccagcgatcaaaaaaacaaggaagcc      <  2:4213836‑M1/140‑96 (MQ=255)
                                                                                                                                    aTTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaaggaagcct     >  1:1309496‑M1/1‑44 (MQ=255)
                                                                                                                                    aTTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaaaagcgatcaaaaaaacaaggaagcct     >  1:1308622‑M1/1‑44 (MQ=255)
                                                                                                                                       tataACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgagattttaagccaaaaaagcgatcaaaaaaacaaggaagcctggg  <  1:2989633‑M1/140‑100 (MQ=255)
                                                                                                                                             |                                                                                                                                      
ACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCATCAATCCTGTTGACAATTAATCATCGGCATAGTATATCGGCATAGTATAATACGACAAGGTGAGGAACTAAACCCAGGAGGCAGATCATGAGTCTGAAAGAAAAA  >  NZ_CP025268/3743150‑3743423

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCATCAATCCTGTTGACAATTAATCATCGGCATAGTATATCGGCATAGTATAATACGACAAGGTGAGGAACTAAACCCA  >  NZ_CP025268/3743140‑3743394
                                                                                                                                                        |                                                                                                        
CGCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC                                                                                                          <  SRR14920166.2316665/151‑1 (MQ=60)
CGCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC                                                                                                          <  SRR14920166.378218/151‑1 (MQ=60)
CGCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAAC                                                                                                          <  SRR14920166.1814637/151‑1 (MQ=60)
 gctcAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCCTTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACC                                                                                                         <  SRR14920166.54412/147‑1 (MQ=60)
 GCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACC                                                                                                         <  SRR14920166.3459344/151‑1 (MQ=60)
 GCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACC                                                                                                         <  SRR14920166.676809/151‑1 (MQ=60)
 GCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTATTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGCTAATTTATAACC                                                                                                         <  SRR14920166.4084871/151‑1 (MQ=60)
    AAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAAT                                                                                                       >  SRR14920166.4173879/1‑151 (MQ=60)
         GACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATT                                                                                                  >  SRR14920166.730635/1‑151 (MQ=60)
         GACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATT                                                                                                  >  SRR14920166.244395/1‑151 (MQ=60)
         GACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATT                                                                                                  <  SRR14920166.2359768/151‑1 (MQ=60)
          ACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTT                                                                                                 >  SRR14920166.1055770/1‑151 (MQ=60)
           CTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTT                                                                                                <  SRR14920166.1415249/151‑1 (MQ=60)
                 CTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTC                                                                                          >  SRR14920166.1184435/1‑151 (MQ=60)
                     CGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTT                                                                                      >  SRR14920166.4955346/1‑151 (MQ=60)
                         TGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAA                                                                                  <  SRR14920166.4703361/151‑1 (MQ=60)
                               CTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATC                                                                            >  SRR14920166.1677049/1‑151 (MQ=60)
                               CTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATC                                                                            >  SRR14920166.4287251/1‑151 (MQ=60)
                                 CTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAA                                                                          >  SRR14920166.3989100/1‑151 (MQ=60)
                                   TTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.3389023/151‑1 (MQ=60)
                                     TTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCag                                                                      >  SRR14920166.2961650/1‑149 (MQ=60)
                                      TCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCaga                                                                     <  SRR14920166.3000504/151‑4 (MQ=60)
                                         TGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattg                                                                  <  SRR14920166.459239/151‑7 (MQ=60)
                                          GCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgc                                                                 <  SRR14920166.680764/151‑8 (MQ=60)
                                                 CTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgtta                                                          >  SRR14920166.1437597/1‑137 (MQ=60)
                                                    GTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaa                                                       >  SRR14920166.4614656/1‑134 (MQ=60)
                                                      TCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaag                                                     <  SRR14920166.3482136/151‑20 (MQ=60)
                                                      TCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaag                                                     <  SRR14920166.4019541/151‑20 (MQ=60)
                                                      TCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaag                                                     <  SRR14920166.2381737/151‑20 (MQ=60)
                                                      TCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaag                                                     <  SRR14920166.4287853/151‑20 (MQ=60)
                                                       CTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagt                                                    <  SRR14920166.2446758/151‑21 (MQ=60)
                                                       CTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagc                                                    <  SRR14920166.1158074/151‑21 (MQ=60)
                                                         CTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtga                                                  >  SRR14920166.5008726/1‑129 (MQ=60)
                                                                TCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatat                                           >  SRR14920166.2190152/1‑122 (MQ=60)
                                                                  CCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatat                                         <  SRR14920166.1977796/151‑32 (MQ=60)
                                                                     TCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATACCCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                      <  SRR14920166.1140080/151‑35 (MQ=60)
                                                                     TCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                      <  SRR14920166.3329453/151‑35 (MQ=60)
                                                                     TCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                      <  SRR14920166.2962015/151‑35 (MQ=60)
                                                                     TCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttg                                      <  SRR14920166.1726685/151‑35 (MQ=60)
                                                                        CTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaa                                   >  SRR14920166.574872/1‑114 (MQ=60)
                                                                          TCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatc                                 >  SRR14920166.3664106/1‑112 (MQ=60)
                                                                                CCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagttt                           <  SRR14920166.4951971/151‑46 (MQ=60)
                                                                                CCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGGTCACGATCTCGACAGATAATTTATAACCAATTGATTTCTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagttt                           <  SRR14920166.4159502/151‑46 (MQ=60)
                                                                                    TTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgca                       <  SRR14920166.1116379/151‑50 (MQ=60)
                                                                                            CGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaa               >  SRR14920166.192129/1‑94 (MQ=60)
                                                                                             GTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATCTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaat              >  SRR14920166.1231417/1‑93 (MQ=60)
                                                                                              TCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaatt             >  SRR14920166.4272510/1‑92 (MQ=60)
                                                                                                  TATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaa         >  SRR14920166.1967215/1‑88 (MQ=60)
                                                                                                    TGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagc       <  SRR14920166.4820993/151‑66 (MQ=60)
                                                                                                    TGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagc       >  SRR14920166.1527802/1‑86 (MQ=60)
                                                                                                         GCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATCagattgccttgttaaaaagtgatcgatatatttgaaatcaagtttcgcatattgaaattttaagccaaaa  >  SRR14920166.356105/1‑81 (MQ=60)
                                                                                                            CAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.1218769/78‑1 (MQ=60)
                                                                                                            CAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.1758557/78‑1 (MQ=60)
                                                                                                            AAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTAGAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.3839162/78‑1 (MQ=60)
                                                                                                                GTGTGCGGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.3596184/1‑74 (MQ=60)
                                                                                                                GTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.513962/1‑74 (MQ=60)
                                                                                                                GTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.4035636/1‑74 (MQ=60)
                                                                                                                GTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.4495153/1‑74 (MQ=60)
                                                                                                                GTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.2392393/1‑74 (MQ=60)
                                                                                                                  GTGCCGTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.3423849/1‑72 (MQ=60)
                                                                                                                       GTAGTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.2245642/1‑67 (MQ=60)
                                                                                                                          GTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.4636980/1‑64 (MQ=60)
                                                                                                                          GTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.1089584/1‑64 (MQ=60)
                                                                                                                          GTTCACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.215293/1‑64 (MQ=60)
                                                                                                                              ACGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.1504951/60‑1 (MQ=60)
                                                                                                                               CGATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.4623564/1‑59 (MQ=60)
                                                                                                                                GATCTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.3345699/1‑58 (MQ=60)
                                                                                                                                   CTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.1312132/1‑55 (MQ=60)
                                                                                                                                   CTCGACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.1311258/1‑55 (MQ=60)
                                                                                                                                       ACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.4084930/51‑1 (MQ=60)
                                                                                                                                       ACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.1488275/51‑1 (MQ=60)
                                                                                                                                       ACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.190794/51‑1 (MQ=60)
                                                                                                                                       ACAGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.3356659/51‑1 (MQ=60)
                                                                                                                                         AGATAATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.232971/49‑1 (MQ=60)
                                                                                                                                              ATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.1810923/1‑44 (MQ=60)
                                                                                                                                              ATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.3672229/1‑44 (MQ=60)
                                                                                                                                              ATTTATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.4963752/1‑44 (MQ=60)
                                                                                                                                                 TATAACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.2995974/41‑1 (MQ=60)
                                                                                                                                                  ATAACCAATTGATCTTTATGTCTTTTGAAATTCATCAATC                                                                        <  SRR14920166.4487037/40‑1 (MQ=54)
                                                                                                                                                    ACCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.3319021/1‑37 (MQ=60)
                                                                                                                                                      CCAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.1363887/1‑36 (MQ=60)
                                                                                                                                                        CAATTGATTTTTATGTCTTTTGAAATTCATCAATC                                                                        >  SRR14920166.182675/1‑35 (MQ=60)
                                                                                                                                                        |                                                                                                        
CGCGAAACAGACTCTGGCTTCCGGCTGGACGCTCTTTTTCCTGCGCCATCTCGTTCTCTTTTTTTCCCATCACTTCTTTCCCCATTTTGTCGCGTCCTGATGGTAGCGCAAAGTGTGCCGTAGTTCACGATCTCGACAGATAATTTATAACAATTGATTTTTATGTCTTTTGAAATTCATCAATCCTGTTGACAATTAATCATCGGCATAGTATATCGGCATAGTATAATACGACAAGGTGAGGAACTAAACCCA  >  NZ_CP025268/3743140‑3743394

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: