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 868,500 A→T 100% S22S (TCA→TCT)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273868,5000AT97.0% 730.5 / 9.0 235S22S (TCA→TCT) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base T (110/118);  total (112/123)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 240 total aligned reads displayed.

ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGA  >  CP009273/868355‑868643
                                                                                                                                                 |                                                                                                                                               
 tttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGATGATAATTTTATTTTTATTGGCAGCACTGTCtttt                                                                                                                                              >  1:2216862/1‑148 (MQ=255)
  ttGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTa                                                                                                                                             >  2:798438/1‑148 (MQ=255)
   tGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                            <  1:3472221/148‑1 (MQ=255)
   tGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAt                                                                                                                                            <  1:1453855/148‑1 (MQ=255)
     gTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTg                                                                                                                                          >  1:1021766/1‑148 (MQ=255)
       gaacaaGCCGCAGCAAGGGACAATTTTGGCCGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAAGATGCTGTGTTTATAAGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGt                                                                                                                                        <  2:1556079/144‑1 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                  >  1:730729/1‑148 (MQ=255)
             gCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttt                                                                                                                                  <  2:113844/148‑1 (MQ=255)
               cGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttttt                                                                                                                                >  2:3306899/1‑148 (MQ=255)
               cGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCttttttt                                                                                                                                >  2:120431/1‑148 (MQ=255)
                      aGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCAtt                                                                                                                         >  2:3553792/1‑148 (MQ=255)
                       gggACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTa                                                                                                                        >  2:937845/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    >  2:1312971/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    >  2:1314400/1‑148 (MQ=255)
                           cAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTa                                                                                                                    <  1:2456095/148‑1 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                >  1:567399/1‑148 (MQ=255)
                               ttttGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCaa                                                                                                                <  2:221828/148‑1 (MQ=255)
                                  tGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTg                                                                                                             <  1:2954585/148‑1 (MQ=255)
                                    cAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGt                                                                                                           <  2:1830168/148‑1 (MQ=255)
                                     aGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTa                                                                                                          <  1:974921/148‑1 (MQ=255)
                                        ggCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCg                                                                                                       <  1:3153427/148‑1 (MQ=255)
                                          cacaGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGa                                                                                                     <  2:1205287/148‑1 (MQ=255)
                                              gCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAgcgc                                                                                                 >  1:2916354/1‑148 (MQ=255)
                                                 ttCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                              <  1:777822/148‑1 (MQ=255)
                                                 ttCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGcgcgcg                                                                                              >  1:1590425/1‑148 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                         <  2:2516960/148‑1 (MQ=255)
                                                      aTAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTaaa                                                                                         <  1:1188968/148‑1 (MQ=255)
                                                         gTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACg                                                                                      >  2:893209/1‑148 (MQ=255)
                                                          ttATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGt                                                                                     >  1:1367846/1‑148 (MQ=255)
                                                            aTTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGc                                                                                   >  1:3213943/1‑148 (MQ=255)
                                                             ttttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGcc                                                                                  >  1:81568/1‑148 (MQ=255)
                                                               ttGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGa                                                                                <  1:1709708/148‑1 (MQ=255)
                                                                tGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAt                                                                               <  2:2901903/148‑1 (MQ=255)
                                                                tGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAATATCAACTGGTATCGGAGCGCGCGGTCAAACGTGCCGAt                                                                               >  1:2460009/1‑148 (MQ=255)
                                                                 gTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATa                                                                              >  2:1476323/1‑148 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAg                                                                             <  1:1892411/148‑1 (MQ=255)
                                                                  ttAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAg                                                                             >  1:2889390/1‑148 (MQ=255)
                                                                     aaTGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCg                                                                          >  1:533002/1‑148 (MQ=255)
                                                                       tGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATATTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCt                                                                        >  2:2166534/1‑148 (MQ=255)
                                                                         tATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCttt                                                                      >  1:1717298/1‑148 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                     >  2:1282010/1‑148 (MQ=255)
                                                                          aTTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTg                                                                     >  1:106238/1‑148 (MQ=255)
                                                                           ttAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGa                                                                    >  1:3076834/1‑148 (MQ=255)
                                                                            tAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGaa                                                                   >  2:2027157/1‑148 (MQ=255)
                                                                             aaCATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAAc                                                                  <  1:471733/148‑1 (MQ=255)
                                                                               cATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACtt                                                                >  2:1533832/1‑148 (MQ=255)
                                                                                 tGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTAt                                                              <  2:656016/148‑1 (MQ=255)
                                                                                  gCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTAATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAAAGCCGCTTTGAACTAAtt                                                             >  1:338070/1‑148 (MQ=255)
                                                                                     gAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAg                                                          >  2:344726/1‑148 (MQ=255)
                                                                                      aGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGa                                                         >  2:837463/1‑148 (MQ=255)
                                                                                       gTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaa                                                        >  1:2162500/1‑148 (MQ=255)
                                                                                         ttATACTAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGaaaa                                                      <  2:1885925/148‑1 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACg                                                   >  2:2198387/1‑148 (MQ=255)
                                                                                            taCGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACg                                                   >  2:404395/1‑148 (MQ=255)
                                                                                              cGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGtt                                                 <  1:1840093/148‑1 (MQ=255)
                                                                                               gAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTg                                                <  1:1860889/148‑1 (MQ=255)
                                                                                                 aaagataaagatAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGc                                              <  1:1122807/148‑1 (MQ=255)
                                                                                                  aagataaagataAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCt                                             >  1:3144704/1‑148 (MQ=255)
                                                                                                   agataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTa                                            <  2:3005513/148‑1 (MQ=255)
                                                                                                    gataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAt                                           <  2:3223268/148‑1 (MQ=255)
                                                                                                     ataaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAtt                                          >  1:2994326/1‑148 (MQ=255)
                                                                                                      taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                         <  1:637937/148‑1 (MQ=255)
                                                                                                      taaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAttt                                         <  2:2083485/148‑1 (MQ=255)
                                                                                                       aaagataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTc                                        >  1:2030982/1‑148 (MQ=255)
                                                                                                          gataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCttt                                     <  1:62634/148‑1 (MQ=255)
                                                                                                           ataaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTa                                    >  1:2287969/1‑148 (MQ=255)
                                                                                                            taaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaa                                   <  2:3170443/148‑1 (MQ=255)
                                                                                                             aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                  >  1:2367348/1‑148 (MQ=255)
                                                                                                             aaGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTaaa                                  >  1:1647950/1‑148 (MQ=255)
                                                                                                              aGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAg                                 <  2:2268151/148‑1 (MQ=255)
                                                                                                               ggCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGa                                >  1:1765994/1‑148 (MQ=255)
                                                                                                                gCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGat                               >  1:3540175/1‑148 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  2:1021766/148‑1 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  1:3306899/148‑1 (MQ=255)
                                                                                                                 cTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGata                              <  1:120431/148‑1 (MQ=255)
                                                                                                                  ttaTAATTTTATTTTTATTGGCATCACTTTATTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGatat                             <  2:2324992/146‑1 (MQ=255)
                                                                                                                     tAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGa                          <  1:2669802/148‑1 (MQ=255)
                                                                                                                       attttattttTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAAc                        >  2:2379116/1‑148 (MQ=255)
                                                                                                                        ttttatttttatTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACg                       >  2:2426132/1‑148 (MQ=255)
                                                                                                                         tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCTATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                      <  1:3112332/148‑1 (MQ=255)
                                                                                                                         tttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGt                      <  2:1910495/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  1:2857273/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  1:1394989/148‑1 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     >  2:3091871/1‑148 (MQ=255)
                                                                                                                          ttatttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGtt                     <  2:1712125/148‑1 (MQ=255)
                                                                                                                            atttttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCg                   >  2:777222/1‑148 (MQ=255)
                                                                                                                               tttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGcc                <  1:174926/148‑1 (MQ=255)
                                                                                                                                ttattGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGccc               <  2:533002/148‑1 (MQ=255)
                                                                                                                                         gcaCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaa      <  1:1994777/148‑1 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa     <  2:1690970/148‑1 (MQ=255)
                                                                                                                                          caCTGTCTTTTATTGGTCTTTATTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTaaa     <  2:2232590/148‑1 (MQ=255)
                                                                                                                                           aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg    <  2:1845730/148‑1 (MQ=255)
                                                                                                                                           aCTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAg    <  2:567399/148‑1 (MQ=255)
                                                                                                                                            cTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTTGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAgg   >  1:2234389/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGTCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  >  2:607959/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGTCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  >  2:606356/1‑148 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  1:1554632/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:1964287/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:3365403/148‑1 (MQ=255)
                                                                                                                                             tGTCTTTTATTGGTCTTTATTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGa  <  2:2682745/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGA  >  CP009273/868355‑868643

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
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCA  >  CP009273/868354‑868647
                                                                                                                                                  |                                                                                                                                                   
AATTTTTTTAAAAAGCCGCATCAAGGGACAATTTTTTCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTT                                                                                                                                                  <  SRR19973066.2897949/150‑1 (MQ=60)
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAACGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTT                                                                                                                                                  <  SRR19973066.1108414/150‑1 (MQ=60)
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTT                                                                                                                                                  <  SRR19973066.1801676/150‑1 (MQ=60)
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTT                                                                                                                                                  >  SRR19973066.204533/1‑150 (MQ=60)
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGATGATAATTTTATTTTTATTGGCAGCACTGTCTTTT                                                                                                                                                  >  SRR19973066.2232206/1‑150 (MQ=60)
 ATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTA                                                                                                                                                 >  SRR19973066.1858933/1‑150 (MQ=60)
  TTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTAT                                                                                                                                                >  SRR19973066.3134927/1‑150 (MQ=60)
    TGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTG                                                                                                                                              >  SRR19973066.1029283/1‑150 (MQ=60)
    TGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTG                                                                                                                                              <  SRR19973066.3495409/150‑1 (MQ=60)
    TGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTG                                                                                                                                              <  SRR19973066.1464164/150‑1 (MQ=60)
        TAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCT                                                                                                                                          >  SRR19973066.3244587/1‑150 (MQ=60)
        TAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCT                                                                                                                                          >  SRR19973066.3244495/1‑150 (MQ=60)
            AAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTT                                                                                                                                      >  SRR19973066.736150/1‑150 (MQ=60)
                CGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCA                                                                                                                                  <  SRR19973066.1190644/150‑1 (MQ=60)
                  CAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATC                                                                                                                                >  SRR19973066.3000301/1‑150 (MQ=60)
                     CAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATT                                                                                                                             <  SRR19973066.1765659/150‑1 (MQ=60)
                        GGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAC                                                                                                                          >  SRR19973066.3546088/1‑150 (MQ=60)
                           ACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTAT                                                                                                                       <  SRR19973066.1815812/150‑1 (MQ=60)
                            CAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATC                                                                                                                      <  SRR19973066.2472834/150‑1 (MQ=60)
                            CAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATC                                                                                                                      <  SRR19973066.2472905/150‑1 (MQ=60)
                              ATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAA                                                                                                                    >  SRR19973066.571418/1‑150 (MQ=60)
                                   TGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGT                                                                                                               <  SRR19973066.2974175/150‑1 (MQ=60)
                                      AGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATC                                                                                                            <  SRR19973066.982094/150‑1 (MQ=60)
                                         GGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGA                                                                                                         <  SRR19973066.3174453/150‑1 (MQ=60)
                                             CAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGC                                                                                                     >  SRR19973066.2935810/1‑150 (MQ=60)
                                             CAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGC                                                                                                     >  SRR19973066.1555119/1‑150 (MQ=60)
                                             CAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGC                                                                                                     >  SRR19973066.1146867/1‑150 (MQ=60)
                                                CGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCG                                                                                                  >  SRR19973066.1601933/1‑150 (MQ=60)
                                                  TTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGT                                                                                                <  SRR19973066.783817/150‑1 (MQ=60)
                                                  TTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGT                                                                                                <  SRR19973066.783571/150‑1 (MQ=60)
                                                   TCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAATAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTA                                                                                               <  SRR19973066.1597533/150‑1 (MQ=60)
                                                    CAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAA                                                                                              <  SRR19973066.2920187/150‑1 (MQ=60)
                                                     AGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAA                                                                                             >  SRR19973066.3451834/1‑150 (MQ=60)
                                                      GATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAA                                                                                            >  SRR19973066.3169887/1‑150 (MQ=60)
                                                       ATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAAC                                                                                           <  SRR19973066.1252337/150‑1 (MQ=60)
                                                       ATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAAC                                                                                           <  SRR19973066.1197304/150‑1 (MQ=60)
                                                         AGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGT                                                                                         >  SRR19973066.1377741/1‑150 (MQ=60)
                                                         AGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGT                                                                                         >  SRR19973066.611993/1‑150 (MQ=60)
                                                           TTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGC                                                                                       >  SRR19973066.3235229/1‑150 (MQ=60)
                                                            TATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCC                                                                                      >  SRR19973066.82493/1‑150 (MQ=60)
                                                            TATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCC                                                                                      <  SRR19973066.1653230/150‑1 (MQ=60)
                                                            TATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCC                                                                                      <  SRR19973066.1652730/150‑1 (MQ=60)
                                                               TTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAAATATCAACTGGTATCGGAGCGCGCGGTCAAACGTGCCGAT                                                                                   >  SRR19973066.2476775/1‑150 (MQ=60)
                                                                TTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATA                                                                                  <  SRR19973066.1722195/150‑1 (MQ=60)
                                                                TTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATA                                                                                  <  SRR19973066.1722110/150‑1 (MQ=60)
                                                                 TGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAG                                                                                 >  SRR19973066.2908713/1‑150 (MQ=60)
                                                                  GTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGC                                                                                <  SRR19973066.951432/150‑1 (MQ=60)
                                                                   TTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCC                                                                               <  SRR19973066.1905903/150‑1 (MQ=60)
                                                                    TAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCG                                                                              >  SRR19973066.536888/1‑150 (MQ=60)
                                                                      AATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCT                                                                            <  SRR19973066.2405666/150‑1 (MQ=60)
                                                                       ATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTT                                                                           >  SRR19973066.1058946/1‑150 (MQ=60)
                                                                        TGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTT                                                                          >  SRR19973066.1729816/1‑150 (MQ=60)
                                                                        TGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTT                                                                          >  SRR19973066.1729745/1‑150 (MQ=60)
                                                                         GTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTG                                                                         >  SRR19973066.107339/1‑150 (MQ=60)
                                                                          TATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGA                                                                        >  SRR19973066.3097386/1‑150 (MQ=60)
                                                                           ATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAA                                                                       >  SRR19973066.1006682/1‑150 (MQ=60)
                                                                              AACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTT                                                                    <  SRR19973066.475309/150‑1 (MQ=60)
                                                                                 ATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTAATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGAAAGCCGCTTTGAACTAATT                                                                 >  SRR19973066.340672/1‑150 (MQ=60)
                                                                                   GCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCA                                                               >  SRR19973066.78433/1‑150 (MQ=60)
                                                                                      GAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAA                                                            >  SRR19973066.2177526/1‑150 (MQ=60)
                                                                                        GTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAA                                                          >  SRR19973066.1098160/1‑150 (MQ=60)
                                                                                        GTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAA                                                          >  SRR19973066.3129609/1‑150 (MQ=60)
                                                                                          TTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAC                                                        <  SRR19973066.334228/150‑1 (MQ=60)
                                                                                               CGAAAAGATAAAGATCAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGG                                                   <  SRR19973066.1853211/150‑1 (MQ=60)
                                                                                               CGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGG                                                   <  SRR19973066.1853241/150‑1 (MQ=60)
                                                                                               CGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGG                                                   <  SRR19973066.2225258/150‑1 (MQ=60)
                                                                                               CGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTATTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGG                                                   >  SRR19973066.2248000/1‑150 (MQ=60)
                                                                                                GAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGC                                                  <  SRR19973066.1874195/150‑1 (MQ=60)
                                                                                                 AAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCT                                                 >  SRR19973066.421746/1‑150 (MQ=60)
                                                                                                 AAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCT                                                 >  SRR19973066.3165679/1‑150 (MQ=60)
                                                                                                 AAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGGTGGCT                                                 >  SRR19973066.485571/1‑150 (MQ=60)
                                                                                                  AAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTA                                                <  SRR19973066.2920784/150‑1 (MQ=60)
                                                                                                  AAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTA                                                <  SRR19973066.1130897/150‑1 (MQ=60)
                                                                                                   AAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTAT                                               >  SRR19973066.2267278/1‑150 (MQ=60)
                                                                                                    AGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCCACTGGTATCGGAGCTCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATT                                              >  SRR19973066.3017214/1‑150 (MQ=60)
                                                                                                    AGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATT                                              >  SRR19973066.2962102/1‑150 (MQ=60)
                                                                                                    AGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATT                                              >  SRR19973066.3014179/1‑150 (MQ=60)
                                                                                                      ATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTC                                            >  SRR19973066.2046132/1‑150 (MQ=60)
                                                                                                      ANAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTC                                            >  SRR19973066.2045163/1‑150 (MQ=60)
                                                                                                       TAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCT                                           <  SRR19973066.642442/150‑1 (MQ=60)
                                                                                                       TAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCT                                           <  SRR19973066.642539/150‑1 (MQ=60)
                                                                                                       TAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCT                                           >  SRR19973066.17286/1‑150 (MQ=60)
                                                                                                          AGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTA                                        <  SRR19973066.2880181/150‑1 (MQ=60)
                                                                                                          AGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTA                                        <  SRR19973066.3416713/150‑1 (MQ=60)
                                                                                                          AGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTA                                        <  SRR19973066.572006/150‑1 (MQ=60)
                                                                                                          AGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTA                                        >  SRR19973066.2303666/1‑150 (MQ=60)
                                                                                                           GATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAA                                       <  SRR19973066.63403/150‑1 (MQ=60)
                                                                                                            ATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAA                                      >  SRR19973066.17516/1‑150 (MQ=60)
                                                                                                            ATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAA                                      >  SRR19973066.1659952/1‑150 (MQ=60)
                                                                                                            ATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAA                                      >  SRR19973066.2383372/1‑150 (MQ=60)
                                                                                                            ATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAA                                      >  SRR19973066.2383566/1‑150 (MQ=60)
                                                                                                              AAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGA                                    >  SRR19973066.1778659/1‑150 (MQ=60)
                                                                                                              AAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGA                                    <  SRR19973066.778281/150‑1 (MQ=60)
                                                                                                               AGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGAT                                   >  SRR19973066.3563697/1‑150 (MQ=60)
                                                                                                               AGGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGAT                                   >  SRR19973066.1502949/1‑150 (MQ=60)
                                                                                                                GGCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATA                                  <  SRR19973066.560166/150‑1 (MQ=60)
                                                                                                                 GCTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATAT                                 <  SRR19973066.2114484/150‑1 (MQ=60)
                                                                                                                  CTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATT                                <  SRR19973066.121649/150‑1 (MQ=60)
                                                                                                                  CTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATT                                <  SRR19973066.3328803/150‑1 (MQ=60)
                                                                                                                  CTGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATT                                <  SRR19973066.121670/150‑1 (MQ=60)
                                                                                                                   TGATAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTG                               <  SRR19973066.944667/150‑1 (MQ=60)
                                                                                                                      TAATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAAC                            <  SRR19973066.2687651/150‑1 (MQ=60)
                                                                                                                       AATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACG                           >  SRR19973066.708647/1‑150 (MQ=60)
                                                                                                                       AATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACG                           >  SRR19973066.2542677/1‑150 (MQ=60)
                                                                                                                        ATTTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGT                          >  SRR19973066.3066671/1‑150 (MQ=60)
                                                                                                                          TTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCTATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTC                        <  SRR19973066.3133119/150‑1 (MQ=60)
                                                                                                                          TTTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCTATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTC                        <  SRR19973066.1905590/150‑1 (MQ=60)
                                                                                                                           TTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCG                       <  SRR19973066.1405043/150‑1 (MQ=60)
                                                                                                                           TTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCG                       <  SRR19973066.1405038/150‑1 (MQ=60)
                                                                                                                           TTATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCG                       <  SRR19973066.2876466/150‑1 (MQ=60)
                                                                                                                            TATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTagatcggaagagcacacgtc                      >  SRR19973066.2869108/1‑130 (MQ=60)
                                                                                                                            TATTTTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGG                      >  SRR19973066.201018/1‑150 (MQ=60)
                                                                                                                                TTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCT                  <  SRR19973066.176473/150‑1 (MQ=60)
                                                                                                                                TTTATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCT                  <  SRR19973066.3138241/150‑1 (MQ=60)
                                                                                                                                   ATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGAC               >  SRR19973066.1211226/1‑150 (MQ=60)
                                                                                                                                   ATTGGCAGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGAC               >  SRR19973066.2738775/1‑150 (MQ=60)
                                                                                                                                         AGCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAA         >  SRR19973066.3204010/1‑150 (MQ=60)
                                                                                                                                          GCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAG        <  SRR19973066.2008803/150‑1 (MQ=60)
                                                                                                                                          GCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAG        <  SRR19973066.1438357/150‑1 (MQ=60)
                                                                                                                                          GCACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAG        <  SRR19973066.2008808/150‑1 (MQ=60)
                                                                                                                                           CACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTTGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGG       >  SRR19973066.2249863/1‑150 (MQ=60)
                                                                                                                                           CACTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTTGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGG       >  SRR19973066.2249807/1‑150 (MQ=60)
                                                                                                                                             CTGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGAC     <  SRR19973066.479495/150‑1 (MQ=60)
                                                                                                                                              TGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACT    <  SRR19973066.1565804/150‑1 (MQ=60)
                                                                                                                                              TGTCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACT    <  SRR19973066.2365110/150‑1 (MQ=60)
                                                                                                                                              TGTCTTTTATTGGTCTTTTTTTCATCATTAAATATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAAAGTTTGCTATTTCTTTAAAGATATTGAAAGTTCGGCCCTGACATTAAAGGACT    >  SRR19973066.1943993/1‑150 (MQ=60)
                                                                                                                                                TCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCA  >  SRR19973066.1780424/1‑150 (MQ=60)
                                                                                                                                                TCTTTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCA  >  SRR19973066.3463280/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                   
AATTTGGTTAAAAAGCCGCAGCAAGGGACAATTTTTGCAGCGGCACAGCGTTCAGATAGTTATTTTGTTAAATGTATTAACATGCTGAGTTTATACGAAAAGATAAAGATAAGGCTGATAATTTTATTTTTATTGGCAGCACTGTCATTTATTGGTCTTTTTTTCATCATTAACTATCAACTGGTATCGGAGCGCGCGGTAAAACGTGCCGATAGCCGCTTTGAACTTATTCAGAAAAACGTTGGCTATTTCTTTAAAGATATTGAACGTTCGGCCCTGACATTAAAGGACTCA  >  CP009273/868354‑868647

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