Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 279,800 G→T 100% E297* (GAA→TAA)  yagF → CP4‑6 prophage; dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273279,8000GT96.9% 387.6 / 4.3 127E297* (GAA→TAA) yagFCP4‑6 prophage; dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (52/71);  total (54/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 127 total aligned reads displayed.

GCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG  >  CP009273/279655‑279943
                                                                                                                                                 |                                                                                                                                               
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:1563292/148‑1 (MQ=255)
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:1563299/148‑1 (MQ=255)
gcTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaa                                                                                                                                               <  2:3130368/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2524859/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2524860/148‑1 (MQ=255)
 cTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaa                                                                                                                                              <  1:2526469/148‑1 (MQ=255)
  tGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTaaaa                                                                                                                                             <  1:1232581/148‑1 (MQ=255)
      tCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAcgcg                                                                                                                                         <  2:1915160/148‑1 (MQ=255)
       cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                        <  1:3325367/148‑1 (MQ=255)
       cTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGa                                                                                                                                        <  1:3208521/148‑1 (MQ=255)
         ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATg                                                                                                                                      <  2:1643514/148‑1 (MQ=255)
         ggCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATg                                                                                                                                      >  2:2706577/1‑148 (MQ=255)
            gcTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTg                                                                                                                                   <  1:2161441/148‑1 (MQ=255)
             cTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGa                                                                                                                                  <  1:2841512/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  1:1926408/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  2:2993661/148‑1 (MQ=255)
               gCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATc                                                                                                                                <  1:1431759/148‑1 (MQ=255)
                 cGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCa                                                                                                                              <  1:2789103/148‑1 (MQ=255)
                      tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACgcgg                                                                                                                         <  1:859295/148‑1 (MQ=255)
                      tCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACgcgg                                                                                                                         <  1:981877/148‑1 (MQ=255)
                        cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCATCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                       <  2:1317823/148‑1 (MQ=255)
                        cgcgCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACgcggcg                                                                                                                       <  1:138919/148‑1 (MQ=255)
                         gcgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGt                                                                                                                      <  2:3130990/148‑1 (MQ=255)
                          cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                     >  1:2741139/1‑148 (MQ=255)
                          cgcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGtt                                                                                                                     <  1:427333/148‑1 (MQ=255)
                           gcTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTc                                                                                                                    <  2:1497848/148‑1 (MQ=255)
                             tGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTcgg                                                                                                                  <  2:252565/148‑1 (MQ=255)
                                      ccGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                         >  2:3363296/1‑148 (MQ=255)
                                      ccGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTccac                                                                                                         >  2:3363294/1‑148 (MQ=255)
                                         ggCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCaa                                                                                                      <  1:1816923/148‑1 (MQ=255)
                                              gCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAc                                                                                                 >  1:1830574/1‑148 (MQ=255)
                                              gCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAc                                                                                                 >  1:2804014/1‑148 (MQ=255)
                                                ggTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActg                                                                                               <  2:1073049/148‑1 (MQ=255)
                                                 gtgtGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTActgc                                                                                              >  1:1972745/1‑148 (MQ=255)
                                                    tgGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACtgctgc                                                                                           >  1:3205092/1‑148 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  2:359931/148‑1 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  2:2128577/148‑1 (MQ=255)
                                                      gCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGcac                                                                                         <  1:2041013/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:2278625/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:2408079/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:645279/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:1071772/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  2:2408130/148‑1 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     >  2:3360287/1‑148 (MQ=255)
                                                          gagaTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTc                                                                                     <  1:1071550/148‑1 (MQ=255)
                                                             aTCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCgg                                                                                  >  1:1556292/1‑148 (MQ=255)
                                                                     ccAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGcc                                                                          <  2:158392/148‑1 (MQ=255)
                                                                       aGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCa                                                                        <  2:2052313/148‑1 (MQ=255)
                                                                                gcgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCggg                                                               <  1:1757082/148‑1 (MQ=255)
                                                                                 cgcgGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGc                                                              <  2:963586/148‑1 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        >  1:609934/1‑148 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        <  2:2779261/148‑1 (MQ=255)
                                                                                       cAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACg                                                        >  1:2372820/1‑148 (MQ=255)
                                                                                         gCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGAt                                                      >  1:1590004/1‑148 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     <  2:600345/148‑1 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     <  1:1023389/148‑1 (MQ=255)
                                                                                          cGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAAGTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATc                                                     >  2:774789/1‑148 (MQ=255)
                                                                                                 gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                              <  1:3174245/148‑1 (MQ=255)
                                                                                                 gATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACg                                              <  1:3174270/148‑1 (MQ=255)
                                                                                                  aTAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGt                                             <  2:1814004/148‑1 (MQ=255)
                                                                                                    aGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTg                                           >  1:1092308/1‑148 (MQ=255)
                                                                                                      ccGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAg                                         >  2:801840/1‑148 (MQ=255)
                                                                                                            cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                   >  2:561949/1‑148 (MQ=255)
                                                                                                            cATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTgg                                   >  1:762285/1‑148 (MQ=255)
                                                                                                             aTCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGa                                  >  2:2898031/1‑148 (MQ=255)
                                                                                                               caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGTTGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                <  2:1409127/148‑1 (MQ=255)
                                                                                                               caccacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcg                                <  2:3124939/148‑1 (MQ=255)
                                                                                                                accacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGAcgc                               >  1:1391642/1‑148 (MQ=255)
                                                                                                                  cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                             <  2:2717558/148‑1 (MQ=255)
                                                                                                                  cacGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACgcgc                             <  2:2717632/148‑1 (MQ=255)
                                                                                                                    cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                           <  2:1830574/148‑1 (MQ=255)
                                                                                                                    cgcgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCAt                           <  1:1607365/148‑1 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  1:3110161/1‑148 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  1:2334077/1‑148 (MQ=255)
                                                                                                                      cgGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCa                         >  2:3022812/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:665398/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:1067703/1‑148 (MQ=255)
                                                                                                                         gATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAAcc                      >  2:665799/1‑148 (MQ=255)
                                                                                                                            atCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTa                   >  1:1698607/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  >  1:713015/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  >  1:713200/1‑148 (MQ=255)
                                                                                                                             tCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaa                  <  2:1216663/148‑1 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  1:3464529/1‑148 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  1:3464482/1‑148 (MQ=255)
                                                                                                                              ccTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTaaa                 >  2:1048007/1‑148 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  2:897448/148‑1 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                >  2:2425154/1‑148 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  1:872997/148‑1 (MQ=255)
                                                                                                                               ctctCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATGTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAg                <  1:873242/148‑1 (MQ=255)
                                                                                                                                tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt               >  1:3282588/1‑148 (MQ=255)
                                                                                                                                tctcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGt               >  1:3282592/1‑148 (MQ=255)
                                                                                                                                  tcCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGc             <  2:3171471/148‑1 (MQ=255)
                                                                                                                                   ccGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGcc            >  2:1314707/1‑148 (MQ=255)
                                                                                                                                       tAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGt        >  1:3118594/1‑148 (MQ=255)
                                                                                                                                        aaaGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTc       >  2:610456/1‑148 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     <  1:1471709/148‑1 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     <  2:1972745/148‑1 (MQ=255)
                                                                                                                                          aGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTg     >  1:411351/1‑148 (MQ=255)
                                                                                                                                             cATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTg  >  1:3075/1‑148 (MQ=255)
                                                                                                                                             cATCTAAAACGCGATGATGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTg  >  1:2810479/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG  >  CP009273/279655‑279943

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  CP009273/279656‑279949
                                                                                                                                                |                                                                                                                                                     
CTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAC                                                                                                                                                  <  SRR19973065.2542332/150‑1 (MQ=60)
CTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAC                                                                                                                                                  <  SRR19973065.2540708/150‑1 (MQ=60)
CTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAAC                                                                                                                                                  <  SRR19973065.2540707/150‑1 (MQ=60)
 TGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACG                                                                                                                                                 <  SRR19973065.1240654/150‑1 (MQ=60)
      CTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATG                                                                                                                                            <  SRR19973065.3345658/150‑1 (MQ=60)
      CTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATG                                                                                                                                            <  SRR19973065.3227828/150‑1 (MQ=60)
          CGCTGCCGCACTACGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGA                                                                                                                                        <  SRR19973065.2160433/150‑1 (MQ=60)
           GCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGAT                                                                                                                                       <  SRR19973065.2175336/150‑1 (MQ=60)
            CTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATC                                                                                                                                      <  SRR19973065.2859005/150‑1 (MQ=60)
              GCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCA                                                                                                                                    <  SRR19973065.1441068/150‑1 (MQ=60)
              GCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCA                                                                                                                                    <  SRR19973065.1939032/150‑1 (MQ=60)
                CGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACG                                                                                                                                  <  SRR19973065.2806391/150‑1 (MQ=60)
                   ACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGG                                                                                                                               >  SRR19973065.3143969/1‑150 (MQ=60)
                     TCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCG                                                                                                                             <  SRR19973065.988393/150‑1 (MQ=60)
                     TCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCG                                                                                                                             <  SRR19973065.864834/150‑1 (MQ=60)
                       CGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTT                                                                                                                           >  SRR19973065.2758186/1‑150 (MQ=60)
                       CGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTT                                                                                                                           <  SRR19973065.140107/150‑1 (MQ=60)
                         CGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCG                                                                                                                         <  SRR19973065.430304/150‑1 (MQ=60)
                                       GGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAAT                                                                                                           <  SRR19973065.779998/150‑1 (MQ=60)
                                        GGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATT                                                                                                          <  SRR19973065.1828891/150‑1 (MQ=60)
                                           AGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAC                                                                                                       >  SRR19973065.2821372/1‑150 (MQ=60)
                                           AGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTAC                                                                                                       >  SRR19973065.1842634/1‑150 (MQ=60)
                                              CGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGC                                                                                                    >  SRR19973065.1985595/1‑150 (MQ=60)
                                                 TGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGC                                                                                                 >  SRR19973065.3224385/1‑150 (MQ=60)
                                                     GCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACAT                                                                                             <  SRR19973065.2054079/150‑1 (MQ=60)
                                                       TGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTC                                                                                           <  SRR19973065.2840073/150‑1 (MQ=60)
                                                         GAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCG                                                                                         <  SRR19973065.1078865/150‑1 (MQ=60)
                                                         GAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCG                                                                                         <  SRR19973065.2292998/150‑1 (MQ=60)
                                                         GAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCG                                                                                         <  SRR19973065.1078640/150‑1 (MQ=60)
                                                          AGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGG                                                                                        >  SRR19973065.1566616/1‑150 (MQ=60)
                                                                      AGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACG                                                                            >  SRR19973065.3484050/1‑150 (MQ=60)
                                                                              CGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGC                                                                    >  SRR19973065.903185/1‑150 (MQ=60)
                                                                               GCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCT                                                                   <  SRR19973065.1768630/150‑1 (MQ=60)
                                                                                    GTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACG                                                              >  SRR19973065.613818/1‑150 (MQ=60)
                                                                                    GTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACG                                                              >  SRR19973065.2387632/1‑150 (MQ=60)
                                                                                      CAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGAT                                                            >  SRR19973065.1600617/1‑150 (MQ=60)
                                                                                         CGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCC                                                         <  SRR19973065.1030167/150‑1 (MQ=60)
                                                                                                GATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTT                                                  <  SRR19973065.3193441/150‑1 (MQ=60)
                                                                                                GATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTT                                                  <  SRR19973065.3193466/150‑1 (MQ=60)
                                                                                                 ATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTG                                                 >  SRR19973065.1099501/1‑150 (MQ=60)
                                                                                                  TAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGA                                                <  SRR19973065.2658342/150‑1 (MQ=60)
                                                                                                         GGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGG                                         >  SRR19973065.767361/1‑150 (MQ=60)
                                                                                                            ATCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACG                                      >  SRR19973065.2444641/1‑150 (MQ=60)
                                                                                                             TCACCACGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGC                                     >  SRR19973065.1400777/1‑150 (MQ=60)
                                                                                                                   CGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCA                               <  SRR19973065.1618124/150‑1 (MQ=60)
                                                                                                                   CGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCA                               >  SRR19973065.2348660/1‑150 (MQ=60)
                                                                                                                   CGCGGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCA                               >  SRR19973065.3129122/1‑150 (MQ=60)
                                                                                                                      GGGATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACC                            >  SRR19973065.1832026/1‑150 (MQ=60)
                                                                                                                         ATATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTA                         >  SRR19973065.1709934/1‑150 (MQ=60)
                                                                                                                          TATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAA                        >  SRR19973065.1629640/1‑150 (MQ=60)
                                                                                                                          TATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAA                        >  SRR19973065.717973/1‑150 (MQ=60)
                                                                                                                          TATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAA                        >  SRR19973065.717787/1‑150 (MQ=60)
                                                                                                                           ATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAA                       >  SRR19973065.3485581/1‑150 (MQ=60)
                                                                                                                           ATCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAA                       >  SRR19973065.3485629/1‑150 (MQ=60)
                                                                                                                            TCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAG                      >  SRR19973065.1117275/1‑150 (MQ=60)
                                                                                                                            TCCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAG                      >  SRR19973065.508394/1‑150 (MQ=60)
                                                                                                                             CCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGT                     >  SRR19973065.3302543/1‑150 (MQ=60)
                                                                                                                             CCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGT                     >  SRR19973065.3302548/1‑150 (MQ=60)
                                                                                                                             CCTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGT                     >  SRR19973065.3302544/1‑150 (MQ=60)
                                                                                                                              CTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTG                    <  SRR19973065.878583/150‑1 (MQ=60)
                                                                                                                              CTCTCCGATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATGTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTG                    <  SRR19973065.878828/150‑1 (MQ=60)
                                                                                                                                    GATAAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGT              >  SRR19973065.3137600/1‑150 (MQ=60)
                                                                                                                                       AAAGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTG           >  SRR19973065.414231/1‑150 (MQ=60)
                                                                                                                                         AGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT         <  SRR19973065.1481136/150‑1 (MQ=60)
                                                                                                                                         AGCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGT         <  SRR19973065.1481119/150‑1 (MQ=60)
                                                                                                                                          GCCATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG        >  SRR19973065.3112/1‑150 (MQ=60)
                                                                                                                                          GCCATCTAAAACGCGATGATGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTG        >  SRR19973065.2827860/1‑150 (MQ=60)
                                                                                                                                            CATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAG      >  SRR19973065.1583924/1‑150 (MQ=60)
                                                                                                                                             ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973065.2602287/1‑150 (MQ=60)
                                                                                                                                             ATCTAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGC     >  SRR19973065.831892/1‑150 (MQ=60)
                                                                                                                                                TAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  <  SRR19973065.2832281/150‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                                     
CTGGGTCTGGCGCTGCCGCACTCCGCGCTGGCGCCGTCCGGGCAGGCGGTGTGGCTGGAGATCGCCCGCCAGTCGGCGCGCGCGGTCAGCGAGCTGGATAGCCGCGGCATCACCACGCGGGATATCCTCTCCGATAAAGCCATCGAAAACGCGATGGTGATCCACGCGGCGTTCGGCGGCTCCACCAATTTACTGCTGCACATTCCGGCCATCGCCCACGCGGCGGGCTGCACGATCCCGGACGTTGAGCACTGGACGCGCATCAACCGTAAAGTGCCGCGTCTGGTGAGCGTG  >  CP009273/279656‑279949

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