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
A2 F1 I1 R1 404 265.1 10478538 96.8% 10143224 150.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,837,790 T→C intergenic (+341/‑210) 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,7900TC97.5% 356.4 / ‑4.4 120intergenic (+341/‑210)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (67/48);  total (68/52)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 121 total aligned reads displayed.

AATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAT  >  CP009974/3837646‑3837937
                                                                                                                                                |                                                                                                                                                   
aaTTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGCTCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCGGCCGGCCAATCTTCACCAGCCGCTCAccg                                                                                                                                            >  2:2348121/1‑150 (MQ=2)
                                                          cTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTCTTCCCGGGTAAACCCGCTCCAGCACACGCATCggtggt                                                                                     >  1:1326978/1‑151 (MQ=11)
                                                             cacacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttgg                                                                                  >  1:5047004/1‑151 (MQ=9)
                                                              acacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                                 >  2:4655247/1‑151 (MQ=255)
                                                              acacGCACCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                                 <  1:1541744/151‑1 (MQ=255)
                                                               cacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCTCACGCAGCGGGGGTTGGTc                                                                                >  1:2670419/1‑151 (MQ=11)
                                                                acGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCCCACGCATCGGTGGTTGGTca                                                                               >  2:4849748/1‑151 (MQ=255)
                                                                  gcggccggggCGGCCGAGGCCGTCGTCGGCGCGGGTCCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                             <  2:5218817/144‑1 (MQ=11)
                                                                  gCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                             >  2:612151/1‑151 (MQ=255)
                                                                   cagcgtgggCGGCACAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAg                                                                            <  1:1603120/145‑1 (MQ=11)
                                                                   cATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAg                                                                            >  2:1404092/1‑151 (MQ=255)
                                                                    aTCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGGGGTTGGTCACAgg                                                                           >  1:4154472/1‑151 (MQ=11)
                                                                       gAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcg                                                                        <  1:2746024/151‑1 (MQ=255)
                                                                       gAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcg                                                                        >  1:2729628/1‑151 (MQ=255)
                                                                        aGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgt                                                                       >  2:4352168/1‑151 (MQ=255)
                                                                         gggCGGCGCAGGCCGTGGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCcgtc                                                                      <  2:4955671/151‑1 (MQ=14)
                                                                            cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                   <  1:4766392/151‑1 (MQ=255)
                                                                            cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                                   >  1:1106779/1‑151 (MQ=255)
                                                                             ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCCCGCTCACGCCCCGGTGGCTGGGCACAGGCCGTCGCAg                                                                  >  2:195752/1‑151 (MQ=11)
                                                                              gcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAgg                                                                 >  2:842508/1‑151 (MQ=255)
                                                                               cacaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGa                                                                >  1:3362563/1‑151 (MQ=255)
                                                                               cacaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGa                                                                <  2:1106779/151‑1 (MQ=255)
                                                                                acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAg                                                               <  1:2329531/151‑1 (MQ=255)
                                                                                acaGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGCAg                                                               >  2:1561242/1‑151 (MQ=255)
                                                                                  aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                             >  1:2060636/1‑151 (MQ=255)
                                                                                   ggCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                            <  2:1981728/151‑1 (MQ=14)
                                                                                   ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                            <  2:3095580/151‑1 (MQ=255)
                                                                                   ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                            >  2:415579/1‑151 (MQ=255)
                                                                                      cGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGtt                                                         >  2:1662721/1‑151 (MQ=255)
                                                                                       gTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTc                                                        >  1:82339/1‑151 (MQ=255)
                                                                                          ttAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                                     <  1:4482402/150‑1 (MQ=255)
                                                                                           tAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAccc                                                    <  1:3276144/151‑1 (MQ=255)
                                                                                            aGGAGCGGGTACACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcg                                                   <  2:332421/151‑1 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  >  1:1933904/1‑151 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  >  1:1734243/1‑151 (MQ=255)
                                                                                             ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  >  1:1933816/1‑151 (MQ=255)
                                                                                             ggAGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  <  1:4484701/151‑1 (MQ=255)
                                                                                               aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                >  2:4539898/1‑151 (MQ=255)
                                                                                               aGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGa                                                <  2:3737219/151‑1 (MQ=255)
                                                                                                gCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGCGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCCaa                                               >  2:2661177/1‑151 (MQ=14)
                                                                                                 cGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaa                                              >  1:2879766/1‑151 (MQ=255)
                                                                                                  gggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGaaaa                                             >  1:1727971/1‑151 (MQ=255)
                                                                                                   ggtaaaCCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                            <  2:1462053/146‑1 (MQ=14)
                                                                                                   ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                            >  1:3132835/1‑151 (MQ=255)
                                                                                                   ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                            <  1:1309753/151‑1 (MQ=255)
                                                                                                   ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                            >  1:1159365/1‑151 (MQ=255)
                                                                                                    gtggagccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGACGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAgg                                           <  2:2263540/145‑1 (MQ=14)
                                                                                                         cccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGt                                      >  1:3900569/1‑151 (MQ=255)
                                                                                                           cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGCGCGGGTTCACCCGCGAACAGGCCGGTAc                                    >  2:4709308/1‑151 (MQ=14)
                                                                                                           cgcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAc                                    <  1:1175574/151‑1 (MQ=255)
                                                                                                            gcgAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACa                                   >  1:782524/1‑151 (MQ=255)
                                                                                                              gAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAgg                                 <  1:3086456/151‑1 (MQ=255)
                                                                                                              gAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  <  2:515672/134‑1 (MQ=255)
                                                                                                              gAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                                  >  1:515672/1‑134 (MQ=255)
                                                                                                               caaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                <  1:3507664/150‑1 (MQ=255)
                                                                                                               aaaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                <  1:5003643/151‑1 (MQ=255)
                                                                                                               aaaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGc                                <  1:4418896/151‑1 (MQ=255)
                                                                                                                aaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCg                               >  2:659539/1‑151 (MQ=255)
                                                                                                                aaGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCg                               <  1:4216980/151‑1 (MQ=255)
                                                                                                                 aGGGCCGGGCCTGCTGGCCAATTGCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGGGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGa                              <  2:2306565/151‑1 (MQ=255)
                                                                                                                   ggccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGCAAAGGCCGGTACAGGCGaaa                            >  2:4558870/1‑151 (MQ=255)
                                                                                                                   ggccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGCTTCACCCGCGAAAAGGCCGGTACAGGCGaaa                            >  2:4058156/1‑151 (MQ=255)
                                                                                                                     ccgggccTGCTGTCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGTCAACCCCGCCCCTGCACACGCATCGCTCGTCGGTCACAGGCCGTCGTAGGACCTGGTTCACCCGCGAAACGGCCCGTCCCAGCGAAAAt                          >  2:3755987/1‑151 (MQ=14)
                                                                                                                     ccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGACGAAAAt                          >  2:3218224/1‑151 (MQ=255)
                                                                                                                      cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGCAAATc                         >  2:3334103/1‑151 (MQ=255)
                                                                                                                      cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                         >  1:2082343/1‑151 (MQ=255)
                                                                                                                       gggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAACTCg                        >  2:3241373/1‑151 (MQ=255)
                                                                                                                       gggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                        <  1:2226269/151‑1 (MQ=255)
                                                                                                                       gggccTGCGGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                        <  2:302803/151‑1 (MQ=255)
                                                                                                                       gcgccTGCTGGCCAATTTCCACCCTCGGTTCACCGTCGCCGCCGCTACCGTATTCGCGGGTGAACCCGCGCCTGCACACGCATCGGTGGTGGGTCACAGGCCGTGGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCg                        <  1:2226318/149‑1 (MQ=17)
                                                                                                                        ggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCgg                       >  2:3653819/1‑151 (MQ=255)
                                                                                                                         gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                      >  1:2711800/1‑151 (MQ=255)
                                                                                                                         gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGCTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc                      >  2:4013095/1‑151 (MQ=255)
                                                                                                                          ccTGCGGGCCAATTTCCCCCCGCGGCTCACCGTTGCGGCGGCGACCGTGTTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTGGGGCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc                     <  2:1396955/151‑1 (MQ=14)
                                                                                                                           ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                    >  1:4086133/1‑151 (MQ=255)
                                                                                                                           ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc                    >  1:892767/1‑151 (MQ=255)
                                                                                                                            tgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcccc                   >  2:1349565/1‑151 (MQ=255)
                                                                                                                             gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGTGGTCGGTCACAGGCCGTCGTAGGAGCCGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAACCGGccccc                  >  2:4365754/1‑151 (MQ=17)
                                                                                                                             gctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccccc                  <  1:1869288/151‑1 (MQ=255)
                                                                                                                              ctgGCCACTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCt                 <  2:1518846/151‑1 (MQ=255)
                                                                                                                               tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                >  2:3407229/1‑151 (MQ=255)
                                                                                                                               tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg                <  1:821614/151‑1 (MQ=255)
                                                                                                                                ggCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGGTACCGTATTGGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGGGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc               <  1:1949510/151‑1 (MQ=18)
                                                                                                                                ggCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc               >  2:2110695/1‑151 (MQ=255)
                                                                                                                                ggCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGc               <  1:4872419/151‑1 (MQ=255)
                                                                                                                                  ccAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtt             >  2:4689175/1‑151 (MQ=255)
                                                                                                                                   cAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCttt            <  2:4875736/151‑1 (MQ=255)
                                                                                                                                    aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGTTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt           <  1:2899447/151‑1 (MQ=255)
                                                                                                                                    aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt           <  2:724072/151‑1 (MQ=255)
                                                                                                                                    aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt           >  1:3090902/1‑151 (MQ=255)
                                                                                                                                    aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCTTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGACAATCGGCCCCCTGCtttt           >  1:1639820/1‑151 (MQ=18)
                                                                                                                                     atttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTc          <  2:2729628/148‑1 (MQ=255)
                                                                                                                                      tttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcc         >  1:4143456/3‑151 (MQ=255)
                                                                                                                                      tttcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcc         <  1:2367061/149‑1 (MQ=255)
                                                                                                                                        tcCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTcccc       >  2:3890398/1‑151 (MQ=255)
                                                                                                                                          cACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCg     <  1:528771/151‑1 (MQ=255)
                                                                                                                                           ccccGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg    <  1:3588564/150‑1 (MQ=255)
                                                                                                                                           aCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg    >  2:1331891/1‑151 (MQ=255)
                                                                                                                                           aCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCgg    >  1:321053/1‑151 (MQ=255)
                                                                                                                                             ccgccgTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGGACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAt  >  1:903283/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAT  >  CP009974/3837646‑3837937

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAT  >  CP009974/3837675‑3837937
                                                                                                                   |                                                                                                                                                   
CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACA‑‑CAT‑‑‑‑‑‑‑‑GCACAGGCCGTCGTAGGAGCGGGTTTACCCGCGAAGAGGACGGGCCTGCTGGCCAATCTCCACCCGCCGTTCACCGTTGCCGCCACCACCGTATTCGCGGGTAAACCCGCTC                                                                                                        >  SRR10918231.2739685/1‑151 (MQ=0)
                             CTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTCTTCCCGGGTAAACCCGCTCCAGCACACGCATCGGTGGT                                                                                     >  SRR10918231.1333338/1‑151 (MQ=0)
                                CACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGG                                                                                  >  SRR10918231.5071568/1‑151 (MQ=3)
                                 ACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGT                                                                                 <  SRR10918231.1547602/151‑1 (MQ=5)
                                 ACACGCACCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGT                                                                                 <  SRR10918231.1548848/151‑1 (MQ=5)
                                     GCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACA                                                                             <  SRR10918231.2676468/151‑1 (MQ=16)
                                      CAGCGTGGGCGGCACAGGCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAG                                                                            <  SRR10918231.1610511/151‑1 (MQ=17)
                                       ATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCCTCGGGGGTTGGTCACAGG                                                                           >  SRR10918231.4174100/1‑151 (MQ=5)
                                         cctaggtgagggaggccggcctcttcgccggtcaACCCGCGAAAGGGCCAGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCC                                                                         <  SRR10918231.3231841/117‑1 (MQ=25)
                                          GAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCG                                                                        >  SRR10918231.2742172/1‑151 (MQ=26)
                                          GAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCG                                                                        <  SRR10918231.2758630/151‑1 (MQ=26)
                                               CGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTA                                                                   >  SRR10918231.1111855/1‑151 (MQ=36)
                                               CGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTA                                                                   <  SRR10918231.4789214/151‑1 (MQ=36)
                                                  CACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGA                                                                >  SRR10918231.3377805/1‑151 (MQ=42)
                                                   ACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAG                                                               <  SRR10918231.2340090/151‑1 (MQ=44)
                                                     AGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCG                                                             >  SRR10918231.2070065/1‑151 (MQ=48)
                                                     AGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCG                                                             >  SRR10918231.2069862/1‑151 (MQ=48)
                                                      ggccggtctcttcgcgggtaaACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGG                                                            <  SRR10918231.177574/130‑1 (MQ=52)
                                                         CGTGGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTT                                                         <  SRR10918231.3695432/151‑1 (MQ=54)
                                                          gtcgtattggcaggtaaACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCTTCGTAGGAGCGGGTTC                                                        <  SRR10918231.1316834/134‑1 (MQ=48)
                                                          GTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTC                                                        >  SRR10918231.82639/1‑151 (MQ=58)
                                                             TTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACC                                                     <  SRR10918231.4503711/151‑1 (MQ=59)
                                                              TAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCC                                                    <  SRR10918231.3291009/151‑1 (MQ=60)
                                                                GGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                                  >  SRR10918231.1942609/1‑151 (MQ=60)
                                                                GGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                                  >  SRR10918231.1942521/1‑151 (MQ=60)
                                                                GGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                                  >  SRR10918231.1742211/1‑151 (MQ=60)
                                                                GGAGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                                  <  SRR10918231.4506024/151‑1 (MQ=58)
                                                                    CGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAA                                              >  SRR10918231.2892927/1‑151 (MQ=60)
                                                                     GGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAA                                             >  SRR10918231.1735900/1‑151 (MQ=60)
                                                                      GGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAG                                            >  SRR10918231.1164756/1‑151 (MQ=60)
                                                                      GGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAG                                            >  SRR10918231.3147140/1‑151 (MQ=60)
                                                                      GGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAG                                            <  SRR10918231.1316053/151‑1 (MQ=60)
                                                                       gtttgcccgcgaagaGGCCGGCCCTTCTGGCCAATTACCACCCGTCGTTCACCATTGCCGCTGCTACCTTATTCGCGGGcaacaccgcccccccgccgcactggggcttgacacccggcggtgacgctcggacgctcgactgtctgaactg                                           >  SRR10918231.1614556/16‑79 (MQ=0)
                                                                       ctcCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGG                                           <  SRR10918231.5196778/148‑1 (MQ=60)
                                                                        TTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGC                                          >  SRR10918231.1869356/1‑151 (MQ=60)
                                                                          CACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCGGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCG                                        <  SRR10918231.3564603/151‑1 (MQ=60)
                                                                            CCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGT                                      >  SRR10918231.3918731/1‑151 (MQ=60)
                                                                              CGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTAC                                    <  SRR10918231.1181067/151‑1 (MQ=60)
                                                                               GCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACA                                   >  SRR10918231.785682/1‑151 (MQ=60)
                                                                                 GAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGG                                 <  SRR10918231.3100547/151‑1 (MQ=60)
                                                                                 GAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCagatcggaagagcacac                                 >  SRR10918231.517497/1‑134 (MQ=60)
                                                                                  CAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGC                                <  SRR10918231.3523644/151‑1 (MQ=60)
                                                                                  AAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGC                                <  SRR10918231.5027898/151‑1 (MQ=60)
                                                                                  AAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGC                                <  SRR10918231.4439853/151‑1 (MQ=60)
                                                                                   AAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCG                               <  SRR10918231.4236893/151‑1 (MQ=60)
                                                                                       GCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAA                           >  SRR10918231.5314/1‑151 (MQ=60)
                                                                                         CGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATC                         >  SRR10918231.2091678/1‑151 (MQ=60)
                                                                                         CGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATC                         <  SRR10918231.4370532/151‑1 (MQ=60)
                                                                                          GGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCG                        <  SRR10918231.2236259/151‑1 (MQ=60)
                                                                                          GCGCCTGCTGGCCAATTTCCACCCTCGGTTCACCGTCGCCGCCGCTACCGTATTCGCGGGTGAACCCGCGCCTGCACACGCATCGGTGGTGGGTCACAGGCCGTGGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCG                        <  SRR10918231.2236308/151‑1 (MQ=60)
                                                                                            GCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGC                      <  SRR10918231.2840408/151‑1 (MQ=60)
                                                                                            GCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGC                      >  SRR10918231.2724277/1‑151 (MQ=60)
                                                                                             CCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCC                     <  SRR10918231.321087/151‑1 (MQ=60)
                                                                                              CTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCC                    >  SRR10918231.4105360/1‑151 (MQ=60)
                                                                                              CTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCC                    >  SRR10918231.896540/1‑151 (MQ=60)
                                                                                               TGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCC                   >  SRR10918231.4889401/1‑151 (MQ=60)
                                                                                                GCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCC                  <  SRR10918231.1877703/151‑1 (MQ=60)
                                                                                                  TGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTG                <  SRR10918231.824981/151‑1 (MQ=60)
                                                                                                   GGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGGTACCGTATTGGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGGGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGC               <  SRR10918231.1958297/151‑1 (MQ=60)
                                                                                                   GGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGC               <  SRR10918231.4895884/151‑1 (MQ=60)
                                                                                                       AATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGTTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTT           <  SRR10918231.2912702/151‑1 (MQ=60)
                                                                                                       AATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTT           >  SRR10918231.3105013/1‑151 (MQ=60)
                                                                                                       AATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCTTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGACAATCGGCCCCCTGCTTTT           >  SRR10918231.1647398/1‑151 (MQ=60)
                                                                                                         TTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCC         <  SRR10918231.2377832/151‑1 (MQ=60)
                                                                                                         TTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCC         >  SRR10918231.4163027/1‑151 (MQ=60)
                                                                                                             CACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCG     <  SRR10918231.530666/151‑1 (MQ=60)
                                                                                                              CCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGG    <  SRR10918231.3604937/151‑1 (MQ=60)
                                                                                                              ACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGG    >  SRR10918231.322206/1‑151 (MQ=60)
                                                                                                                CCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGGACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAT  >  SRR10918231.907131/1‑151 (MQ=60)
                                                                                                                   |                                                                                                                                                   
CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCCCCCGGAT  >  CP009974/3837675‑3837937

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