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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,836,192 T→A 90.7% I67I (ATA→ATT)  ynjH ← DUF1496 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,836,1920TA90.7% 380.1 / 24.0 128I67I (ATA→ATT) ynjHDUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base A (50/66);  total (52/76)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGT  >  CP009273/1836050‑1836337
                                                                                                                                              |                                                                                                                                                 
gtatggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTTACATTCAATTATTAGCGGCGTTTATGTCTTTACGCGCCGACAGACCAGCGGATTAGTGCTAAGGGTTTTGTCTTCACGCTGGCATTGCAGCAGAATTcc                                                                                                                                              <  2:1679305/145‑1 (MQ=255)
 tctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCt                                                                                                                                             <  1:2589870/148‑1 (MQ=255)
 tctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCt                                                                                                                                             <  1:885027/148‑1 (MQ=255)
 tctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTATTGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTGATCACGCTGGCATTGCAGCAGAATTCCt                                                                                                                                             >  1:110277/1‑148 (MQ=255)
     gtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCgg                                                                                                                                         >  2:2632633/1‑148 (MQ=255)
     gtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCgg                                                                                                                                         <  1:658536/148‑1 (MQ=255)
       gCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCt                                                                                                                                       <  2:1450750/148‑1 (MQ=255)
        cTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCtt                                                                                                                                      <  1:690357/148‑1 (MQ=255)
         tCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTTTGGCTTTAAGCGCCGCCCCACCAGCGGCTTAGTGCTAAGGGGTTTGTGATCACGCAGGCAATGGAGCAGAAATCCTTCGGCttt                                                                                                                                     >  1:3147649/1‑148 (MQ=255)
          cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                                    <  1:1580933/148‑1 (MQ=255)
          cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                                    <  1:2852380/148‑1 (MQ=255)
          cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTa                                                                                                                                    <  1:2054808/148‑1 (MQ=255)
           aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTaa                                                                                                                                   >  2:2100861/1‑148 (MQ=255)
           aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTaa                                                                                                                                   >  2:1653686/1‑148 (MQ=255)
           aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTaa                                                                                                                                   <  1:1085375/148‑1 (MQ=255)
            ccTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAt                                                                                                                                  <  2:1944871/148‑1 (MQ=255)
              tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCCGACAAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTACGTCGGCGTTAATTa                                                                                                                                >  2:3119851/1‑148 (MQ=255)
              tGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTa                                                                                                                                <  1:1975108/148‑1 (MQ=255)
               gAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTCCAGCAGAATTCCTTCGGCTTTAATTAc                                                                                                                               <  1:1805089/148‑1 (MQ=255)
                 agagGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCg                                                                                                                             <  1:966573/148‑1 (MQ=255)
                  gagGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGc                                                                                                                            <  2:2848725/148‑1 (MQ=255)
                  gagGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGc                                                                                                                            <  2:2850634/148‑1 (MQ=255)
                    ggCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGccc                                                                                                                          <  2:3049581/148‑1 (MQ=255)
                    ggCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGGGTTTTGGTTTAACCGGCGCCAGACCAGCGGCTTAGTGGAAAGGGGGTTGTGATCACGCGTGCCTTTCAGCAGAATTACTTCGGGTTTAATTaaacccc                                                                                                                          >  1:2659327/1‑142 (MQ=255)
                     gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTATGGCGTTTCGCGCCGCCAGACCAGCGGCTTAGTGGTAAGGGGTTTGTGATCACGCGGGGATTGCAGCAGTATTGCTTCGGGGTTTATTAAcccccc                                                                                                                         >  1:1185402/1‑148 (MQ=255)
                        cTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGAGTAGTGCTAAGGGTTTTGTCATCACGCTGGCCTTGCACCAGTATTCCTTCGGCTTTAATTACCGCCCCTTc                                                                                                                      >  2:2926308/1‑148 (MQ=255)
                          gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTTTGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTAACGCCCCTTccg                                                                                                                    >  2:1446634/1‑146 (MQ=255)
                           cAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCGTTAATTACCGCCCCTTCAGa                                                                                                                   >  2:797218/1‑148 (MQ=255)
                            aaTATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGaa                                                                                                                  <  2:660245/148‑1 (MQ=255)
                              tatCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGaata                                                                                                                <  2:813899/148‑1 (MQ=255)
                               atCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTAGGCTTTAATTACCGCCCCGTCAGaaaaa                                                                                                               >  1:1623646/1‑148 (MQ=255)
                                tCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTTTGGCTTTTCGCGCCGCCAGACCAGCGGATTAGTGGTAAGGGTTTTGTCATAACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTTCCGCCCCGTCAGaataat                                                                                                              >  2:397718/1‑148 (MQ=255)
                                 ccGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                                                                             <  1:2541499/148‑1 (MQ=255)
                                 ccGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                                                                             <  1:1805355/148‑1 (MQ=255)
                                 ccGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                                                                             <  2:2995658/148‑1 (MQ=255)
                                 ccGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAtt                                                                                                             <  1:386647/148‑1 (MQ=255)
                                  cGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAttt                                                                                                            <  2:1796324/148‑1 (MQ=255)
                                     tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGa                                                                                                         <  2:2465586/148‑1 (MQ=255)
                                       gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTTAta                                                                                                       >  2:3265468/1‑147 (MQ=255)
                                       gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                                                                       <  1:1815737/148‑1 (MQ=255)
                                       gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                                                                                       <  2:3206873/148‑1 (MQ=255)
                                       gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCGTTACGCGCCGCCACACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCGGCAGAATTACTTCGGCTTTAATTACCGCCCCAACAGAATAATTTTGATc                                                                                                       >  1:2943369/1‑148 (MQ=255)
                                          ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTg                                                                                                    <  2:2627123/148‑1 (MQ=255)
                                          ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCACACCAGCGGATTAGTGCTAAGGGTTTTGTGATCACGCGGGCATTGCAGCAGAAATACTTCGGCTTTAATTACCCCCCCTTCAGAAGAAATTTGATCCTg                                                                                                    >  1:2907283/1‑148 (MQ=255)
                                              ctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATaa                                                                                                <  2:73643/148‑1 (MQ=255)
                                                 aCATTCCATTATTAGAGGCGTTTAAGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGGTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTGCGGGGTTAATTACCGCCCCGTCAGCATAAATTGGTTCGTGCTTAACg                                                                                             >  2:2257222/1‑148 (MQ=255)
                                                  cATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgc                                                                                            <  1:3291975/148‑1 (MQ=255)
                                                  cATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgc                                                                                            <  1:1904548/148‑1 (MQ=255)
                                                  cATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgc                                                                                            <  1:2371980/148‑1 (MQ=255)
                                                           attaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAg                                                                                   <  1:2100861/148‑1 (MQ=255)
                                                           attaGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCGTTAATTACCGCCCCATCAGAATAATTTTGATCCTGATAAAAGCAGCACagag                                                                                   >  2:1927251/1‑148 (MQ=255)
                                                               gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACa                                                                               <  1:286271/148‑1 (MQ=255)
                                                                 ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAgg                                                                             <  1:2109905/148‑1 (MQ=255)
                                                                  gCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAggg                                                                            <  2:909333/148‑1 (MQ=255)
                                                                  gCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAggg                                                                            <  2:1670115/148‑1 (MQ=255)
                                                                   cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGc                                                                           <  2:1067874/148‑1 (MQ=255)
                                                                   cGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGc                                                                           <  2:2719518/148‑1 (MQ=255)
                                                                      ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgc                                                                        >  1:2485941/1‑148 (MQ=255)
                                                                       tATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTgcg                                                                       <  2:2781330/148‑1 (MQ=255)
                                                                           gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTATGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTAAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGa                                                                   <  1:18919/148‑1 (MQ=255)
                                                                           gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGa                                                                   <  2:1428382/148‑1 (MQ=255)
                                                                           gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGa                                                                   <  2:1420146/148‑1 (MQ=255)
                                                                           gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGa                                                                   >  2:197911/1‑148 (MQ=255)
                                                                           gCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGa                                                                   <  2:1420148/148‑1 (MQ=255)
                                                                            cTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAc                                                                  >  1:2140795/1‑148 (MQ=255)
                                                                              ttACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                                >  1:1203803/1‑148 (MQ=255)
                                                                              ttACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg                                                                >  1:1410102/1‑148 (MQ=255)
                                                                               tACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGc                                                               <  2:2668870/148‑1 (MQ=255)
                                                                                          aGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCATCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaact                                                    >  1:1782777/1‑148 (MQ=255)
                                                                                           gACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactg                                                   >  1:505254/1‑148 (MQ=255)
                                                                                            aCCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCATGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCctgaactga                                                  >  2:729570/1‑148 (MQ=255)
                                                                                              cAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCtgaactgaac                                                <  2:2334792/148‑1 (MQ=255)
                                                                                               aGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAaca                                               <  2:2039082/148‑1 (MQ=255)
                                                                                               aGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAaca                                               <  2:2039014/148‑1 (MQ=255)
                                                                                                     ttAGTGATAAGGTTTTTGTCATCACGCTGGCATTGCAGCATAATACCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTTATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTg                                         <  1:2772452/148‑1 (MQ=255)
                                                                                                      tAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTgg                                        <  2:1768578/148‑1 (MQ=255)
                                                                                                          gCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGt                                    >  2:408684/1‑148 (MQ=255)
                                                                                                          gCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGt                                    >  2:280369/1‑148 (MQ=255)
                                                                                                          gCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGt                                    >  1:2357938/1‑148 (MQ=255)
                                                                                                                ggTTTTGTCATCACGCTGGCATTGCAGCAGAATACCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGTCGGTACGttt                              >  2:2852383/1‑148 (MQ=255)
                                                                                                                 gTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTa                             >  1:2011295/1‑148 (MQ=255)
                                                                                                                  ttttGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTAc                            <  1:513996/148‑1 (MQ=255)
                                                                                                                      gTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTcc                        <  1:2845536/148‑1 (MQ=255)
                                                                                                                         atcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGGACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACg                     >  2:2873245/1‑148 (MQ=255)
                                                                                                                          tcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTCCGTTTACCTCCACGt                    >  1:1650997/1‑148 (MQ=255)
                                                                                                                          tcaCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGt                    >  1:1911728/1‑148 (MQ=255)
                                                                                                                          tcaCGCTGGCAGTGCAGCATAATTCCTTCGGCTTTAATTACCCCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGACCTGAACACTTCTGGCGGTACGTTTACCTCCACGt                    <  1:3087223/148‑1 (MQ=255)
                                                                                                                           caCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc                   >  1:3018529/1‑148 (MQ=255)
                                                                                                                           caCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc                   >  1:2482233/1‑148 (MQ=255)
                                                                                                                           caCGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTc                   >  1:2214302/1‑148 (MQ=255)
                                                                                                                              gCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCgg                <  2:1447984/148‑1 (MQ=255)
                                                                                                                              gCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCgg                >  2:335667/1‑148 (MQ=255)
                                                                                                                                  gCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGa            >  1:192314/1‑148 (MQ=255)
                                                                                                                                     ttGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATaa         >  1:2350177/1‑148 (MQ=255)
                                                                                                                                     ttGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATaa         >  1:2350166/1‑148 (MQ=255)
                                                                                                                                       gcagcaGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATg       <  2:924648/148‑1 (MQ=255)
                                                                                                                                       gcagcaGAATTCCTTCGGATTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATg       <  1:2087745/148‑1 (MQ=255)
                                                                                                                                         agcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCg     <  2:631023/148‑1 (MQ=255)
                                                                                                                                          gcagAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCgg    <  2:838826/148‑1 (MQ=255)
                                                                                                                                            agTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGt  <  1:2181863/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGT  >  CP009273/1836050‑1836337

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGC  >  CP009273/1836043‑1836341
                                                                                                                                                     |                                                                                                                                                     
GGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGA                                                                                                                                                       <  SRR19973070.2527996/150‑1 (MQ=60)
  CGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAT                                                                                                                                                     <  SRR19973070.872119/150‑1 (MQ=60)
  CGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGACCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAAT                                                                                                                                                     <  SRR19973070.1724259/150‑1 (MQ=60)
   GCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATT                                                                                                                                                    <  SRR19973070.1090295/150‑1 (MQ=60)
   GCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTAGCATTGCAGCAGAATT                                                                                                                                                    <  SRR19973070.807601/150‑1 (MQ=60)
      GGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTATTGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTGATCACGCTGGCATTGCAGCAGAATTCCT                                                                                                                                                 >  SRR19973070.111301/1‑150 (MQ=60)
        TCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTC                                                                                                                                               <  SRR19973070.2476623/150‑1 (MQ=60)
        TCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTC                                                                                                                                               <  SRR19973070.891265/150‑1 (MQ=60)
        TCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTC                                                                                                                                               <  SRR19973070.2606794/150‑1 (MQ=60)
            GTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCT                                                                                                                                           <  SRR19973070.663227/150‑1 (MQ=60)
             TGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTTTGGCTTTAAGCGCCGCCAGACCCCCGGCTTAGTGCTAATGGTTTTGTCATGACGCTGGCATTGCAGCAGAATTCCTTCGGGTT                                                                                                                                          >  SRR19973070.92078/1‑150 (MQ=60)
              GCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTTTGGCTTTAAGCGCCGCCCCACCAGCGGCTTAGTGCTAAGGGGTTTGTGATCACGCAGGCAATGGAGCAGAAATCCTTCGGCTTT                                                                                                                                         >  SRR19973070.3167989/1‑150 (MQ=60)
               CTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTA                                                                                                                                        <  SRR19973070.695288/150‑1 (MQ=60)
                 CACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAT                                                                                                                                      <  SRR19973070.2870854/150‑1 (MQ=60)
                 CACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAT                                                                                                                                      <  SRR19973070.2870829/150‑1 (MQ=60)
                 CACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAT                                                                                                                                      <  SRR19973070.1592017/150‑1 (MQ=60)
                 CACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAT                                                                                                                                      <  SRR19973070.2068692/150‑1 (MQ=60)
                 CACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAAT                                                                                                                                      <  SRR19973070.2068619/150‑1 (MQ=60)
                  ACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATT                                                                                                                                     <  SRR19973070.1092904/150‑1 (MQ=60)
                     TGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACC                                                                                                                                  <  SRR19973070.1988466/150‑1 (MQ=60)
                      GAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTCCAGCAGAATTCCTTCGGCTTTAATTACCG                                                                                                                                 <  SRR19973070.1817549/150‑1 (MQ=60)
                        AGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCC                                                                                                                               <  SRR19973070.973478/150‑1 (MQ=60)
                          AGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTGTATGGCGTTTCGCGCCGCCAGACCAGCGGCTTAGTGGTAAGGGGTTTGTGATCACGCGGGGATTGCAGCAGTATTGCTTCGGggtttattaacccccc                                                                                                                             >  SRR19973070.1193736/1‑134 (MQ=60)
                                    ATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTAGGCTTTAATTACCGCCCCGTCAGAAAAA                                                                                                                   >  SRR19973070.1634925/1‑150 (MQ=60)
                                        CCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTT                                                                                                               <  SRR19973070.2558120/150‑1 (MQ=60)
                                        CCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTT                                                                                                               <  SRR19973070.1817816/150‑1 (MQ=60)
                                        CCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTT                                                                                                               <  SRR19973070.389483/150‑1 (MQ=60)
                                        CCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTT                                                                                                               <  SRR19973070.1819743/150‑1 (MQ=60)
                                            TGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCGTTACGCGCCGCCACACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCGGCAGAATTACTTCGGCTTTAATTACCGCCCCAACAGAATAATTTTGATC                                                                                                           >  SRR19973070.2962333/1‑150 (MQ=60)
                                              GCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCT                                                                                                         <  SRR19973070.1828250/150‑1 (MQ=60)
                                               CCCCTGCTGACATTCCATTATTAGAGGCGTTTTTTGCTTTAAGCGCCCCCCAACCAGCCGAGTAGTTGTAAGGGTTTTTTGATAACCCTTGCAATGCAGCACAATACCATCGGCGTTAAATAACGCCCCGTCAGAATAATTTTGATCCTG                                                                                                        >  SRR19973070.1581797/1‑150 (MQ=14)
                                               CCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCACACCAGCGGATTAGTGCTAAGGGTTTTGTGATCACGCGGGCATTGCAGCAGAAATACTTCGGCTTTAATTACCCCCCCTTCAGAAGAAATTTGATCCTG                                                                                                        >  SRR19973070.2926084/1‑150 (MQ=60)
                                                CCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGA                                                                                                       <  SRR19973070.1665090/150‑1 (MQ=60)
                                                        ACATTCCATTATTAGAGGCGTTTTTTGGTTTACGCGCCGCCAGACCCCCGGATTAGTGCTAAGGGTTTTGTGATAACGCTGGGATTGCAGCAGAATTACATCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATACTGATAAAAGCA                                                                                               >  SRR19973070.2403461/1‑150 (MQ=60)
                                                         CATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAG                                                                                              <  SRR19973070.2387669/150‑1 (MQ=60)
                                                         CATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAG                                                                                              <  SRR19973070.1917429/150‑1 (MQ=60)
                                                         CATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAG                                                                                              <  SRR19973070.3313087/150‑1 (MQ=60)
                                                                  ATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTA                                                                                     <  SRR19973070.2114958/150‑1 (MQ=60)
                                                                      GAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGG                                                                                 <  SRR19973070.288191/150‑1 (MQ=60)
                                                                        GGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGC                                                                               <  SRR19973070.2124058/150‑1 (MQ=60)
                                                                           GTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGC                                                                            >  SRR19973070.2502098/1‑150 (MQ=60)
                                                                                 GGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGAC                                                                      >  SRR19973070.2155068/1‑150 (MQ=60)
                                                                                  GCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTATGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTAAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACT                                                                     <  SRR19973070.19195/150‑1 (MQ=60)
                                                                                   CTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG                                                                    >  SRR19973070.1212256/1‑150 (MQ=60)
                                                                                   CTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG                                                                    >  SRR19973070.1419711/1‑150 (MQ=60)
                                                                                         GCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCC                                                              <  SRR19973070.1928266/150‑1 (MQ=60)
                                                                                           GCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTG                                                            <  SRR19973070.2356184/150‑1 (MQ=60)
                                                                                               CCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCATCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACT                                                        >  SRR19973070.1795111/1‑150 (MQ=60)
                                                                                                CAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTG                                                       >  SRR19973070.508751/1‑150 (MQ=60)
                                                                                                        CGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTC                                               >  SRR19973070.1307378/1‑150 (MQ=60)
                                                                                                            TTAGTGATAAGGTTTTTGTCATCACGCTGGCATTGCAGCATAATACCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTTATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGC                                           <  SRR19973070.2790210/150‑1 (MQ=60)
                                                                                                               GTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGT                                        >  SRR19973070.2373558/1‑150 (MQ=60)
                                                                                                                      GGGTTTTGTCATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTA                                 >  SRR19973070.2024908/1‑150 (MQ=60)
                                                                                                                         TTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCT                              <  SRR19973070.517522/150‑1 (MQ=60)
                                                                                                                             GTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCAC                          <  SRR19973070.2863960/150‑1 (MQ=60)
                                                                                                                               CATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTCCGTTTACCTCCACGT                        >  SRR19973070.1662389/1‑150 (MQ=60)
                                                                                                                               CATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGT                        >  SRR19973070.1924645/1‑150 (MQ=60)
                                                                                                                                ATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTC                       >  SRR19973070.2498381/1‑150 (MQ=60)
                                                                                                                                ATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTC                       >  SRR19973070.2228852/1‑150 (MQ=60)
                                                                                                                                ATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTC                       >  SRR19973070.3037773/1‑150 (MQ=60)
                                                                                                                                ATCACGCTGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTC                       >  SRR19973070.578384/1‑150 (MQ=60)
                                                                                                                                 TCACGCTGGCAGTGCAGCATAATTCCTTCGGCTTTAATTACCCCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGACCTGAACACTTCTGGCGGTACGTTTACCTCCACGTCC                      <  SRR19973070.3107082/150‑1 (MQ=60)
                                                                                                                                       TGGCATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGA                >  SRR19973070.193795/1‑150 (MQ=60)
                                                                                                                                         GCTTTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATA              <  SRR19973070.2649799/150‑1 (MQ=60)
                                                                                                                                          CATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAA             >  SRR19973070.3271599/1‑150 (MQ=60)
                                                                                                                                          CATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAA             >  SRR19973070.2365759/1‑150 (MQ=60)
                                                                                                                                          CATTGCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAA             >  SRR19973070.2365748/1‑150 (MQ=60)
                                                                                                                                              GCAGCAGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCG         >  SRR19973070.2134875/1‑150 (MQ=60)
                                                                                                                                              GCAGCAGAATTCCTTCGGATTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCG         <  SRR19973070.2101789/150‑1 (MQ=60)
                                                                                                                                                   AGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTA    <  SRR19973070.2196296/150‑1 (MQ=60)
                                                                                                                                                   AGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTATAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTA    >  SRR19973070.501299/1‑150 (MQ=60)
                                                                                                                                                   AGAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTA    >  SRR19973070.3041944/1‑150 (MQ=60)
                                                                                                                                                    GAATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAG   >  SRR19973070.1090138/1‑150 (MQ=60)
                                                                                                                                                     AATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGC  >  SRR19973070.655105/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
GGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTGCCCGCCTGAACTGAACACTTCTGGCGGTACGTTTACCTCCACGTCCGGACGATAATGCGGGTTAGC  >  CP009273/1836043‑1836341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑