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 2,322,609 C→A 100% W223L (TGG→TTG)  yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,322,6090CA99.3% 458.6 / ‑4.1 140W223L (TGG→TTG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (73/66);  total (73/67)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 141 total aligned reads displayed.

GTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCAC  >  CP009273/2322465‑2322751
                                                                                                                                                |                                                                                                                                              
      ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                       <  1:1967863/148‑1 (MQ=255)
      ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                       <  1:712094/148‑1 (MQ=255)
      ttcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                       <  1:712633/148‑1 (MQ=255)
       tcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtgt                                                                                                                                      <  1:1591720/148‑1 (MQ=255)
         ttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAc                                                                                                                                    <  1:2824134/148‑1 (MQ=255)
          tcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAcc                                                                                                                                   <  2:2573517/148‑1 (MQ=255)
            aGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAt                                                                                                                                 <  1:888199/148‑1 (MQ=255)
            aGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAt                                                                                                                                 <  1:2646045/148‑1 (MQ=255)
            aGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAt                                                                                                                                 >  1:642087/1‑148 (MQ=255)
               ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  1:2549674/148‑1 (MQ=255)
               ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                              <  1:2376654/148‑1 (MQ=255)
                gATAGCGGGCAAACAGCTCCACCTCTACACAGTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGt                                                                                                                             <  2:2403127/148‑1 (MQ=255)
                    gCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAc                                                                                                                         >  2:166236/1‑148 (MQ=255)
                     cGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAcc                                                                                                                        >  1:2150175/1‑148 (MQ=255)
                     cGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAcc                                                                                                                        >  1:2150257/1‑148 (MQ=255)
                       ggCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccag                                                                                                                      >  1:3152112/1‑148 (MQ=255)
                        gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagc                                                                                                                     <  1:2523309/148‑1 (MQ=255)
                        gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagc                                                                                                                     <  2:2458935/148‑1 (MQ=255)
                         cAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcc                                                                                                                    <  2:2009412/148‑1 (MQ=255)
                             cAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCa                                                                                                                >  1:587874/1‑148 (MQ=255)
                               gCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaa                                                                                                              >  1:1622940/1‑148 (MQ=255)
                                cTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaaa                                                                                                             <  1:964287/148‑1 (MQ=255)
                                cTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaaa                                                                                                             <  1:2718144/148‑1 (MQ=255)
                                   cACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacca                                                                                                          <  2:2984960/148‑1 (MQ=255)
                                     ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                        <  2:3379386/148‑1 (MQ=255)
                                     ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                        >  2:2308814/1‑148 (MQ=255)
                                      ctctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                       <  2:2305080/148‑1 (MQ=255)
                                       tctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCt                                                                                                      >  2:1311147/1‑148 (MQ=255)
                                        ctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTg                                                                                                     <  2:413948/148‑1 (MQ=255)
                                          acacATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAc                                                                                                   >  1:270792/1‑148 (MQ=255)
                                          acacATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAc                                                                                                   <  1:2297206/148‑1 (MQ=255)
                                           cacaTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACg                                                                                                  >  1:2242888/1‑148 (MQ=255)
                                             caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGa                                                                                                <  2:3450153/148‑1 (MQ=255)
                                             caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGa                                                                                                <  2:443428/148‑1 (MQ=255)
                                                    cACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTg                                                                                         <  1:1069879/148‑1 (MQ=255)
                                                     aCTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTTCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGt                                                                                        >  1:1127817/1‑148 (MQ=255)
                                                      cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGtt                                                                                       >  2:2655669/1‑148 (MQ=255)
                                                      cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGtt                                                                                       <  2:470268/148‑1 (MQ=255)
                                                      cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGtt                                                                                       >  2:890548/1‑148 (MQ=255)
                                                        ggCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTa                                                                                     >  2:348318/1‑148 (MQ=255)
                                                            aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCACAACCACCCTGACGGATTTCCTGTTTatcat                                                                                 >  2:2709687/1‑148 (MQ=255)
                                                            aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcat                                                                                 >  1:2285649/1‑148 (MQ=255)
                                                            aCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcat                                                                                 >  2:2911031/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  1:371479/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  1:371448/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:825492/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  1:841480/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  1:371436/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:1624771/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:401693/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:825530/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:1403578/1‑148 (MQ=255)
                                                             ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                >  2:1403284/1‑148 (MQ=255)
                                                                 cGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCa                                                                            <  1:1338689/148‑1 (MQ=255)
                                                                    cGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAc                                                                         <  1:2830168/148‑1 (MQ=255)
                                                                    cGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAc                                                                         >  2:390782/1‑148 (MQ=255)
                                                                     ggCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAAcc                                                                        <  2:3046294/148‑1 (MQ=255)
                                                                      gCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCa                                                                       <  2:452149/148‑1 (MQ=255)
                                                                           ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTcc                                                                  <  2:535198/148‑1 (MQ=255)
                                                                              aaaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAgg                                                               <  2:2209821/148‑1 (MQ=255)
                                                                                aaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                             >  2:1426032/1‑148 (MQ=255)
                                                                                aaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                             >  2:1426371/1‑148 (MQ=255)
                                                                                 aCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGttt                                                            >  1:1245131/1‑148 (MQ=255)
                                                                                  cccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtttt                                                           >  1:1539987/1‑148 (MQ=255)
                                                                                  cccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtttt                                                           <  1:1660315/148‑1 (MQ=255)
                                                                                    cGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCa                                                         >  1:1668112/1‑148 (MQ=255)
                                                                                    cGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCa                                                         >  1:1258334/1‑148 (MQ=255)
                                                                                    cGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCa                                                         >  1:1258408/1‑148 (MQ=255)
                                                                                      aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                       <  2:1609047/148‑1 (MQ=255)
                                                                                      aTTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAcc                                                       <  2:1608042/148‑1 (MQ=255)
                                                                                       ttAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCAccc                                                      >  2:1483033/1‑148 (MQ=255)
                                                                                           gCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTAc                                                  <  1:1480015/148‑1 (MQ=255)
                                                                                              ccGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCg                                               >  2:231427/1‑148 (MQ=255)
                                                                                              ccGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCg                                               >  2:1722526/1‑148 (MQ=255)
                                                                                               cGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCACCAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGa                                              <  1:2911031/148‑1 (MQ=255)
                                                                                                 tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACt                                            >  2:2512882/1‑148 (MQ=255)
                                                                                                 tCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACt                                            >  2:1963839/1‑148 (MQ=255)
                                                                                                  cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGTTTAACCAGCCAGCAAAACCACCCTGACTGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTa                                           <  1:1078954/148‑1 (MQ=255)
                                                                                                  cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTa                                           <  1:1078665/148‑1 (MQ=255)
                                                                                                    cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCANCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAcc                                         <  1:1784070/148‑1 (MQ=255)
                                                                                                    cATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTAcc                                         >  1:295876/1‑148 (MQ=255)
                                                                                                     aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCg                                        >  2:1351182/1‑148 (MQ=255)
                                                                                                      tAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGa                                       >  2:116952/1‑148 (MQ=255)
                                                                                                       aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                      >  1:1084574/1‑148 (MQ=255)
                                                                                                       aaCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                      >  1:1851738/1‑148 (MQ=255)
                                                                                                        aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAGTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGc                                     >  2:2148040/1‑148 (MQ=255)
                                                                                                           ggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATg                                  >  1:341231/1‑148 (MQ=255)
                                                                                                            gTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATgg                                 <  2:1539987/148‑1 (MQ=255)
                                                                                                            gTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATgg                                 <  2:2150257/148‑1 (MQ=255)
                                                                                                             ttAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTAACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGa                                <  1:2652155/148‑1 (MQ=255)
                                                                                                               aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATg                              <  2:898944/148‑1 (MQ=255)
                                                                                                                 tACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                            >  1:982747/1‑148 (MQ=255)
                                                                                                                   ccGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGac                          <  2:1063869/148‑1 (MQ=255)
                                                                                                                       tttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                      >  2:3161960/1‑148 (MQ=255)
                                                                                                                       tttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                      >  2:632392/1‑148 (MQ=255)
                                                                                                                           cTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAttt                  <  1:1645889/148‑1 (MQ=255)
                                                                                                                                gCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGcc             <  1:1311166/148‑1 (MQ=255)
                                                                                                                                    cTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGt         <  2:3148854/148‑1 (MQ=255)
                                                                                                                                         aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGc    >  1:776820/1‑148 (MQ=255)
                                                                                                                                          gTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGca   >  1:3421958/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCAC  >  CP009273/2322465‑2322751

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
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  CP009273/2322462‑2322754
                                                                                                                                                   |                                                                                                                                                 
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAAC                                                                                                                                                 <  SRR19973065.2968465/150‑1 (MQ=60)
         TTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTA                                                                                                                                        <  SRR19973065.1980700/150‑1 (MQ=60)
         TTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTA                                                                                                                                        <  SRR19973065.717401/150‑1 (MQ=60)
         TTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTA                                                                                                                                        <  SRR19973065.716857/150‑1 (MQ=60)
          TCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAC                                                                                                                                       <  SRR19973065.1602353/150‑1 (MQ=60)
            TTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCA                                                                                                                                     <  SRR19973065.2841571/150‑1 (MQ=60)
             TCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCAT                                                                                                                                    >  SRR19973065.646300/1‑150 (MQ=60)
               AGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCG                                                                                                                                  <  SRR19973065.893837/150‑1 (MQ=60)
               AGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCG                                                                                                                                  <  SRR19973065.2662360/150‑1 (MQ=60)
                  GGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTT                                                                                                                               <  SRR19973065.2391503/150‑1 (MQ=60)
                  GGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTT                                                                                                                               <  SRR19973065.2565634/150‑1 (MQ=60)
                      AGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACC                                                                                                                           >  SRR19973065.2164086/1‑150 (MQ=60)
                      AGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACC                                                                                                                           >  SRR19973065.2164004/1‑150 (MQ=60)
                      AGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACC                                                                                                                           <  SRR19973065.680612/150‑1 (MQ=60)
                        CGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAG                                                                                                                         >  SRR19973065.3171237/1‑150 (MQ=60)
                           GCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCA                                                                                                                      <  SRR19973065.2539151/150‑1 (MQ=60)
                              AACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCA                                                                                                                   >  SRR19973065.591526/1‑150 (MQ=60)
                                CAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAA                                                                                                                 >  SRR19973065.1633820/1‑150 (MQ=60)
                                   CTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACC                                                                                                              <  SRR19973065.2735051/150‑1 (MQ=60)
                                   CTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACC                                                                                                              <  SRR19973065.970684/150‑1 (MQ=60)
                                           CTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAC                                                                                                      >  SRR19973065.272623/1‑150 (MQ=60)
                                           CTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGAC                                                                                                      <  SRR19973065.564251/150‑1 (MQ=60)
                                            TACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACG                                                                                                     >  SRR19973065.2257135/1‑150 (MQ=60)
                                             ACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGG                                                                                                    <  SRR19973065.2311668/150‑1 (MQ=60)
                                              CACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGA                                                                                                   <  SRR19973065.1603526/150‑1 (MQ=60)
                                                CATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATT                                                                                                 <  SRR19973065.33515/150‑1 (MQ=60)
                                                 ATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTTCCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTT                                                                                                <  SRR19973065.2947405/150‑1 (MQ=60)
                                                    GCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTGAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCT                                                                                             <  SRR19973065.2552684/150‑1 (MQ=60)
                                                      CCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTTCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGT                                                                                           >  SRR19973065.1135152/1‑150 (MQ=60)
                                                       CACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTT                                                                                          <  SRR19973065.1076960/150‑1 (MQ=60)
                                                         CTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTA                                                                                        <  SRR19973065.1840116/150‑1 (MQ=60)
                                                             CAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCAT                                                                                    >  SRR19973065.1998580/1‑150 (MQ=60)
                                                             CAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCAT                                                                                    >  SRR19973065.2300057/1‑150 (MQ=60)
                                                              AACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATC                                                                                   >  SRR19973065.374047/1‑150 (MQ=60)
                                                              AACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATC                                                                                   >  SRR19973065.846956/1‑150 (MQ=60)
                                                              AACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATC                                                                                   >  SRR19973065.374036/1‑150 (MQ=60)
                                                              AACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATC                                                                                   >  SRR19973065.374048/1‑150 (MQ=60)
                                                              AACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATC                                                                                   >  SRR19973065.374079/1‑150 (MQ=60)
                                                                    CGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAA                                                                             <  SRR19973065.1347552/150‑1 (MQ=60)
                                                                       CGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCA                                                                          <  SRR19973065.2847624/150‑1 (MQ=60)
                                                                                  AAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTT                                                               >  SRR19973065.1253325/1‑150 (MQ=60)
                                                                                   AACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTT                                                              >  SRR19973065.1550160/1‑150 (MQ=60)
                                                                                     CCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCA                                                            >  SRR19973065.1679275/1‑150 (MQ=60)
                                                                                     CCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCA                                                            <  SRR19973065.1671437/150‑1 (MQ=60)
                                                                                     CCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCA                                                            >  SRR19973065.1266698/1‑150 (MQ=60)
                                                                                     CCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCA                                                            >  SRR19973065.1266624/1‑150 (MQ=60)
                                                                                           TAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTGTTCACCCGTA                                                      >  SRR19973065.2165791/1‑150 (MQ=60)
                                                                                              GCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCC                                                   <  SRR19973065.1489472/150‑1 (MQ=60)
                                                                                                  CGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCACCAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACT                                               <  SRR19973065.2928750/150‑1 (MQ=60)
                                                                                                     CGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGTTTAACCAGCCAGCAAAACCACCCTGACTGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACC                                            <  SRR19973065.1086080/150‑1 (MQ=60)
                                                                                                     CGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACC                                            <  SRR19973065.1085789/150‑1 (MQ=60)
                                                                                                     CGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACC                                            >  SRR19973065.297777/1‑150 (MQ=60)
                                                                                                       CATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCANCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGA                                          <  SRR19973065.1795737/150‑1 (MQ=60)
                                                                                                       CATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGA                                          >  SRR19973065.2826905/1‑150 (MQ=60)
                                                                                                        ATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAG                                         >  SRR19973065.1091728/1‑150 (MQ=60)
                                                                                                        ATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAG                                         >  SRR19973065.1091729/1‑150 (MQ=60)
                                                                                                        ATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAG                                         >  SRR19973065.1863937/1‑150 (MQ=60)
                                                                                                            CGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATG                                     >  SRR19973065.343615/1‑150 (MQ=60)
                                                                                                                TTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATG                                 <  SRR19973065.1423498/150‑1 (MQ=60)
                                                                                                                TTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATG                                 <  SRR19973065.2323319/150‑1 (MQ=60)
                                                                                                                TTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTAACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATG                                 <  SRR19973065.2668524/150‑1 (MQ=60)
                                                                                                                  AGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCG                               >  SRR19973065.989273/1‑150 (MQ=60)
                                                                                                                     ACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACA                            <  SRR19973065.916787/150‑1 (MQ=60)
                                                                                                                          TTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGCTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAT                       >  SRR19973065.782793/1‑150 (MQ=60)
                                                                                                                              CTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTC                   <  SRR19973065.1656983/150‑1 (MQ=60)
                                                                                                                              CTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTC                   <  SRR19973065.1656937/150‑1 (MQ=60)
                                                                                                                                   GCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTC              <  SRR19973065.1319825/150‑1 (MQ=60)
                                                                                                                                   GCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTC              <  SRR19973065.1319806/150‑1 (MQ=60)
                                                                                                                                          CCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGC       >  SRR19973065.781984/1‑150 (MQ=60)
                                                                                                                                           CAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCA      >  SRR19973065.3442857/1‑150 (MQ=60)
                                                                                                                                            AGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCAC     >  SRR19973065.2722398/1‑150 (MQ=60)
                                                                                                                                              TAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGAACCA   <  SRR19973065.123701/150‑1 (MQ=60)
                                                                                                                                               AACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  SRR19973065.50775/1‑150 (MQ=60)
                                                                                                                                               AACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  SRR19973065.3418803/1‑150 (MQ=60)
                                                                                                                                               AACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  SRR19973065.3418348/1‑150 (MQ=60)
                                                                                                                                                   |                                                                                                                                                 
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCAC  >  CP009273/2322462‑2322754

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