Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 699,725 A→G 47.6% 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,7250AG47.6% 5.5 / 94.2 75intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (17/18);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTAT  >  NC_002947/699581‑699830
                                                                                                                                                 |                                                                                                         
tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa                                                                                                      <  2:1100612/150‑1 (MQ=39)
 cAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAac                                                                                                     >  2:57119/1‑150 (MQ=38)
    cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAACGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACAcc                                                                                                  >  1:1602763/1‑150 (MQ=255)
     ggTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCa                                                                                                 <  1:1707378/150‑1 (MQ=39)
      gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                                                >  2:764871/1‑150 (MQ=255)
       ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt                                                                                               <  1:394840/150‑1 (MQ=255)
       ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt                                                                                               >  2:537414/1‑150 (MQ=255)
       ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGt                                                                                               <  1:1147251/150‑1 (MQ=255)
       ttCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAAAACCAGt                                                                                               >  2:846030/1‑150 (MQ=255)
         cGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGttt                                                                                             <  1:537414/150‑1 (MQ=255)
          gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                                                                                            >  1:131909/1‑150 (MQ=38)
          gATACCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTAATTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTACTTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTAGGATATGTGATAGCTATACACTCACCACCAGTTTc                                                                                            >  2:66218/1‑150 (MQ=11)
           aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                                                                                           <  1:1110269/150‑1 (MQ=38)
           aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa                                                                                           >  1:1569806/1‑150 (MQ=39)
            tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATTTTATAGATATAGGCTGAACACCAGTTTCAc                                                                                          <  1:876211/150‑1 (MQ=255)
            tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc                                                                                          <  1:1720457/150‑1 (MQ=39)
             cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCTTATATGGGTTAGAGATAGACTGCAAAACCGTTTctct                                                                                         >  2:804740/1‑150 (MQ=14)
             cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct                                                                                         >  1:360480/1‑150 (MQ=39)
                 cTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActgctg                                                                                     >  2:1741735/1‑150 (MQ=39)
                  ttGGCTCCAGCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtgg                                                                                    >  2:963394/1‑150 (MQ=38)
                  ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtgg                                                                                    <  1:804740/150‑1 (MQ=39)
                   tgtcTCCACCACTTTCGCCTTACGCAATTATTTGTCAAAACTTAGAAATGAACATTCAGCTAGAATGTTTATTTCTGACGTTTGTCAGATCATCTTTTAAAAATTCGGATATGCAACAAATATACTCTTAACAACAGTTGCACTGCttgt                                                                                   <  2:300729/147‑1 (MQ=2)
                   tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt                                                                                   <  1:782397/150‑1 (MQ=255)
                   tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                                                   <  2:1166579/150‑1 (MQ=38)
                     gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                                                                 >  2:516954/1‑150 (MQ=38)
                       tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGAt                                                                               >  2:62714/1‑150 (MQ=255)
                        ccaccaCTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGacc                                                                              >  1:1052703/1‑148 (MQ=39)
                          accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAg                                                                            >  2:16666/1‑150 (MQ=38)
                          accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                                            <  2:360480/150‑1 (MQ=39)
                             acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtg                                                                                  <  2:884147/141‑1 (MQ=255)
                             acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtg                                                                                  >  1:884147/1‑141 (MQ=255)
                             acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCt                                                                         >  1:1603360/1‑150 (MQ=38)
                                     cGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaaa                                                                 <  1:66218/150‑1 (MQ=39)
                                      gTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaaaa                                                                <  1:1741735/150‑1 (MQ=39)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  1:1332656/150‑1 (MQ=255)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  2:1710071/150‑1 (MQ=37)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  2:1603360/150‑1 (MQ=255)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  2:1569806/150‑1 (MQ=39)
                                           cAGTTACTTGTACGAACTTAGAAATGAACATTCAGCTCTAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  2:903492/150‑1 (MQ=255)
                                                  ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                                                    >  1:1355344/1‑150 (MQ=25)
                                                  ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                                                    >  1:1028048/1‑150 (MQ=35)
                                                  ttGTCAGAACTCAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                                                    >  2:143208/1‑150 (MQ=25)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           >  2:1827083/1‑140 (MQ=25)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           >  2:1826935/1‑140 (MQ=25)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  1:1827083/140‑1 (MQ=37)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                                           <  1:1826935/140‑1 (MQ=37)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                                                 >  1:727251/1‑150 (MQ=35)
                                                      cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc                                                <  1:16666/150‑1 (MQ=25)
                                                      cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc                                                <  1:1311233/150‑1 (MQ=35)
                                                       aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                                               >  1:1756740/1‑150 (MQ=35)
                                                          aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                                            <  2:1355344/150‑1 (MQ=25)
                                                           gccagaaATCAACATTCGAGCTGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAG                                           <  2:1162842‑M2/147‑1 (MQ=255)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGA                                          >  2:773221‑M2/1‑150 (MQ=255)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct                                                                                         >  1:176526/1‑103 (MQ=38)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct                                                                                         <  2:176526/103‑1 (MQ=38)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                                           <  2:1513319/117‑1 (MQ=37)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                                           >  1:1513319/1‑117 (MQ=38)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGACTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGA                                          >  1:433120‑M2/1‑150 (MQ=255)
                                                            ttAGAAATGAAAATACAGCTCGAATTTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGA                                          >  1:1197500‑M2/1‑150 (MQ=255)
                                                               gAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCA                                       >  2:131909‑M2/1‑150 (MQ=255)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                                                       <  1:438861/117‑1 (MQ=255)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                                                       >  2:438861/1‑117 (MQ=255)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                                      <  2:1368483/150‑1 (MQ=35)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                                      >  2:628292/1‑150 (MQ=35)
                                                                   tGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc                                     <  2:30019/148‑1 (MQ=21)
                                                                   tGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc                                     >  1:30019/1‑148 (MQ=25)
                                                                    gAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTa                                  <  1:393599/150‑1 (MQ=25)
                                                                     aaCATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTCTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc                                 >  2:1591540/1‑150 (MQ=255)
                                                                       cATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa                                                                                                      >  1:110558/1‑79 (MQ=39)
                                                                       cATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa                                                                                                      <  2:110558/79‑1 (MQ=40)
                                                                             aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt                          >  1:364751/1‑150 (MQ=32)
                                                                              gCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTc                         <  1:963394/150‑1 (MQ=31)
                                                                                              cTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCTTTCGAGCAgaactcacaaattttaccttagcctgatcc         <  2:567715‑M2/150‑31 (MQ=255)
                                                                                              cTGACTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCTTTCGAGCAgaactcacaaattttaccttagcctgatcc         >  1:1102032‑M2/1‑120 (MQ=255)
                                                                                                     ttGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCTTTCGAGCAgaactcacaaattttaccttagcctgatccaccagca  >  1:567715‑M2/1‑113 (MQ=255)
                                                                                                                                                 |                                                                                                         
TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCGGCGAATGTCGTCTTCACAGTAT  >  NC_002947/699581‑699830

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