Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,551,503 T→C 8.6% intergenic (‑65/‑124) rtcB ← / → rtcR RNA‑splicing ligase/sigma 54‑dependent transcriptional regulator of rtcBA expression

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,551,5030TC 8.6% 665.4 / 42.6 209intergenic (‑65/‑124)rtcB/rtcRRNA‑splicing ligase/sigma 54‑dependent transcriptional regulator of rtcBA expression
Reads supporting (aligned to +/- strand):  ref base T (90/101);  new base C (2/16);  total (92/117)
Fisher's exact test for biased strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 210 total aligned reads displayed.

GCTGACGCGCATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCT  >  CP009273/3551354‑3551651
                                                                                                                                                     |                                                                                                                                                    
 cTGACGCGCATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATa                                                                                                                                                     <  1:2609201/150‑1 (MQ=255)
 cTGACGCGCATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATa                                                                                                                                                     >  1:2105052/1‑150 (MQ=255)
     cgcgCATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAagca                                                                                                                                                 <  2:1728038/150‑1 (MQ=255)
         cATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCtt                                                                                                                                             >  2:1678747/1‑150 (MQ=255)
         cATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCtt                                                                                                                                             <  1:530857/150‑1 (MQ=255)
          aTCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTa                                                                                                                                            <  1:2078796/150‑1 (MQ=255)
           tCGGCCTCGACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTaa                                                                                                                                           >  2:2550325/1‑150 (MQ=255)
            cGGCCTCTACCGGCTCGGCTTTGGTCCACATTTTGACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCCGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTaag                                                                                                                                          <  1:1450476/150‑2 (MQ=255)
            cGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAAt                                                                                                                                          <  2:4257704/150‑1 (MQ=255)
               ccTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAAtttt                                                                                                                                       <  1:834942/150‑1 (MQ=255)
                 tctACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAAtttttt                                                                                                                                     >  1:1793167/1‑150 (MQ=255)
                   tACCGGCACGCCTTTGGTCCACATTATTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTaa                                                                                                                                   <  1:424408/150‑1 (MQ=255)
                     ccGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAAtt                                                                                                                                 <  2:4461093/150‑1 (MQ=255)
                      cGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTc                                                                                                                                >  2:2359922/1‑150 (MQ=255)
                      cGGCACGCCTTTGGTCCACATTTTGACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAATTTTTTAATTc                                                                                                                                <  2:1351811/150‑1 (MQ=255)
                       ggCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCa                                                                                                                               >  1:198744/1‑150 (MQ=255)
                           cGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGt                                                                                                                           <  2:3932393/150‑1 (MQ=255)
                            gCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGtt                                                                                                                          <  1:834314/150‑1 (MQ=255)
                              cTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAATTTTTTAATTCATTGtttt                                                                                                                        <  2:4619342/150‑1 (MQ=255)
                               tttGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTa                                                                                                                       <  1:2300354/150‑1 (MQ=255)
                                ttGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTaa                                                                                                                      <  2:2636271/150‑1 (MQ=255)
                                ttGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTaa                                                                                                                      >  2:4719032/1‑150 (MQ=255)
                                 tGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTaaa                                                                                                                     <  2:872854/150‑1 (MQ=255)
                                     ccACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAg                                                                                                                 >  2:248887/1‑150 (MQ=255)
                                       acaTTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGAt                                                                                                               >  2:3570672/1‑150 (MQ=255)
                                        caTTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGAtt                                                                                                              >  1:2842625/1‑150 (MQ=255)
                                         aTTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATta                                                                                                             <  1:2778462/150‑1 (MQ=255)
                                          tttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtat                                                                                                            <  1:246731/150‑1 (MQ=255)
                                          tttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtat                                                                                                            <  2:2997153/150‑1 (MQ=255)
                                           ttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtata                                                                                                           >  1:3412596/1‑150 (MQ=255)
                                           ttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtata                                                                                                           <  2:309219/150‑1 (MQ=255)
                                           ttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtata                                                                                                           >  1:1013395/1‑150 (MQ=255)
                                           ttttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATtata                                                                                                           >  1:1719346/1‑150 (MQ=255)
                                            tttACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTatat                                                                                                          <  2:3582752/150‑1 (MQ=255)
                                              tACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCt                                                                                                        >  1:325247/1‑150 (MQ=255)
                                               aCCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCtt                                                                                                       <  2:1691847/150‑1 (MQ=255)
                                                  ggggCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTAc                                                                                                    <  2:1332697/150‑1 (MQ=255)
                                                  ggggCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTAc                                                                                                    <  2:3663012/150‑1 (MQ=255)
                                                   gggCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAACATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACg                                                                                                   >  1:928433/1‑150 (MQ=255)
                                                       aTTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCg                                                                                               <  1:137558/150‑1 (MQ=255)
                                                                 gTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGgaga                                                                                     <  1:2135797/150‑1 (MQ=255)
                                                                       aattcgtaattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTccc                                                                               <  2:2816609/150‑1 (MQ=255)
                                                                            gtaattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCaaa                                                                          <  2:739949/150‑1 (MQ=255)
                                                                             taattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAg                                                                         <  2:1107473/150‑1 (MQ=255)
                                                                             taattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAg                                                                         <  2:1108539/150‑1 (MQ=255)
                                                                             taattcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAg                                                                         >  1:830413/1‑150 (MQ=255)
                                                                                ttcATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCa                                                                      >  2:2540509/1‑150 (MQ=255)
                                                                                  cATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAAt                                                                    >  2:2187181/1‑150 (MQ=255)
                                                                                   aTTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTCTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAAtt                                                                   <  2:4157794/150‑1 (MQ=255)
                                                                                     tttGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTa                                                                 <  2:4226878/150‑1 (MQ=255)
                                                                                        gTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCg                                                              >  1:554467/1‑150 (MQ=255)
                                                                                        gTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCg                                                              >  1:2943382/1‑150 (MQ=255)
                                                                                         tttttCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGa                                                             <  1:3975794/150‑1 (MQ=255)
                                                                                             tctcTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAAt                                                         <  2:4043970/150‑1 (MQ=255)
                                                                                             tctcTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAAt                                                         <  2:1456734/150‑1 (MQ=255)
                                                                                                 ttttCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGaaa                                                     >  1:412041/1‑150 (MQ=255)
                                                                                                  tttCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAat                                                    >  2:1208227/1‑150 (MQ=255)
                                                                                                  tttCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAat                                                    >  1:4647742/1‑150 (MQ=255)
                                                                                                    tCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAatat                                                  <  1:3820168/150‑1 (MQ=255)
                                                                                                       ttgttTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGt                                               <  2:3238921/150‑1 (MQ=255)
                                                                                                         gttTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATCAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCt                                             <  2:3884210/150‑1 (MQ=255)
                                                                                                           ttGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCttt                                           <  1:2012590/150‑1 (MQ=255)
                                                                                                           ttGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCttt                                           >  2:2314881/1‑150 (MQ=255)
                                                                                                            tGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCtttt                                          >  2:697050/1‑150 (MQ=255)
                                                                                                              cTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTtat                                        <  1:3044851/150‑1 (MQ=255)
                                                                                                               tGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTtata                                       <  1:4211150/150‑1 (MQ=255)
                                                                                                                 tCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAACAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAg                                     <  1:4367145/150‑1 (MQ=255)
                                                                                                                   cTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGa                                   >  1:1513600/1‑150 (MQ=255)
                                                                                                                   cTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGa                                   >  1:3718938/1‑150 (MQ=255)
                                                                                                                   cTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGa                                   >  1:1732711/1‑150 (MQ=255)
                                                                                                                   cTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGa                                   >  1:1308177/1‑150 (MQ=255)
                                                                                                                   cTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGa                                   >  1:4137705/1‑150 (MQ=255)
                                                                                                                     gATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGata                                 <  1:377175/150‑1 (MQ=255)
                                                                                                                      aTAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGatat                                >  2:1617112/1‑150 (MQ=255)
                                                                                                                       tAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATc                               <  1:484407/150‑1 (MQ=255)
                                                                                                                        aagtGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCt                              <  2:4670939/147‑1 (MQ=255)
                                                                                                                        aaTTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCt                              <  2:4671000/150‑1 (MQ=255)
                                                                                                                        aaTTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCt                              <  1:3570672/150‑1 (MQ=255)
                                                                                                                         aTTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTa                             <  2:2105061/150‑1 (MQ=255)
                                                                                                                         aTTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTa                             <  1:2119538/150‑1 (MQ=255)
                                                                                                                             cAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGAt                         >  1:1103426/1‑150 (MQ=255)
                                                                                                                                   tcgtGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTaa                   >  1:612381/1‑150 (MQ=255)
                                                                                                                                       gCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACa               <  1:2209473/150‑1 (MQ=255)
                                                                                                                                       gCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACa               <  1:791755/150‑1 (MQ=255)
                                                                                                                                         cAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGt             >  1:2605168/1‑150 (MQ=255)
                                                                                                                                          aGAAAATTAAACAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtg            <  2:374100/150‑1 (MQ=255)
                                                                                                                                            aaaaTTAAATAAGCAGCTTAATTTTTTAATACATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggc          <  1:3305322/150‑1 (MQ=255)
                                                                                                                                             aaaTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggct         >  2:1013518/1‑150 (MQ=255)
                                                                                                                                             aaaTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggct         >  1:1572207/1‑150 (MQ=255)
                                                                                                                                              aaTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggctt        >  2:4622585/1‑150 (MQ=255)
                                                                                                                                               aTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttt       <  2:3077317/150‑1 (MQ=255)
                                                                                                                                               aTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttt       <  2:2324173/150‑1 (MQ=255)
                                                                                                                                               aTTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttt       <  1:2260382/150‑1 (MQ=255)
                                                                                                                                                ttAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggctttt      >  1:3461452/1‑150 (MQ=255)
                                                                                                                                                 tAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttg     <  2:1989065/150‑1 (MQ=255)
                                                                                                                                                  aaaTAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttgg    >  2:2286791/1‑150 (MQ=255)
                                                                                                                                                  aaaTAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttgg    >  2:520717/1‑150 (MQ=255)
                                                                                                                                                  aaaTAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttgg    >  2:1000891/1‑150 (MQ=255)
                                                                                                                                                    aTAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttggct  >  2:1979658/1‑150 (MQ=255)
                                                                                                                                                    aTAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGtggcttttggct  >  1:2799311/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTGACGCGCATCGGCCTCTACCGGCACGCCTTTGGTCCACATTTTTACCGGGGCATTTTCAGTGGTCAGTAATTCGTAATTCATTTTGTTTTTCTCTTTTCGTTGTTTGCTGTCCTGATAATTGCAACCGTCGTGCCAGAAAATTAAATAAGCAGCTTAATTTTTTAATTCATTGTTTTAAAAAAGATTATATCTTTACGTCCGTAACCGGAGATTTCCCGCAAAGCCAATTTACCGATAATGAAATATCGTCTTTTATAAGGATATCTAAGATGCGTAAAACAGTGGCTTTTGGCT  >  CP009273/3551354‑3551651

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: