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,781 C→T intergenic (+332/‑219) 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,7810CT95.7% 350.0 / ‑2.8 115intergenic (+332/‑219)RPPX_17095/RPPX_17100proteophosphoglycan precursor/cupin
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (63/47);  total (66/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

GGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTT  >  CP009974/3837634‑3837928
                                                                                                                                                   |                                                                                                                                                   
ggccTGCTGGCCATTTTTCACCGGCCGTGCACCGTTGCTGCCGCTACCGTAGTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcc                                                                                                                                               <  2:3480885/151‑1 (MQ=2)
 gccTGCTGGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAcca                                                                                                                                              >  1:4726002/1‑150 (MQ=2)
 gccTGCTGGCCAATTTTCACCAGCCGCTCACCGTTGCTGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGGTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTTCAccc                                                                                                                                              >  2:1762978/1‑151 (MQ=9)
     gctgGCCAATTTTCACCAGCCGTTCACCGTTGCTGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACTGATCTGGACTGCGCGGTCCTCTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTTCAccagccg                                                                                                                                          >  1:4043833/1‑146 (MQ=2)
            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)
                                                                          acacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                        <  1:1540499/151‑1 (MQ=255)
                                                                          acacGCACCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGtggttggt                                                                        <  1:1541744/151‑1 (MQ=255)
                                                                           cacGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCTCACGCAGCGGGGGTTGGTc                                                                       >  1:2670419/1‑151 (MQ=11)
                                                                            acGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCCCACGCATCGGTGGTTGGTca                                                                      >  2:4849748/1‑151 (MQ=255)
                                                                             cGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACAGTTGCCGCCGCTACCGTATTCGCGGGTGAACCAGCTCCTGGACACGCATCGGTGGCTGGTcca                                                                     >  2:3136998/1‑149 (MQ=11)
                                                                              gcggccggggCGGCCGAGGCCGTCGTCGGCGCGGGTCCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                    <  2:5218817/144‑1 (MQ=11)
                                                                              gCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTcaca                                                                    >  2:612151/1‑151 (MQ=255)
                                                                               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)
                                                                                        cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                          <  1:4766392/151‑1 (MQ=255)
                                                                                        cggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTa                                                          >  1:1106779/1‑151 (MQ=255)
                                                                                         ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCCCACGACTCGGTGGTTGGTCCCAGGCcgtcgtcg                                                         >  2:3615442/1‑149 (MQ=14)
                                                                                         ggcACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCCCGCTCACGCCCCGGTGGCTGGGCACAGGCCGTCGCAg                                                         >  2:195752/1‑151 (MQ=11)
                                                                                           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)
                                                                                             caGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGc                                                     >  2:4474875/1‑151 (MQ=255)
                                                                                              aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                    >  1:2060838/1‑151 (MQ=255)
                                                                                              aGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCg                                                    >  1:2060636/1‑151 (MQ=255)
                                                                                               ggCCGTTGTAGGCGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                   <  2:1981728/151‑1 (MQ=14)
                                                                                               ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                   >  2:415579/1‑151 (MQ=255)
                                                                                               ggCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCgg                                                   <  2:3095580/151‑1 (MQ=255)
                                                                                                 ccGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGt                                                 >  2:1418241/1‑151 (MQ=255)
                                                                                                  cGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGtt                                                >  2:1662761/1‑151 (MQ=255)
                                                                                                  cGTGGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGtt                                                <  1:3678577/151‑1 (MQ=255)
                                                                                                   gTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTc                                               >  1:82339/1‑151 (MQ=255)
                                                                                                      gTAGGCGCGGGTTGACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAcc                                            <  2:1052279/151‑1 (MQ=255)
                                                                                                       tAGGAGCTGGTGCACCCGCGAAAGGGCCGGGCCTGCGTGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCGACCGTATTCGCGGGTAATCCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAccc                                           <  2:590134/151‑1 (MQ=14)
                                                                                                       tAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCAccc                                           <  1:3276144/151‑1 (MQ=255)
                                                                                                         ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                         >  1:1933904/1‑151 (MQ=255)
                                                                                                         ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                         >  1:1933816/1‑151 (MQ=255)
                                                                                                         ggAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                         >  1:1734243/1‑151 (MQ=255)
                                                                                                         ggAGCGGGTGCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                         <  1:4484701/151‑1 (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:1159365/1‑151 (MQ=255)
                                                                                                               ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                   <  1:1309753/151‑1 (MQ=255)
                                                                                                               ggTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAg                                   >  1:3132835/1‑151 (MQ=255)
                                                                                                                gtggagccGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGACGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAgg                                  <  2:2263540/145‑1 (MQ=14)
                                                                                                                 ttCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGc                                 >  1:1860971/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)
                                                                                                                          gAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCcgc                                         <  2:515672/134‑1 (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)
                                                                                                                                gccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGaaaa                  >  1:5293/1‑151 (MQ=255)
                                                                                                                                 ccgggccTGCTGGGCAATCTCGACTCGTCGTTCCACGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAAt                 <  2:3986931/151‑1 (MQ=18)
                                                                                                                                 ccgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGACGAAAAt                 >  2:3218224/1‑151 (MQ=255)
                                                                                                                                  cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGCAAATc                >  2:3334103/1‑151 (MQ=255)
                                                                                                                                  cgggccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATc                <  1:4349891/151‑1 (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:2827459/151‑1 (MQ=255)
                                                                                                                                     gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc             >  1:2711800/1‑151 (MQ=255)
                                                                                                                                     gccTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGCTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGc             >  2:4013095/1‑151 (MQ=255)
                                                                                                                                      ccTGCTGGCCAATTTCCCCCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc            <  1:319941/151‑1 (MQ=255)
                                                                                                                                      ccTGCGGGCCAATTTCCCCCCGCGGCTCACCGTTGCGGCGGCGACCGTGTTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTGGGGCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcc            <  2:1396955/151‑1 (MQ=14)
                                                                                                                                       ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc           >  1:4086133/1‑151 (MQ=255)
                                                                                                                                       ctgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGccc           >  1:892767/1‑151 (MQ=255)
                                                                                                                                        tgctgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGcccc          >  1:4865977/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:4757691/1‑151 (MQ=255)
                                                                                                                                           tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg       <  1:821614/151‑1 (MQ=255)
                                                                                                                                           tgGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTg       >  2:3407229/1‑151 (MQ=255)
                                                                                                                                            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  >  1:3090902/1‑151 (MQ=255)
                                                                                                                                                aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCtttt  <  2:724072/151‑1 (MQ=255)
                                                                                                                                                aaTTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCCCTTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGACAATCGGCCCCCTGCtttt  >  1:1639820/1‑151 (MQ=18)
                                                                                                                                                   |                                                                                                                                                   
GGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTT  >  CP009974/3837634‑3837928

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
CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCC  >  CP009974/3837675‑3837930
                                                                                                          |                                                                                                                                                     
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.2758630/151‑1 (MQ=26)
                                          GAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCG                                                                 >  SRR10918231.2742172/1‑151 (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.2069862/1‑151 (MQ=48)
                                                     AGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCG                                                      >  SRR10918231.2070065/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.1742211/1‑151 (MQ=60)
                                                                GGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                           >  SRR10918231.1942609/1‑151 (MQ=60)
                                                                GGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGC                                           >  SRR10918231.1942521/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.3147140/1‑151 (MQ=60)
                                                                      GGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAG                                     >  SRR10918231.1164756/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.2724277/1‑151 (MQ=60)
                                                                                            GCCTGCTGGCCAATTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGC               <  SRR10918231.2840408/151‑1 (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.4163027/1‑151 (MQ=60)
                                                                                                         TTTCCACCCGCCGTTCACCGTTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCC  <  SRR10918231.2377832/151‑1 (MQ=60)
                                                                                                          |                                                                                                                                                     
CGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACGCATCGAGGGCGGCACAGGCCGTTGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGTCGTTCACCATTGCCGCCGCTACCGTATTCGCGGGTAAACCCGCTCCTGCACACGCATCGGTGGTTGGTCACAGGCCGTCGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGTACAGGCGAAAATCGGCCCCCTGCTTTTCC  >  CP009974/3837675‑3837930

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: