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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,409 C→T 100% intergenic (+448/+166) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,4090CT90.7% 229.2 / 16.0 109intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (5/5);  new base T (43/53);  total (50/59)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 111 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552
                                                                                                                                                |                                                                                                                                               
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                              <  1:749842/148‑1 (MQ=37)
 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                             <  1:1207133/148‑1 (MQ=37)
   gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTg                                                                                                                                           <  1:2139031/148‑1 (MQ=37)
   gggTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTg                                                                                                                                           <  1:2138977/148‑1 (MQ=37)
    gggataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGc                                                                                                                                          <  2:3452917/145‑1 (MQ=255)
     gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGcc                                                                                                                                         <  2:1526375/148‑1 (MQ=37)
      tataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGccc                                                                                                                                        <  1:443323/148‑1 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                     <  1:845757/148‑1 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                                     >  2:2044908/1‑148 (MQ=37)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                                                                                                                                     >  2:2044993/1‑148 (MQ=37)
               gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAg                                                                                                                               <  2:2074020/148‑1 (MQ=37)
               gTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAg                                                                                                                               <  2:2073933/148‑1 (MQ=37)
                 ggACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCa                                                                                                                             >  2:2337754/1‑148 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                          <  2:468290/148‑1 (MQ=37)
                    cAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCaaaa                                                                                                                          <  2:2244780/148‑1 (MQ=37)
                     aGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCATTTGGGTATAGTGTAGTGGAAAGTCATTAATTTTTCTTCCCCTCcaaaaacaaaaa                                                                                                                         >  2:2102984/1‑148 (MQ=11)
                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAAcc                                                                                                                       >  1:2068598/1‑148 (MQ=255)
                        cattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                      <  2:2805671/148‑1 (MQ=37)
                        cagtcatCTTTCTGCCCCTCCAAACGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAccc                                                                                                                      <  2:1871892/145‑1 (MQ=255)
                            catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATGTTTGTGGCCCTCCAAAAGCAAAAACcagccg                                                                                                                  >  2:668417/1‑148 (MQ=255)
                              tCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                <  2:285938/148‑1 (MQ=37)
                                    tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                          <  2:825150/148‑1 (MQ=37)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                         <  1:596045/148‑1 (MQ=37)
                                     gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtt                                                                                                         <  1:2384855/148‑1 (MQ=37)
                                      ccccTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  1:3371255/148‑1 (MQ=12)
                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  1:3371262/148‑1 (MQ=37)
                                       cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                       >  1:278796/1‑148 (MQ=37)
                                         cTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGGCCCTGCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTa                                                                                                     >  1:1958685/1‑148 (MQ=11)
                                             aaaaGAACAAACCCGCCGCAGCGGGTTATTACGTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTa                                                                                                 <  2:1818358/148‑1 (MQ=11)
                                              taaGTAACACCCCGCCGAAGCGGGTTTTTACGTAACACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTCTAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  2:2096905/147‑1 (MQ=11)
                                              aaaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  1:416275/148‑1 (MQ=37)
                                              aaaCTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  1:416395/148‑1 (MQ=25)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                               >  2:2953017/1‑148 (MQ=255)
                                               aaGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaa                                                                                               <  2:2821523/148‑1 (MQ=37)
                                                    aaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                          >  2:12479/1‑148 (MQ=12)
                                                     aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGt                                                                                         <  2:1513895/148‑1 (MQ=37)
                                                     aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCTTGGCCAGTCCGTCATTGTTATGCCCCGCCCAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTcaaacgggg                                                                                         >  2:137441/1‑142 (MQ=11)
                                                       aCCCGCCGAAGCGGGTTTTTACGTACAACCGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                                                                       <  1:448094/148‑1 (MQ=255)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCCAGCAATTTGCTCCCCCACCAAAAGCAGAATCCCGACGTAGCGGGGTTTTCAGTAAACAAGGTgt                                                                                       >  2:1408458/1‑147 (MQ=11)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGa                                                                                       <  2:278796/148‑1 (MQ=17)
                                                             cttatcGGTTTTTTAAGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                 <  2:841431/143‑1 (MQ=11)
                                                                 gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGa                                                                             >  2:2109857/1‑148 (MQ=255)
                                                                    ggTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTATGCCCCGCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATaa                                                                          >  1:1315007/1‑148 (MQ=11)
                                                                     gTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTTAACCAGGTGAAACTGCCCGCTAAg                                                                         >  2:2885713/1‑148 (MQ=255)
                                                                      tttttACGTACAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGc                                                                        <  2:151460/148‑1 (MQ=255)
                                                                      tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                        >  1:240304/1‑148 (MQ=12)
                                                                         ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                     <  2:3273678/148‑1 (MQ=38)
                                                                           aCGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctt                                                                   >  2:1741323/1‑148 (MQ=255)
                                                                            cGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGcttt                                                                  <  2:711337/148‑1 (MQ=255)
                                                                              tacatcAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                <  2:3183250/143‑1 (MQ=12)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcgttct                                                                >  1:3179427/1‑148 (MQ=255)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttct                                                                <  1:641266/148‑1 (MQ=255)
                                                                               aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttctt                                                               <  2:2153618/148‑1 (MQ=255)
                                                                                aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                              <  1:355037/148‑1 (MQ=255)
                                                                                aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                              <  1:2147413/148‑1 (MQ=255)
                                                                                 caCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:219659/147‑1 (MQ=255)
                                                                                 aaCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCACAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:1744765/148‑1 (MQ=255)
                                                                                 aaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:220449/148‑1 (MQ=38)
                                                                                 aaCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                             <  2:342198/148‑1 (MQ=16)
                                                                                   cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAGCAGGTTAAACTGACCGATAAGCCGGTTTCTTTTgg                                                           >  2:151414/1‑148 (MQ=255)
                                                                                     ggTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTgggg                                                         >  2:708967/1‑147 (MQ=255)
                                                                                       tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtat                                                       <  2:926729/148‑1 (MQ=25)
                                                                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                      >  2:3318317/1‑148 (MQ=2)
                                                                                         aaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCTGGTTTTTACGTTAACCAGGTGAAACAGACCCATAAGCCGCGTTATTTTGGGtattg                                                     >  2:431780/1‑146 (MQ=11)
                                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                  >  2:3198406/1‑148 (MQ=37)
                                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgt                                                  <  1:1044662/148‑1 (MQ=37)
                                                                                                ccGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTAAATTTTGTTCCCCACCAAAAACAAAAAAACCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAACAGACCCATAATCCGCTTTGTTTTGGGTGTAGCGTCGTg                                              >  1:1235600/1‑148 (MQ=2)
                                                                                                 cGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                             <  1:703004/148‑1 (MQ=37)
                                                                                                   aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTATGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGagtcattgaa                                           >  2:2312499/1‑142 (MQ=11)
                                                                                                    tAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTATGCCCCACCAAAATCAAAAAAACGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCCCTAAGCCGCTTTCTTTTGGGTATAGCGTCGTggtac                                          >  2:2703706/1‑145 (MQ=11)
                                                                                                     aaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGAACGATAAGCCGCTTTCTTTTTGGTATAGCGTCGTGGCCAg                                         >  1:275841/1‑148 (MQ=11)
                                                                                                       gccgccTACTTTTGTGTACAGTGTGGAGGACAGACATTCAACCTTCAGCCCCCCCAACCGCACTAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                       <  2:761209/148‑1 (MQ=11)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                       <  1:2337754/148‑1 (MQ=37)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTATGCCCCTCCAAAATCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAAAAGCCGCTTTGTTTTGGGTATAGCGTCGTGGGCCGtca                                      >  2:2480091/1‑148 (MQ=11)
                                                                                                        ccGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtca                                      >  2:3189492/1‑148 (MQ=255)
                                                                                                         cGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAAAAGACCGATAACCCGCTTTCTTTTTGGTATATCGTCGGGGagaggtaa                                     >  1:3332236/1‑144 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAATGAACGATAAGCCGCGTTCTGTTGGGGATAGCGTATGGGACAGgcaat                                    >  2:3243189/1‑146 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTCACCCACCAAAAGCAAAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTGAAACTTACCGATAAGCCGCTTAGTTGTGGGTATAGCGTCGTGGCCAGtcctt                                    >  2:3242980/1‑148 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       <  1:1804563/65‑1 (MQ=35)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       >  2:1804563/1‑65 (MQ=35)
                                                                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTTGAAAGtcaatcaaaa                               >  1:2980481/1‑145 (MQ=11)
                                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                              <  1:2044993/148‑1 (MQ=37)
                                                                                                                ttttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCtt                              <  1:2044908/148‑1 (MQ=37)
                                                                                                                 tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                             <  1:2109857/148‑1 (MQ=37)
                                                                                                                 tttGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCttt                             <  1:2109847/148‑1 (MQ=255)
                                                                                                                     ggTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATCtttaagc                         >  2:1955826/1‑144 (MQ=255)
                                                                                                                       tataGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAGCGTCGTGGCACGTCAATCATAATTCTGgcc                       >  2:2422766/1‑148 (MQ=11)
                                                                                                                        ataGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAACAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAATCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCAATCAAATTTGTGcccc                      >  1:60830/1‑148 (MQ=11)
                                                                                                                         gcgtgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCt                     <  1:331288/146‑1 (MQ=37)
                                                                                                                          aGTGTCGTGGACAGTCATTCATTTGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTc                    <  1:2441435/148‑1 (MQ=255)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaa           >  1:883087/1‑148 (MQ=37)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATAATTCTGCCCCTCCAAAAGCaa           >  2:1150348/1‑148 (MQ=255)
                                                                                                                                    aCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaa          >  1:3432054/1‑148 (MQ=37)
                                                                                                                                    aCAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGCGTCGTGGAAAGTCATTCATCTTTATGCCCCTCCAAAAGCaaa          >  2:1974627/1‑148 (MQ=255)
                                                                                                                                      aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaac        >  2:2897636/1‑147 (MQ=37)
                                                                                                                                      aGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTCATCTTTCTGCCCCGCCAAAAGCaaaaa        >  1:3179184/1‑148 (MQ=255)
                                                                                                                                         cattcattTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCCTCGTGGACAGTAATAAAACTTTATGCCCCGCAAAAAGCAAAAAacc     >  1:1178427/1‑148 (MQ=11)
                                                                                                                                           ttcatTTTTCTGCCCCTCCACACGCAACAACCCGCCGAAGCGGGTTTTTACGTCTTCCCGGTGACACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc   <  2:1672982/148‑1 (MQ=11)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc  >  2:290724/1‑148 (MQ=37)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATAAAATTTCTGCCCCACCAAAATCAAAAACccgcc  >  1:833932/1‑148 (MQ=11)
                                                                                                                                                |                                                                                                                                               
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCG  >  CP009273/3263265‑3263558
                                                                                                                                                |                                                                                                                                                     
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCT                                                                                                                                                  <  SRR19973066.755407/150‑1 (MQ=32)
 ccccccccccccccccccacccccacccaacccaccccccccccaaaagcaaaaacccgccgaatcccgatatatccttggtagtgactgcatttcagacgtgtgcccttccgttcGGTATAGTGTCGTGGACAGTCATTCATCTTTCTG                                                                                                                                                 <  SRR19973066.2082038/34‑1 (MQ=0)
 TTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTG                                                                                                                                                 <  SRR19973066.1215555/150‑1 (MQ=32)
   GGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCC                                                                                                                                               <  SRR19973066.2153813/150‑1 (MQ=32)
   GGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCC                                                                                                                                               <  SRR19973066.2153868/150‑1 (MQ=32)
      TATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCT                                                                                                                                            <  SRR19973066.446749/150‑1 (MQ=32)
         AGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCA                                                                                                                                         <  SRR19973066.2490406/150‑1 (MQ=32)
         AGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCA                                                                                                                                         <  SRR19973066.851892/150‑1 (MQ=32)
                   ACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTAGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAA                                                                                                                               >  SRR19973066.3534414/1‑150 (MQ=31)
                     AGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACC                                                                                                                             >  SRR19973066.2082919/1‑150 (MQ=30)
                     AGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGTGTCGTTGACAGTCATTCATTTTTCTGCCCCTCCAAAACCAAAAACC                                                                                                                             >  SRR19973066.2082898/1‑150 (MQ=29)
                               CTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGGCCTTCATTTTTCTGCCCCGCCCAAAGCAAAAACCCGCCGACGCG                                                                                                                   >  SRR19973066.3140182/1‑150 (MQ=18)
                                     GCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTT                                                                                                             >  SRR19973066.280827/1‑150 (MQ=32)
                                     GCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTT                                                                                                             <  SRR19973066.2400931/150‑1 (MQ=32)
                                     GCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTT                                                                                                             <  SRR19973066.600169/150‑1 (MQ=32)
                                      CCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTT                                                                                                            <  SRR19973066.3393747/150‑1 (MQ=7)
                                      CCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTT                                                                                                            <  SRR19973066.3393754/150‑1 (MQ=32)
                                       CCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGGCCCTGCAAAAGCAAAAACCCGCCGCAGCGGGGTTTTA                                                                                                           >  SRR19973066.1972538/1‑150 (MQ=28)
                                         CTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAGTCATTTTTATGCCCCACCAAAAGCAAAAACCCGCCGAACCGGGTTTTTACG                                                                                                         >  SRR19973066.1859840/1‑150 (MQ=28)
                                              AAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGGCATTCATTTTTCTGCCCCTCCAAAAACAAAAAacccccccagcgggtttttaacgaaac                                                                                                    >  SRR19973066.3074929/1‑123 (MQ=60)
                                              AAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC                                                                                                    <  SRR19973066.419525/150‑1 (MQ=34)
                                              AAACTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAC                                                                                                    <  SRR19973066.419646/150‑1 (MQ=25)
                                                       ACCCGCCGTAACGGGTTTTTACGTACAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA                                                                                           <  SRR19973066.2258092/150‑1 (MQ=6)
                                                       ACCCGCCGAAGCGGGTTTTTACGTACAACCGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAA                                                                                           <  SRR19973066.451546/150‑1 (MQ=29)
                                                          CGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTG                                                                                        <  SRR19973066.2060043/150‑1 (MQ=32)
                                                                  CGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTATGCCCCGCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTTCGTAAACCAGGTGAAACTGACCGATAA                                                                                >  SRR19973066.1324542/1‑150 (MQ=57)
                                                                    GGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCGCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGTAAACCAGGTGAAACTGACCGATAAGC                                                                              >  SRR19973066.242188/1‑150 (MQ=30)
                                                                            CGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCT                                                                      >  SRR19973066.3200579/1‑150 (MQ=19)
                                                                            CGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAAATGACCGATAAGCCGCTTTCT                                                                      >  SRR19973066.2669126/1‑150 (MQ=24)
                                                                              TAAAACAGGTGAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTT                                                                    <  SRR19973066.3446060/150‑1 (MQ=41)
                                                                              TAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTT                                                                    <  SRR19973066.645904/150‑1 (MQ=42)
                                                                               AAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTT                                                                   <  SRR19973066.1447425/150‑1 (MQ=40)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTTCGGAAAACAGGTGAAACTGACCGGTAAGCCGCggtctgttg                                                                  >  SRR19973066.3433819/1‑141 (MQ=14)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTG                                                                  <  SRR19973066.357796/150‑1 (MQ=38)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAATCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAACCACGTGAAACTGACCGCTAAGCCGggttgtgtgg                                                                  >  SRR19973066.3433700/1‑140 (MQ=16)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTG                                                                  <  SRR19973066.2162316/150‑1 (MQ=37)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCACCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAATCCGCGTTCTTTTG                                                                  >  SRR19973066.3433796/1‑150 (MQ=31)
                                                                                AAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTAGAGTGTCGTGGACAGTCATTCATTTTTATGCCCCTCCAAAATCAAAAAACCGCCGAAGCGGGTTTTTAAGTAAAAAAGGTGTAACTGAACGATAAGCCGCTgtgttttt                                                                  >  SRR19973066.3407927/1‑142 (MQ=13)
                                                                                          AACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTGTTCCCCGCCAAAAGCAAAAAACCCCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAAAATCCGCTTTCTTTTTGGTTTAGCGG                                                        >  SRR19973066.2717767/1‑150 (MQ=16)
                                                                                            CTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCG                                                      <  SRR19973066.678465/150‑1 (MQ=32)
                                                                                            CTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCG                                                      <  SRR19973066.1052341/150‑1 (MQ=32)
                                                                                              GACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGTCAGTCATTAAATTTTGTTCCCCACCAAAAACAAAAAAACCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAACAGACCCATAATCCGCTTTGTTTTGGGTGTAGCGTCGTG                                                    >  SRR19973066.1244388/1‑150 (MQ=3)
                                                                                                 CGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAC                                                 <  SRR19973066.708192/150‑1 (MQ=32)
                                                                                                   ATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTATGCCCCACCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGAACGATAAGCCGCTTTCTTTTTGGTATAGCGTCGTGGCCAG                                               >  SRR19973066.277864/1‑150 (MQ=27)
                                                                                                    TAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGT                                              <  SRR19973066.610030/150‑1 (MQ=32)
                                                                                                     AAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGTCAATAATTTTTATGCCCCGCCAAAACCAAAAAACCCCCGAAGCGGGGTTTTACGTAAAACAGGTGAAACTtaacgataaaccgcggtatttttgggatagagtggttggaaagc                                             >  SRR19973066.23220/1‑106 (MQ=0)
                                                                                                     AAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTC                                             <  SRR19973066.1713927/150‑1 (MQ=32)
                                                                                                      AGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTATGCCCCCCCAAAAACAAAAACCCCCCCAAACCGGTTTTTACGTAAACCACGTGAAAAAGACCGATAATACCCGTTATTTTGGGTATAGCGTCCTGGTCAGTCA                                            >  SRR19973066.2838469/1‑150 (MQ=14)
                                                                                                       GCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCCCCCAACCGGGTTTTTAAGTAAAACAGGTGAAAAAGACCGATAACCCGCTTTCTTTTTGGTATATCGTCGgggagaggtaa                                           >  SRR19973066.3354391/1‑139 (MQ=0)
                                                                                                       GCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAT                                           <  SRR19973066.2353679/150‑1 (MQ=32)
                                                                                                             TTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTtgaaagtcaatcaaaa                                     >  SRR19973066.3000163/1‑134 (MQ=7)
                                                                                                                TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTC                                  <  SRR19973066.2059134/150‑1 (MQ=32)
                                                                                                                TTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTC                                  <  SRR19973066.2059220/150‑1 (MQ=32)
                                                                                                                 ccTGGGTATAGTCACGTGGACAGTCATTCATATTACTGCCCCAACAAAATCAAAAACCCGGCGTATCCGCTTTTTTCTTAAACCAGCTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                 <  SRR19973066.3248904/148‑1 (MQ=16)
                                                                                                                 TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                 <  SRR19973066.2688631/150‑1 (MQ=32)
                                                                                                                 TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                 <  SRR19973066.2124424/150‑1 (MQ=32)
                                                                                                                 TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTCCCCCTCCACAAGCAAAAACCCGCCGAATCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                 <  SRR19973066.888306/150‑1 (MQ=29)
                                                                                                                 TTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCGGCCCCTCAAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCT                                 <  SRR19973066.2124414/150‑1 (MQ=30)
                                                                                                                      GTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAACAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAATCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCAATCAAATTTGTGCCCC                            >  SRR19973066.61583/1‑150 (MQ=5)
                                                                                                                        ATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGCGTCGTGGGCAGTCAatcatattttttccccgc                          >  SRR19973066.1634155/1‑132 (MQ=7)
                                                                                                                        ATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAACCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGgaactcaataaaaattgttcccccc                          >  SRR19973066.1385293/1‑125 (MQ=33)
                                                                                                                         gcGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCC                         <  SRR19973066.333813/148‑1 (MQ=32)
                                                                                                                          AGTGTCGTGTACAGTCATTCATGTTTCTGCCCCTCCAAAACCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCA                        <  SRR19973066.2117465/150‑1 (MQ=30)
                                                                                                                          AGTGTCGTGGACAGTCATTCATTTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTACACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCA                        <  SRR19973066.1041619/150‑1 (MQ=30)
                                                                                                                          AGTGTCGTGGACAGTCATTCATTTGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCA                        <  SRR19973066.2458063/150‑1 (MQ=31)
                                                                                                                             GTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGGCAGTCctgcctcgttctggcccgcccaca                     >  SRR19973066.1763820/1‑126 (MQ=32)
                                                                                                                                 TGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAA                 >  SRR19973066.889362/1‑150 (MQ=32)
                                                                                                                                  GGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTATGCCCCACCAAAAGCAAA                >  SRR19973066.3455015/1‑150 (MQ=30)
                                                                                                                                    ACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCATTCATCTTTCTGCCCCGCCAAAAGCAAAAA              >  SRR19973066.3200336/1‑150 (MQ=30)
                                                                                                                                     CAGTCATTAATTTTTCTGCCACTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAC             <  SRR19973066.3181455/150‑1 (MQ=7)
                                                                                                                                       GTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCCTCGTGGACAGTaataaaactttatgccccgcaaaaagcaaaaaacc           >  SRR19973066.1186722/1‑115 (MQ=30)
                                                                                                                                        TCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCG          <  SRR19973066.1817375/150‑1 (MQ=31)
                                                                                                                                          ATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATAAAATTTCTGCCCCACCAAAATCAAAAACCCGCC        >  SRR19973066.840010/1‑150 (MQ=27)
                                                                                                                                          ATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGTCCGGTAAGCCGCGTTGTTTTGGGGTTAGCGTCGTtggaagtgattaatgattctggcccacccaaagcaaaaaaccccc        >  SRR19973066.3590599/1‑105 (MQ=10)
                                                                                                                                           TTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCG       <  SRR19973066.2290591/150‑1 (MQ=32)
                                                                                                                                              ATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTTGGTTTAGCGTCGTGGAaatttattcaaccttgatccccaacaatatcaaaaaccccccgacg    >  SRR19973066.2413245/1‑104 (MQ=48)
                                                                                                                                              ATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTGATCTTTCTGCCCCGCCAAAAGCACAAACCCCCCGAAG    >  SRR19973066.2413527/1‑150 (MQ=29)
                                                                                                                                                TTTTCTGCCCATCCAAATGCAAAATTCCGCCGAAGGGTGTTTTTACGTAAAGCTGGTGAATCTGATGGATATGCCGCTTTGTTTTGGGTTTAGCGTCGTGGTCTGTCATTCTTCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCG  <  SRR19973066.791603/150‑1 (MQ=13)
                                                                                                                                                CTGTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATGTGGTGTTACTGAACGATATGCGGGTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCG  <  SRR19973066.791636/150‑1 (MQ=0)
                                                                                                                                                |                                                                                                                                                     
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCG  >  CP009273/3263265‑3263558

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