Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,229,926 C→A 20.2% L10I (CTC→ATC)  ygjQ → DUF218 superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,229,9260CA20.2% 873.8 / 143.3 304L10I (CTC→ATC) ygjQDUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (159/83);  new base A (28/34);  total (187/117)
Fisher's exact test for biased strand distribution p-value = 3.53e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

GACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAAC  >  CP009273/3229779‑3230075
                                                                                                                                                   |                                                                                                                                                     
 aCTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTcc                                                                                                                                                    <  2:4527628/150‑1 (MQ=255)
   tgtgGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGt                                                                                                                                                  >  1:3960343/1‑150 (MQ=255)
   tgtgGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGt                                                                                                                                                  >  1:1782647/1‑150 (MQ=255)
     tgGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTAt                                                                                                                                                >  2:4502996/1‑150 (MQ=255)
       gCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTACTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTAttt                                                                                                                                              <  1:2056040/150‑1 (MQ=255)
           aCGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGttt                                                                                                                                          <  2:1966361/150‑1 (MQ=255)
            cGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGtttt                                                                                                                                         <  1:2299551/150‑1 (MQ=255)
             gCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGttttt                                                                                                                                        >  1:2257502/1‑150 (MQ=255)
               aGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTtct                                                                                                                                      >  1:2013536/1‑150 (MQ=255)
               aGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTtct                                                                                                                                      >  2:3969042/1‑150 (MQ=255)
                gTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTtctc                                                                                                                                     <  1:2817356/150‑1 (MQ=255)
                 ttGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCg                                                                                                                                    <  1:1333225/150‑1 (MQ=255)
                 ttGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCg                                                                                                                                    >  1:2336759/1‑150 (MQ=255)
                 ttGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCg                                                                                                                                    <  1:1333058/150‑1 (MQ=255)
                   gTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGcc                                                                                                                                  <  1:2526226/150‑1 (MQ=255)
                   gTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGcc                                                                                                                                  >  2:4314309/1‑150 (MQ=255)
                   gTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGcc                                                                                                                                  >  2:245997/1‑150 (MQ=255)
                      gCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGa                                                                                                                               <  1:3658393/150‑1 (MQ=255)
                        ttGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGAct                                                                                                                             <  1:1879272/150‑1 (MQ=255)
                        ttGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGAct                                                                                                                             <  2:2265205/150‑1 (MQ=255)
                          gggCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGActct                                                                                                                           <  1:3686527/150‑1 (MQ=255)
                          gggCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGActct                                                                                                                           >  2:3629671/1‑150 (MQ=255)
                          gggCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGActct                                                                                                                           >  2:2681644/1‑150 (MQ=255)
                            gCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTaa                                                                                                                         <  1:402924/150‑1 (MQ=255)
                            gCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTaa                                                                                                                         >  2:742161/1‑150 (MQ=255)
                              aGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAg                                                                                                                       >  2:1373161/1‑150 (MQ=255)
                               gAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGa                                                                                                                      <  1:101724/150‑1 (MQ=255)
                                 aCAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGTATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGAtt                                                                                                                    <  1:2564699/150‑1 (MQ=255)
                                 aCAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGAtt                                                                                                                    <  2:2816719/150‑1 (MQ=255)
                                   aaaaTCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGc                                                                                                                  >  2:4130045/1‑150 (MQ=255)
                                    aaaTCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGcc                                                                                                                 >  2:1557225/1‑150 (MQ=255)
                                    aaaTCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGcc                                                                                                                 >  2:794360/1‑150 (MQ=255)
                                    aaaTCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGcc                                                                                                                 <  2:3885568/150‑1 (MQ=255)
                                     aaTCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCt                                                                                                                >  2:439507/1‑150 (MQ=255)
                                          gAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTa                                                                                                           <  1:1100662/150‑1 (MQ=255)
                                           aaaaTGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTAc                                                                                                          >  2:3237965/1‑150 (MQ=255)
                                           aaaaTGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTAc                                                                                                          >  2:183514/1‑150 (MQ=255)
                                           aaaaTGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTAc                                                                                                          <  1:702131/150‑1 (MQ=255)
                                                tatgTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTa                                                                                                     >  1:4636178/2‑150 (MQ=255)
                                                gatgTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTa                                                                                                     <  1:4461241/150‑1 (MQ=255)
                                                     cAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTtgt                                                                                                >  1:1234583/1‑150 (MQ=255)
                                                          gTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAg                                                                                           >  1:3882268/1‑150 (MQ=255)
                                                          gTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCGTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAg                                                                                           <  1:162594/150‑1 (MQ=255)
                                                           ttATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGGATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAgg                                                                                          >  2:3504109/1‑150 (MQ=255)
                                                            tATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGa                                                                                         <  2:3340714/150‑1 (MQ=255)
                                                            tATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGa                                                                                         <  2:1017165/150‑1 (MQ=255)
                                                              tcttttctttCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGc                                                                                       <  2:4079722/150‑1 (MQ=255)
                                                                 tttctttcTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGAc                                                                                    >  1:3499843/1‑150 (MQ=255)
                                                                   tctttcttTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCa                                                                                  <  2:3960343/150‑1 (MQ=255)
                                                                   tctttcttTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCa                                                                                  >  1:4314381/1‑150 (MQ=255)
                                                                         tttAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTa                                                                            >  2:1408928/1‑150 (MQ=255)
                                                                         tttAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTa                                                                            >  2:65324/1‑150 (MQ=255)
                                                                           tAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATc                                                                          >  1:3158496/1‑150 (MQ=255)
                                                                            aGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCg                                                                         >  2:2773585/1‑150 (MQ=255)
                                                                               gcCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCg                                                                      >  1:3770179/1‑150 (MQ=255)
                                                                                ccTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGa                                                                     >  2:4064142/1‑150 (MQ=255)
                                                                                ccTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGa                                                                     >  2:612384/1‑150 (MQ=255)
                                                                                 cTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGAt                                                                    >  1:872187/1‑150 (MQ=255)
                                                                                  tAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATc                                                                   >  2:327840/1‑150 (MQ=255)
                                                                                  tAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATAGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATc                                                                   >  2:2645522/1‑150 (MQ=255)
                                                                                  tAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATc                                                                   <  1:95984/150‑1 (MQ=255)
                                                                                   aaaaTCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATcg                                                                  <  1:3136136/150‑1 (MQ=255)
                                                                                    aaaTCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATcgc                                                                 >  2:1879499/1‑150 (MQ=255)
                                                                                     aaTCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATcgcg                                                                <  1:2181947/150‑1 (MQ=255)
                                                                                     aaTCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATcgcg                                                                >  2:3098875/1‑150 (MQ=255)
                                                                                          aCCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATg                                                           >  2:4625958/1‑150 (MQ=255)
                                                                                           ccTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATgg                                                          <  2:1460563/150‑1 (MQ=255)
                                                                                             tCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTc                                                        <  1:1673726/150‑1 (MQ=255)
                                                                                                   cTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGcc                                                  >  2:3078284/1‑150 (MQ=255)
                                                                                                    tttCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCACCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGccc                                                 <  1:4206376/150‑2 (MQ=255)
                                                                                                            aGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAg                                         <  2:3172102/150‑1 (MQ=255)
                                                                                                             ggATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGc                                        <  1:3791299/150‑1 (MQ=255)
                                                                                                             ggATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGc                                        >  1:565818/1‑150 (MQ=255)
                                                                                                             ggATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGc                                        >  1:2303330/1‑150 (MQ=255)
                                                                                                               aTAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTg                                      >  2:2973231/1‑150 (MQ=255)
                                                                                                               aTAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTg                                      <  1:4541181/150‑1 (MQ=255)
                                                                                                                 aaCGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGAc                                    <  1:1498058/150‑1 (MQ=255)
                                                                                                                 aaCGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGAc                                    <  1:3266237/150‑1 (MQ=255)
                                                                                                                       cATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAg                              >  2:3717408/1‑150 (MQ=255)
                                                                                                                       cATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAg                              >  2:1932674/1‑150 (MQ=255)
                                                                                                                        aTGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCATATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGc                             >  2:1828570/1‑150 (MQ=255)
                                                                                                                           ctgcGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGAt                          <  2:3018933/150‑1 (MQ=255)
                                                                                                                            tgcGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATg                         >  2:1696922/1‑150 (MQ=255)
                                                                                                                             gcgcgcATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGt                        <  1:2636679/150‑1 (MQ=255)
                                                                                                                              cgcgcATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTc                       <  1:1301771/150‑1 (MQ=255)
                                                                                                                               gcgcATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCa                      <  1:4370222/150‑1 (MQ=255)
                                                                                                                               gcgcATTTGCCCGCCTTCTTATCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCa                      >  2:4221702/1‑150 (MQ=255)
                                                                                                                                cgcATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCaa                     <  2:2883687/150‑1 (MQ=255)
                                                                                                                                 gcATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAAc                    <  1:1935540/150‑1 (MQ=255)
                                                                                                                                 gcATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAAc                    >  2:1175722/1‑150 (MQ=255)
                                                                                                                                   aTTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGc                  >  2:3375793/1‑150 (MQ=255)
                                                                                                                                      tGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGt               >  1:3588938/1‑150 (MQ=255)
                                                                                                                                        cccGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTc             >  2:754024/1‑150 (MQ=255)
                                                                                                                                        cccGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTc             >  2:2888895/1‑150 (MQ=255)
                                                                                                                                          cGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCg           >  1:2754180/1‑150 (MQ=255)
                                                                                                                                              ttcttaTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGgcgc       >  1:3500441/1‑150 (MQ=255)
                                                                                                                                                cttctCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGt     <  2:1747834/150‑1 (MQ=255)
                                                                                                                                                 ttctCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTa    >  2:2129983/1‑150 (MQ=255)
                                                                                                                                                 ttctCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCATGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTa    <  2:1201955/150‑1 (MQ=255)
                                                                                                                                                   ctcCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAAc  >  2:1737238/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACTGTGGCTAACGCAGTTGTCGCTTGGGCAGAACAAAATCTGAAAATGATGTCAGGCGTTATCTTTTCTTTCTTTAGCGCCTAAAATCGACCTCCCCCCTTTCGTTAAGGATAACGACCATGCTGCGCGCATTTGCCCGCCTTCTTCTCCGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAAC  >  CP009273/3229779‑3230075

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: ‑