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
A0 F0 I1 R1 336 367.5 14491858 96.8% 14028118 150.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,837,678 C→T intergenic (+229/‑322) RPPX_17095 → / → RPPX_17100 proteophosphoglycan precursor/cupin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,837,6780CT91.2% 174.7 / 6.9 115intergenic (+229/‑322)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base T (45/59);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGT  >  CP009974/3837528‑3837793
                                                                                                                                                      |                                                                                                                    
gtTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGGGGCGGCACACACGGCGGTAGGAGCGCGGGCACCCCCGAAAAGACCCGGCCAGCGGGCCAATTTCGCCCCGCGGTTCACCCCCgccgccgcc                                                                                                                      >  2:3874511/1‑151 (MQ=9)
   cACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGGCTGCGCGGTCCTCTGTAGGAGCGGGTGCACCCGCGCAAGGGCCGGGCCGGCCCGCCAATTGTCACCAGCCGGTCACCGGTGCGGCCGCCAccctc                                                                                                                >  2:4534666/1‑148 (MQ=9)
   cACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGGTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCCGCCGCCCACCGTAGCTGCCGCCACCGTa                                                                                                                >  2:667946/1‑151 (MQ=9)
   cACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTAcc                                                                                                                   <  1:6081616/151‑1 (MQ=255)
    aCCGTTGGCGCCGCTACCGTATGCGCGGGTAAACCCGCTACTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCg                                                                                                                  <  2:592369/151‑1 (MQ=255)
    aCCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCg                                                                                                                  <  1:612182/151‑1 (MQ=255)
     ccATTGCCGCCGCTACTGTATTCGCGGGTCAACCCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                 <  1:2724719/151‑1 (MQ=255)
     ccATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGt                                                                                                                 <  1:2724447/151‑1 (MQ=255)
       aTTGCCGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTACACACGCGTCGAGGGGGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAt                                                                                                               <  1:4826913/151‑1 (MQ=255)
       aTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGGGGCCGCACCGCCCGCTGTAGGAGCGGGGCCACCCCCGACCGGGCCGCGCCCCCCGGCCAATTACCACCCGCCGGACACCCGTGCCGCCGCTACAGtct                                                                                                               >  2:6996364/1‑149 (MQ=2)
       aTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGCTGCCGCCGCCACCGtct                                                                                                               >  2:6380148/1‑149 (MQ=255)
        gtGCCGCCGCTACCGGATTCGCGGCCAAACCCGCCCCTACACCCGCGGCGAGGGCGGCACAGGCCGTTGTCGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAtt                                                                                                              <  2:3916097/150‑1 (MQ=11)
         tGCCGCCGCGACCGTACTCGCGTGTAAACCCGCTCCTACACACGCATCGAGTGCGGCAGAGGCCCTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTc                                                                                                             <  2:5785543/151‑1 (MQ=11)
         ggccgccGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTc                                                                                                             <  1:34977/150‑1 (MQ=37)
         cgccgccGCTAGCGTATTTGCGGTCAGACCCGCTGCTACGCTCGCAGCGTGGGCGGCACAGGCCGTTGTAGGGGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTc                                                                                                             <  2:1564502/150‑1 (MQ=11)
          gccgccGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGATAACGTATTcg                                                                                                            >  1:2212811/1‑151 (MQ=255)
           ccgccgCTGCCGTCGTGGCGGGCAAACCCGCTCCTGCACACGCGTCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcgc                                                                                                           <  2:2002286/151‑1 (MQ=11)
           ccgccgCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGCCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTcgc                                                                                                           >  2:145255/1‑151 (MQ=255)
           ccgccgCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCGGGCCAATTTCGGCCCGCGGTTCACCCTTGCCCCCGCCACCGTATTcgc                                                                                                           >  2:5429941/1‑151 (MQ=14)
            cgccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGGCGCTACCGTATTtgcg                                                                                                          >  1:4948163/1‑151 (MQ=255)
            cgccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCCCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCCCCGTTGCCGCCGCTCCCCTATTcgcc                                                                                                          >  2:3825879/1‑150 (MQ=11)
             gccgcTACGGTACTCGCGGGTAAACCCGCTCCTACACCCGCGTCGAGGGCGGCACAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                         <  2:6240240/151‑1 (MQ=11)
             gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                         <  1:6109592/151‑1 (MQ=38)
             gccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAACTTCCACCCGCCGTTCACCGTTGCCGCCCCTACCGTATCcgccc                                                                                                         >  2:1868589/1‑149 (MQ=255)
             gccgcGACCGGCTTCGCGGGTAGACCGGCCCCTACAGCCGCGTCGAGGGCGGCACAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCgg                                                                                                         <  2:4919572/151‑1 (MQ=11)
              gcgcgACCGTGTTCGCGGGTCACCCCGCGCCTACGCACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCggg                                                                                                        <  2:5619905/150‑1 (MQ=11)
              gcgcTACGGTGTTCGCGGGTAAACCGGCTCCTACACACGCGTCGAGGGCGGCAGAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCggg                                                                                                        <  2:1985851/150‑1 (MQ=11)
              ccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCggg                                                                                                        <  1:665868/151‑1 (MQ=37)
              ccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCggg                                                                                                        >  2:419778/1‑151 (MQ=37)
              ccgcTACCGTATTCGCGGGTAAACCCGCTCCTACACTGCTCTGGCCTGCGCGGTCCTCTGTCGCAGCGGGTTCACCCCCGAAAGCGCCGGGCCTGCTCGCCCATTTTCACCAGCCGCTCCCCGTTGCTGCCGCCCCCGTATTCCCGGGCaa                                                                                                     >  2:1690578/1‑151 (MQ=9)
              ccgcTACCGTATTCGCGGGTAAACCCGCTCCGCCTCCCGCATCGAGGGCGGCCCAGGCCGGGGTAGGAGCGGGTTCCCCCGCGAAAGGGCCGGGCCTGCCGGGCAATTTCCACCCGCCGTCCACCGTTGGCGCCGCGGCCGTATGTGCggg                                                                                                        >  2:1662855/1‑151 (MQ=14)
               cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                       >  1:3150660/1‑151 (MQ=38)
               cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                       >  1:3437865/1‑151 (MQ=38)
               cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGt                                                                                                       >  1:3739869/1‑151 (MQ=38)
               cgcTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTCGAAGGAGCGGGTTCACCCGCGAGAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGCTCCCCGTTGCCGCCGCTACCGTCTTCGCGGGt                                                                                                       >  1:3739848/1‑151 (MQ=11)
                 cTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTaa                                                                                                     >  2:352510/1‑151 (MQ=25)
                 cTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGCaa                                                                                                     >  1:2066109/1‑151 (MQ=25)
                     cGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAccc                                                                                                 >  1:6029902/1‑151 (MQ=25)
                     cGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAccc                                                                                                 <  1:1134982/151‑1 (MQ=25)
                      gTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                >  1:2264430/1‑151 (MQ=21)
                      gTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                >  1:2048370/1‑151 (MQ=21)
                      gTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCg                                                                                                <  1:3808026/151‑1 (MQ=25)
                          tCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                            <  1:6691975/151‑1 (MQ=18)
                          tCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                            >  1:703470/1‑151 (MQ=25)
                          gcgcgGGTAGACCCGCTGCTACACCCGGAGCGAGGGCGGCACAGGCCGTTGTAGGGGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                            <  2:5782163/150‑1 (MQ=11)
                          gcgcgGGTAGACCCGCGCCTACGCGCGCGGCGGGGGCGGCAGAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                            <  2:2930639/150‑1 (MQ=11)
                          ccgcgGGTAGACCCGATCGTACACAGGAATCGAAGGCGGAGCAGGCCGTTGTCGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTcc                                                                                            <  1:5166651/150‑1 (MQ=11)
                           ggcgGGTAAACCCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCt                                                                                           <  1:6354057/150‑1 (MQ=25)
                           cgcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAACCCCGccccc                                                                                           >  2:7228749/1‑150 (MQ=255)
                           cgcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTccc                                                                                           >  2:3277168/1‑150 (MQ=25)
                           cgcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGGAAACCCGCTCCt                                                                                           >  2:6332381/1‑151 (MQ=21)
                            gtgggTAGACCCGCTGCGACACACGCGTCGGGGGCGGCACAGGCCGTTGTCGGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtg                                                                                          <  2:5074338/149‑2 (MQ=11)
                            gcgGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGAGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTGGCCGCGGCTACCGTATTCGCGGGTAAACGCGCTCCtg                                                                                          >  2:5664017/1‑150 (MQ=11)
                             ggggTAAACCCGCTCCTGCGCCGGCATCGTGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCtgc                                                                                         <  1:4941941/150‑3 (MQ=11)
                              gggTCAACCCGCCCCTACACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgca                                                                                        <  1:6469116/151‑1 (MQ=11)
                                gTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTCAACCCGCTcccgcacc                                                                                      >  2:6771479/1‑145 (MQ=255)
                                 tAAACCGGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcacac                                                                                     <  1:7122654/151‑1 (MQ=255)
                                 tAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcacac                                                                                     >  2:6371442/1‑151 (MQ=18)
                                 tAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTgcacac                                                                                     >  1:1069935/1‑151 (MQ=17)
                                   agcccGCACCTACAAACGCATCAAGGGCGGCAAAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGc                                                                                   <  1:997845/149‑1 (MQ=11)
                                   aaCCCGCTCCTGCGCACGCGGCGAGGGCGGCACAGGACGTTGTAGGAGCGGGCTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGc                                                                                   <  2:6379938/151‑1 (MQ=11)
                                       cGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTACACCCGCTCCTGCCCACGCATcc                                                                               >  2:3661188/1‑150 (MQ=255)
                                       cGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCg                                                                               >  2:3943138/1‑151 (MQ=12)
                                       cGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCg                                                                               >  1:5966978/1‑151 (MQ=12)
                                       cGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCACCg                                                                               >  2:3100928/1‑151 (MQ=12)
                                          tcatACACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCggtg                                                                            <  1:2709785/148‑4 (MQ=11)
                                             tACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTCCACACGCATCggtggtt                                                                         >  2:3999330/1‑145 (MQ=9)
                                                                                       tgcACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTGCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAgggccgggc                           <  2:4541906/149‑1 (MQ=2)
                                                                                          ccccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATCCGCGGGTAACCCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctg                        <  1:1533915/150‑1 (MQ=2)
                                                                                           cccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgc                       <  1:954749/151‑1 (MQ=2)
                                                                                            ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                      <  2:2466324/151‑1 (MQ=11)
                                                                                            ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                      <  2:6350226/151‑1 (MQ=2)
                                                                                            ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                      <  1:5036672/151‑1 (MQ=2)
                                                                                            ccGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTCCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgct                      <  2:2319301/151‑1 (MQ=2)
                                                                                             cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                     <  1:4872579/151‑1 (MQ=2)
                                                                                             cgcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCATACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                     <  2:310282/151‑1 (MQ=2)
                                                                                              gcgAAAGGGCCGGGCCTGCTGGCCAGTGTTCACCGGCCGCGCACCGTTGCTGCCGCCACCGGATTCGCGGGTAAACCCGCTCCGACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTgg                    <  2:4764041/151‑1 (MQ=2)
                                                                                              gcgAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTgg                    <  2:422535/151‑1 (MQ=2)
                                                                                                gAAAGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGcc                  >  1:5310161/1‑151 (MQ=2)
                                                                                                 aaaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCa                 >  1:6613035/1‑151 (MQ=2)
                                                                                                 aaaGGGCCGGGCCTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCa                 <  1:4193960/151‑1 (MQ=2)
                                                                                                    gggccgggccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtt              >  1:3268388/1‑150 (MQ=2)
                                                                                                     ggccgggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATGCGCGGGTAACCCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttt             <  2:1822290/151‑3 (MQ=2)
                                                                                                      gcggggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAACCCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtttt            <  1:5417885/148‑4 (MQ=2)
                                                                                                      gccgggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtttt            <  2:4483641/151‑4 (MQ=2)
                                                                                                       cggggcgTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACCCTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAttttc           <  2:6737841/149‑5 (MQ=2)
                                                                                                        cgggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCGTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCACTCTTCa          >  1:4679083/1‑151 (MQ=2)
                                                                                                         gggccTGCTGGCCAGTTTTCACCGGCCGTTCACCGTTGCTGCCGCTACCGTATCCGCGGGTGGCGCCGCGCCTACACTGCTCGGGAGTGCGCGGGCCGCTGTAGGGGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAc         <  2:3039169/151‑1 (MQ=2)
                                                                                                          ggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcc        >  1:6179779/1‑151 (MQ=2)
                                                                                                          ggccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcc        <  2:3395752/151‑1 (MQ=2)
                                                                                                           gccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccc       >  1:5171361/1‑151 (MQ=2)
                                                                                                            ccTGCTGGCCAATTTTCACCAGCCGTGCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccag      <  2:3813885/151‑3 (MQ=2)
                                                                                                            ccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccag      >  1:3541201/1‑149 (MQ=2)
                                                                                                              tgctgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGCTCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCCATTTTCAccagcc    >  1:2231078/1‑147 (MQ=2)
                                                                                                              tgctgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagcc    >  1:3376517/1‑147 (MQ=2)
                                                                                                               gctgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccg   >  1:3731712/1‑146 (MQ=2)
                                                                                                                ctgGCCAATGTCCACCAGCCGTCCACCGTTGCTGCCGCCCCCGTAGTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccgc  <  1:5190233/151‑7 (MQ=2)
                                                                                                                   gCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtttt            >  2:1051431/1‑135 (MQ=2)
                                                                                                                   gCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtttt            <  1:1051431/138‑4 (MQ=2)
                                                                                                                                               ctgccgccACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTgg                    <  1:4090654/100‑1 (MQ=2)
                                                                                                                                                      |                                                                                                                    
GTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGT  >  CP009974/3837528‑3837793

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGC  >  CP009974/3837528‑3837722
                                                                                                                                                      |                                            
GTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCT                                              >  SRR10918235.6259752/1‑151 (MQ=51)
   CACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACC                                           <  SRR10918235.6109804/151‑1 (MQ=44)
    ACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCG                                          <  SRR10918235.614397/151‑1 (MQ=40)
     CCATTGCCGCCGCTACTGTATTCGCGGGTCAACCCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGT                                         <  SRR10918235.2736436/151‑1 (MQ=35)
     CCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGT                                         <  SRR10918235.2736164/151‑1 (MQ=38)
     CCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCAGGGGCGGCACAGGCCCTTGTAGAAGCCAGTTCACCCGCGCAAGGGCCGCGGCTGCAGGCCAATTTCAACCCCCAGTTCACCCTTGCCGCCCCTACCGT                                         >  SRR10918235.1092389/1‑151 (MQ=21)
       ATTGCCGCCGCTCCCGTATTCGCGGGTAAACCCGCTCCTACACACGCGTCGAGGGGGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAT                                       <  SRR10918235.4848279/151‑1 (MQ=31)
         GGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTC                                     <  SRR10918235.35101/151‑1 (MQ=31)
          GCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGATAACGTATTCG                                    >  SRR10918235.2222723/1‑151 (MQ=36)
            CGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGGCGCTACCGTATTTGCG                                  >  SRR10918235.4970077/1‑151 (MQ=30)
             GCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGG                                 <  SRR10918235.6137919/151‑1 (MQ=31)
              CCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGG                                <  SRR10918235.668272/151‑1 (MQ=31)
               CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                               >  SRR10918235.3164365/1‑151 (MQ=29)
               CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                               >  SRR10918235.3453099/1‑151 (MQ=29)
               CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGT                               >  SRR10918235.3756514/1‑151 (MQ=29)
               CGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTCGAAGGAGCGGGTTCACCCGCGAGAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGCTCCCCGTTGCCGCCGCTACCGTCTTCGCGGGT                               >  SRR10918235.3756493/1‑151 (MQ=25)
                 CTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGCAA                             >  SRR10918235.2075397/1‑151 (MQ=25)
                   ACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAAC                           >  SRR10918235.2826971/1‑151 (MQ=21)
                     CGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCC                         >  SRR10918235.6057834/1‑151 (MQ=21)
                     CGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCC                         >  SRR10918235.4655807/1‑151 (MQ=21)
                     CGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCC                         <  SRR10918235.1139449/151‑1 (MQ=21)
                     CGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGCTAAACCC                         >  SRR10918235.865363/1‑151 (MQ=20)
                      GTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCG                        <  SRR10918235.3824939/151‑1 (MQ=21)
                      GTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCG                        >  SRR10918235.2057598/1‑151 (MQ=21)
                      GTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCG                        >  SRR10918235.2274546/1‑151 (MQ=21)
                       TATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGCGTAAACCCGCTCCTACACAC             >  SRR10918235.2337182/1‑151 (MQ=0)
                       TATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGGTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACAC             >  SRR10918235.2334289/1‑151 (MQ=0)
                          TCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCC                    >  SRR10918235.705975/1‑151 (MQ=21)
                          TCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCC                    <  SRR10918235.6723345/151‑1 (MQ=21)
                          CCGCGGGTAGACCCGATCGTACACAGGAATCGAAGGCGGAGCAGGCCGTTGTCGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCC                    <  SRR10918235.5189753/151‑1 (MQ=16)
                           GGCGGGTAAACCCGCTCCTACACCCGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCT                   <  SRR10918235.6383607/151‑1 (MQ=20)
                             GGGGTAAACCCGCTCCTGCGCCGGCATCGTGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGC                 <  SRR10918235.4963822/151‑1 (MQ=14)
                             CGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACatgcac       >  SRR10918235.2541579/1‑145 (MQ=0)
                              gagtaaacccgggcctacgagcgcaccgagggcggcacgggtcgtggtaggtgcggttccgCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCA                <  SRR10918235.6339151/90‑1 (MQ=0)
                              GGGTCAACCCGCCCCTACACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCA                <  SRR10918235.6499273/151‑1 (MQ=10)
                              GGGTAAACCCCGTCCTACACGCGCATCGAGGGCGGCTCAGGCCGATGTGGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCA                <  SRR10918235.6340099/151‑1 (MQ=9)
                               gccaaaccggtccctacagctgcagcgagtgctgcccatgccgtcctatGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCAC               <  SRR10918235.1622246/102‑1 (MQ=0)
                               GGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACatgcacag     >  SRR10918235.518849/1‑143 (MQ=0)
                               GGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTCCTACACACatgcacag     >  SRR10918235.5888621/1‑143 (MQ=0)
                                 TAAACCGGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACAC             <  SRR10918235.7156739/151‑1 (MQ=7)
                                 TAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACAC             >  SRR10918235.1074074/1‑151 (MQ=7)
                                   AGCCCGCACCTACAAACGCATCAAGGGCGGCAAAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGC           <  SRR10918235.1001600/151‑1 (MQ=6)
                                   AACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGC           >  SRR10918235.6020673/1‑151 (MQ=7)
                                      ccgccccagcacacgaggcgaaGGCGGCACAGGCCGTTGTCGGCGCGGGTAAACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATC        <  SRR10918235.851075/129‑1 (MQ=6)
                                      CCCCTCCTACGCACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATC        <  SRR10918235.5121724/151‑1 (MQ=7)
                                       CGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCG       >  SRR10918235.5994609/1‑151 (MQ=7)
                                          TCATACACACGCAGCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtg    <  SRR10918235.2721451/151‑4 (MQ=1)
                                          GCCTACAGACGCATCGGGGGCGGCACGGGCCGTTGTCGGGGCGGGTTCACCTGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtg    <  SRR10918235.4645106/151‑4 (MQ=0)
                                           CCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtgg   <  SRR10918235.3882132/151‑5 (MQ=0)
                                           CCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtgg   <  SRR10918235.7120058/151‑5 (MQ=0)
                                           CCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtgg   >  SRR10918235.2457310/1‑147 (MQ=0)
                                            CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGgtggt  >  SRR10918235.3681891/1‑146 (MQ=0)
                                                                                                                                                      |                                            
GTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGC  >  CP009974/3837528‑3837722

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