Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,7250AG58.9% 29.4 / 81.8 74intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base G (23/20);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTT  >  NC_002947/699581‑699805
                                                                                                                                                  |                                                                                
tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGaa                                                                             <  1:1150354/150‑1 (MQ=255)
   gCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAacac                                                                          <  2:313238/150‑1 (MQ=255)
   gCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAacac                                                                          >  2:721506/1‑150 (MQ=38)
    cGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACAcc                                                                         <  1:559958/150‑1 (MQ=255)
        tCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGtt                                                                     <  2:1128559/150‑1 (MQ=255)
           aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCa                                                                  >  2:389555/1‑150 (MQ=38)
            tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCAc                                                                 >  1:802363/1‑150 (MQ=39)
               ggCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCActgc                                                              >  2:311096/2‑150 (MQ=37)
               cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCActgc                                                              >  2:22726/1‑150 (MQ=38)
               cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCActgc                                                              >  2:993224/1‑150 (MQ=39)
               cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCActgc                                                              >  2:717384/1‑150 (MQ=39)
                gCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCActgct                                                             <  1:662885/150‑1 (MQ=255)
                gCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCActgct                                                             >  2:398030/1‑150 (MQ=38)
                 cTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCActgctg                                                            <  1:717384/150‑1 (MQ=38)
                  ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                           >  1:848169/1‑150 (MQ=38)
                   tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                          <  1:398030/150‑1 (MQ=38)
                    ggCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                         <  1:599384/150‑1 (MQ=38)
                       tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTTATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGAt                                                      >  1:1062883/1‑150 (MQ=255)
                         caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa                                                    >  1:647145/1‑150 (MQ=255)
                           ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                  >  1:211029/1‑150 (MQ=38)
                           ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGTTGGTGGATCAgg                                                  >  2:1166802/1‑150 (MQ=255)
                               tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGTTGGTGGATCAGGCTaa                                              <  1:1166802/150‑1 (MQ=255)
                               tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                              <  1:993224/150‑1 (MQ=39)
                               tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCATTCTTTAAAAACTCAGCTATTAGTGGTAAATAAGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                              <  1:243684/150‑1 (MQ=11)
                                ttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                             <  1:1274028/150‑1 (MQ=255)
                                ttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                             <  2:647145/150‑1 (MQ=255)
                                ttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                             <  2:1168185/150‑1 (MQ=39)
                                    ccGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTaa                                         >  2:1311202/1‑150 (MQ=38)
                                       tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                          >  1:570547/1‑130 (MQ=40)
                                       tACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                          <  2:570547/130‑1 (MQ=40)
                                         cGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                                    <  1:305082/150‑1 (MQ=255)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  <  1:389555/150‑1 (MQ=255)
                                           cAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  <  1:22726/150‑1 (MQ=255)
                                             gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGaa                                                                             <  1:975621/105‑1 (MQ=255)
                                             gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGaa                                                                             >  2:975621/1‑105 (MQ=40)
                                             gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  <  2:847890/148‑1 (MQ=38)
                                             gTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  >  1:847890/1‑148 (MQ=38)
                                                 cTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTCCCGATGGATCAGTCTAAGGTAAAACTTGTCAGTTc                            >  2:134448/1‑150 (MQ=11)
                                                  ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                           <  1:134448/150‑1 (MQ=34)
                                                    gTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGc                         >  2:489156/1‑150 (MQ=35)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                        >  1:1093572/1‑150 (MQ=25)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  >  1:1224394/1‑140 (MQ=37)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  <  2:1224394/140‑1 (MQ=37)
                                                     tCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                        >  1:570541/1‑150 (MQ=35)
                                                      cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGttt                                                                    <  2:641400/105‑1 (MQ=38)
                                                      cAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGttt                                                                    >  1:641400/1‑105 (MQ=38)
                                                       aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                      >  2:820098/1‑150 (MQ=25)
                                                       aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                      >  1:1360580/1‑150 (MQ=25)
                                                        gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGt                                                                      >  2:809973/1‑101 (MQ=255)
                                                        gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGt                                                                      <  1:809973/101‑1 (MQ=255)
                                                        gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                <  2:845783/139‑1 (MQ=37)
                                                        gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                >  1:845783/1‑139 (MQ=37)
                                                         aaCTTAGAAATGAACATTCAGCGCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGttt                                                                    <  1:657005/102‑1 (MQ=255)
                                                         aaCTTAGAAATGAACATTCAGCGCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGttt                                                                    >  2:657005/1‑102 (MQ=255)
                                                          aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTg                                                                               <  1:686061/90‑1 (MQ=255)
                                                          aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTg                                                                               >  2:686061/1‑90 (MQ=255)
                                                          aCTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                   <  2:570541/150‑1 (MQ=35)
                                                            ttAGAAATGAACATTCCGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                             <  1:851908/122‑1 (MQ=255)
                                                            ttAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                             >  2:851908/1‑122 (MQ=25)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTGAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCGCTGCTGGTGGATCAGGCTAAGGTAAAATTtg                                   >  1:1335039/1‑128 (MQ=255)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGTATCAGGCTAAGGTAAAATTtg                                   <  2:1335039/128‑1 (MQ=255)
                                                                aaaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa             >  2:879838/1‑150 (MQ=25)
                                                                 aaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtg                                                         >  1:1283479/1‑105 (MQ=255)
                                                                 aaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCtggtg                                                         <  2:1283479/105‑1 (MQ=255)
                                                                 aaTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc            <  1:1048506/150‑1 (MQ=35)
                                                                  aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                   <  1:1154804/142‑1 (MQ=25)
                                                                  aTGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                   >  2:1154804/1‑142 (MQ=25)
                                                                    gAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  <  1:791030/125‑1 (MQ=25)
                                                                    gAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                  >  2:791030/1‑125 (MQ=37)
                                                                        aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                              <  1:1093547/109‑1 (MQ=255)
                                                                        aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                              >  2:1093547/1‑109 (MQ=38)
                                                                        aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa     >  1:1230811/1‑150 (MQ=25)
                                                                        aTTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGaa     >  1:350122/1‑150 (MQ=25)
                                                                           cAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAttt  >  1:1208537/1‑150 (MQ=32)
                                                                                                                                                  |                                                                                
TCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATAT‑‑GTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTT  >  NC_002947/699581‑699805

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