Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 699,725 A→G 58.5% intergenic (+115/‑129) PP_t06 → / → PP_23SD tRNA‑Ala/23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,7250AG58.5% 33.1 / 98.2 89intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (22/15);  new base G (28/24);  total (50/39)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG  >  NC_002947/699579‑699808
                                                                                                                                                     |                                                                                   
ggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTg                                                                                  <  1:259082/150‑1 (MQ=38)
  tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGAc                                                                                    <  2:1178855/146‑1 (MQ=39)
  tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGAc                                                                                    >  1:1178855/1‑146 (MQ=39)
  tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGaa                                                                                <  2:1374131/150‑1 (MQ=39)
   cAGCGTTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAac                                                                               >  2:835327/1‑150 (MQ=255)
   cAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAac                                                                               <  2:1949830/150‑1 (MQ=255)
       ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCa                                                                           <  2:1669665/150‑1 (MQ=255)
        gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCt                                                                                   <  1:1108960/141‑1 (MQ=38)
        gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCt                                                                                   >  2:1108960/1‑141 (MQ=38)
        gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAg                                                                          >  2:1366977/1‑150 (MQ=39)
         ttagATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGt                                                                         >  2:1133832/4‑150 (MQ=255)
         ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGt                                                                         >  1:856084/1‑150 (MQ=38)
            gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGt                                                                         >  1:2384360/1‑147 (MQ=37)
            gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGt                                                                         <  2:2384360/147‑1 (MQ=38)
            gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTc                                                                      <  2:2405130/150‑1 (MQ=38)
             aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCa                                                                           <  2:316503/144‑1 (MQ=38)
             aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCa                                                                           >  1:316503/1‑144 (MQ=38)
                ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCActg                                                                  >  1:984347/1‑150 (MQ=38)
                  gCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCActgct                                                                >  1:1754792/1‑150 (MQ=39)
                    ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                              <  1:924386/150‑1 (MQ=255)
                    ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                              >  1:69351/1‑150 (MQ=255)
                     tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                             <  1:1366977/150‑1 (MQ=37)
                         tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt                                                         >  1:2092589/1‑150 (MQ=39)
                         tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt                                                         >  1:862922/1‑150 (MQ=39)
                         tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt                                                         <  1:359250/150‑1 (MQ=39)
                           caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:1754792/150‑1 (MQ=39)
                            accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGa                                                                                     >  1:2056665/1‑119 (MQ=255)
                            accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGa                                                                                     <  2:2056665/119‑1 (MQ=255)
                            accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAg                                                      >  1:2274645/1‑150 (MQ=255)
                            accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                      >  2:1381876/1‑150 (MQ=38)
                             ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                     >  2:480466/1‑150 (MQ=255)
                             ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                     <  2:856084/150‑1 (MQ=37)
                               acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCt                                                   >  1:2261971/1‑150 (MQ=255)
                                cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTa                                                  <  2:1222167/150‑1 (MQ=255)
                                    cGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGt                                              <  1:480466/150‑1 (MQ=255)
                                       cGTACGCAGTTACTTGTCAGAACTTAGGAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GTTAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaaa                                           >  1:2131792/1‑150 (MQ=255)
                                       cGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGCTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGTCTAAGGTaaa                                           >  2:982990/1‑150 (MQ=255)
                                        gTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCGTGCTAAGGTaaaa                                          >  2:441801/1‑150 (MQ=255)
                                         tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAATACCAGTTTGAATGCTGGTGGATCAGGCTAAGGTAAAAt                                         >  2:864581/1‑150 (MQ=255)
                                            gCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCa                                                                           >  1:1276079/1‑113 (MQ=255)
                                            gCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCa                                                                           <  2:1276079/113‑1 (MQ=255)
                                            gCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg                                      >  1:1890822/1‑150 (MQ=255)
                                             cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  1:732177/150‑1 (MQ=38)
                                             cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  1:2215152/150‑1 (MQ=255)
                                               gTTACTTGTCAGAACTTAGAAATGGACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTCTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:322784/148‑1 (MQ=255)
                                               gTTACTTGTCAGAACTTAGAAATGGACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTCTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     >  1:322784/1‑148 (MQ=255)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                   <  1:1079773/150‑1 (MQ=38)
                                               gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                   <  2:1797026/150‑1 (MQ=37)
                                                ttACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                  <  2:723367/150‑1 (MQ=25)
                                                   cTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTc                               >  2:1256920/1‑150 (MQ=25)
                                                   cTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTc                               >  2:1844567/1‑150 (MQ=25)
                                                      gTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc                            >  2:1782737/1‑150 (MQ=25)
                                                      gTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACAAGTTTCACTGATGGTGGATCAGGCTAAGATAAAATTTGTGAGTTCTGc                            >  2:875195/1‑150 (MQ=255)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAgg                                               <  1:2116993/130‑1 (MQ=25)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAgg                                               >  2:2116993/1‑130 (MQ=37)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           <  1:1844567/150‑1 (MQ=21)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           <  2:174641/150‑1 (MQ=25)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  2:1439652/1‑150 (MQ=35)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTT‑‑TCCCTCGAATTGGCCAGAGATATATAGACTGAACACCAGTTTCACTGCTGCTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  2:1289406/1‑150 (MQ=11)
                                                       tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCCTTAAAAATTCGGTTATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  2:885167/1‑150 (MQ=255)
                                                       tCAGAACTTAGAAATGAACATTCAACTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  1:2407796/1‑150 (MQ=21)
                                                         aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                              <  1:1766073/113‑1 (MQ=255)
                                                         aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                              >  2:1766073/1‑113 (MQ=255)
                                                         aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                         >  2:2123049/1‑150 (MQ=25)
                                                          gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGa                        >  2:1404051/1‑150 (MQ=25)
                                                           aaCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGa                                                                                     <  2:410710/88‑1 (MQ=38)
                                                           aaCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGa                                                                                     >  1:410710/1‑88 (MQ=39)
                                                            aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                      <  2:2274645/150‑1 (MQ=25)
                                                            aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                      <  1:1782737/150‑1 (MQ=25)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    <  1:1404051/150‑1 (MQ=25)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc                            >  1:1454229/1‑142 (MQ=35)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc                            <  2:1454229/142‑1 (MQ=35)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  2:825279/1‑150 (MQ=35)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  2:476586/1‑150 (MQ=35)
                                                              ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  2:2120608/1‑150 (MQ=35)
                                                               tAGAAATGAACATTCAGCTCGAATGTTGATTTCTGTCTGTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag                   <  1:384065/150‑1 (MQ=12)
                                                               tAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGATAAAATTTGTGAGTTCTGCTCGAAagag                   <  1:1172986/150‑1 (MQ=255)
                                                               caGAAATCAACATTCGAGCTGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag                   >  1:1051068/2‑150 (MQ=11)
                                                                aGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGc                  >  1:1776869/1‑150 (MQ=25)
                                                                  aaaTAAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                <  1:875195/150‑1 (MQ=21)
                                                                    aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg              >  2:2102178/1‑150 (MQ=25)
                                                                    aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg              >  1:1430134/1‑150 (MQ=25)
                                                                    aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGGTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTAAGTTTTGGTAGAAAGAGCAACg              >  1:288388/1‑150 (MQ=11)
                                                                          aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTTAGTTCTGCTCGAAAGAGCAACGTACGaa        >  1:1362690/1‑150 (MQ=25)
                                                                            tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGt             >  2:2291487/1‑143 (MQ=18)
                                                                            tCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGt             <  1:2291487/143‑1 (MQ=18)
                                                                               aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt    >  2:1825560/1‑150 (MQ=17)
                                                                                gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGTTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTc   <  1:441801/150‑1 (MQ=14)
                                                                                gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                                <  2:2012849/105‑1 (MQ=14)
                                                                                gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                                >  1:2012849/1‑105 (MQ=14)
                                                                                gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTc   >  1:485581/1‑150 (MQ=31)
                                                                                 cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGGCTGAACACCCGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg  <  2:1883539/150‑1 (MQ=11)
                                                                                 cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg  <  1:476586/150‑1 (MQ=11)
                                                                                                                                                     |                                                                                   
GGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGT‑‑GATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG  >  NC_002947/699579‑699808

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