Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947699,7250AG49.6% 41.3 / inf 267intergenic (+115/‑129)PP_t06/PP_23SDtRNA‑Ala/23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (62/72);  new base G (68/64);  total (130/137)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 269 total aligned reads displayed.

GGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG  >  NC_002947/699576‑699808
                                                                                                                                                      |                                                                                   
ggaggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGa                                                                                     <  1:2965414/150‑1 (MQ=37)
  aggTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACt                                                                                   <  1:5037588/150‑1 (MQ=37)
   ggTCAGCGGCTCGATCCCGCTTGGCTCCACCACTTACGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTTATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTg                                                                                  <  1:3440572/150‑1 (MQ=255)
     tCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGaa                                                                                <  2:1698435/150‑1 (MQ=37)
       aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAaca                                                                              <  2:1374827/150‑1 (MQ=37)
        gCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGCTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAacac                                                                             <  1:1750454/150‑1 (MQ=37)
           gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                          >  1:4179059/1‑150 (MQ=37)
           gTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAg                                                                          >  1:152058/1‑150 (MQ=37)
             tCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACGAGtt                                                                        >  2:3040785/1‑150 (MQ=255)
               gATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTc                                                                      >  1:1003921/1‑150 (MQ=37)
                aTCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCa                                                                     >  1:2222631/1‑150 (MQ=37)
                 tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc                                                                    <  2:3155703/150‑1 (MQ=37)
                   ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActg                                                                  >  2:2741949/1‑150 (MQ=37)
                   ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActg                                                                  >  2:88148/1‑150 (MQ=37)
                    cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgc                                                                 <  1:281668/150‑1 (MQ=37)
                       ttGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtgg                                                              >  2:3935811/1‑150 (MQ=37)
                        tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGGTTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATCTGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                             <  2:1311182/150‑1 (MQ=255)
                        tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt                                                             <  1:88148/150‑1 (MQ=37)
                        tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt                                                             <  1:2484818/150‑1 (MQ=37)
                        tGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt                                                             <  1:3376284/150‑1 (MQ=37)
                          gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                                           >  2:347092/1‑150 (MQ=37)
                          gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtgg                                                           <  2:4772305/150‑1 (MQ=37)
                            tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGAt                                                         >  2:1177951/1‑150 (MQ=37)
                            tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGAt                                                         >  1:2680033/1‑150 (MQ=37)
                            tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt                                                         <  1:3537318/150‑1 (MQ=37)
                            tCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt                                                         >  2:1535759/1‑150 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:2376487/150‑1 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:1145390/150‑1 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa                                                       >  1:4332145/1‑150 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:4751814/150‑1 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:3474168/150‑1 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  1:976134/150‑1 (MQ=37)
                              caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTATTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCa                                                       <  2:2693998/150‑1 (MQ=255)
                               accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                      <  2:4067482/150‑1 (MQ=37)
                               accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAg                                                      >  2:652204/1‑150 (MQ=37)
                               accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACAGTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGGTCTGTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAg                                                      >  1:4112717/1‑150 (MQ=255)
                                ccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg                                                     <  1:3594976/150‑1 (MQ=37)
                                 cacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGc                                                    >  2:2359353/1‑150 (MQ=37)
                                  acTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCt                                                   >  2:1012028/1‑150 (MQ=37)
                                    tttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTaa                                                 <  1:3779080/150‑1 (MQ=37)
                                     ttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                                >  1:1012708/1‑150 (MQ=37)
                                     ttCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAg                                                <  1:139499/150‑1 (MQ=37)
                                        gCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTa                                             >  2:2153573/1‑150 (MQ=37)
                                             aCGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAtt                                        >  2:598757/1‑150 (MQ=37)
                                              cGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATCTGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                                       <  2:3345630/150‑1 (MQ=37)
                                              cGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAAttt                                       <  1:164307/150‑1 (MQ=37)
                                               gCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtg                                      <  2:2811774/150‑1 (MQ=37)
                                                ctgTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:3689485/148‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:4259712/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:778040/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  1:4886327/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  1:320526/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:1065199/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:3402101/150‑1 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     >  2:4373197/1‑150 (MQ=37)
                                                cAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTCGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt                                     <  2:4512654/150‑1 (MQ=37)
                                                  gTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa                                   <  2:1644794/150‑1 (MQ=37)
                                                   ttACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCGTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACAGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAg                                  >  2:3871639/1‑150 (MQ=25)
                                                       ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                              >  2:4883416/1‑150 (MQ=21)
                                                       ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggtg                                                            <  1:2283260/120‑1 (MQ=35)
                                                       ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                              >  2:5043656/1‑150 (MQ=34)
                                                       ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                              >  1:1828281/1‑150 (MQ=34)
                                                       ttGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCt                              >  2:4202132/1‑150 (MQ=34)
                                                        tGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg                             <  1:796041/150‑1 (MQ=34)
                                                        tGTCAGAACTTAGAAATGAACAGTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTg                             <  2:2046506/150‑1 (MQ=21)
                                                          tCAGAACTTAGAAATGAACATTC‑GGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  2:4578598/1‑150 (MQ=21)
                                                          tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATATATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAATGTAAAATTTGTGAGTTCTGCt                           <  1:3885377/150‑1 (MQ=21)
                                                          tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  1:3626976/1‑150 (MQ=21)
                                                          tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  1:4491423/1‑150 (MQ=34)
                                                          tCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCt                           >  1:23239/1‑150 (MQ=34)
                                                            aGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCg                         >  1:4655716/1‑150 (MQ=21)
                                                               aCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTTAAATTTCGGTTATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGTTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaat                      <  1:837305/150‑2 (MQ=11)
                                                               aCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                      <  2:4692311/150‑1 (MQ=34)
                                                               aCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                      <  1:2112943/150‑1 (MQ=34)
                                                               aCTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa                      <  1:4781968/150‑1 (MQ=34)
                                                                 ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  1:4007243/1‑150 (MQ=21)
                                                                 ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  1:2747317/1‑150 (MQ=34)
                                                                 ttAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAaga                    >  2:400551/1‑150 (MQ=34)
                                                                  tAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag                   >  2:1769567/1‑150 (MQ=21)
                                                                   aGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGc                  >  2:1927909/1‑150 (MQ=34)
                                                                    gAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCa                 >  1:4401165/1‑150 (MQ=34)
                                                                     aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCac                >  2:979442/1‑149 (MQ=21)
                                                                     aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                >  2:4289355/1‑150 (MQ=21)
                                                                     aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                <  2:152058/150‑1 (MQ=21)
                                                                     aaaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                <  2:92975/150‑1 (MQ=34)
                                                                      aaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc               <  1:1290499/150‑1 (MQ=21)
                                                                      aaTGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAAc               >  1:1143981/1‑150 (MQ=21)
                                                                         gAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTa            >  2:2307693/1‑150 (MQ=21)
                                                                         caaCATTCGAGCT‑GAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTATCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGATTTCTGCTCGAAAGAGCAACGTa            >  1:1678000/2‑150 (MQ=11)
                                                                          aaCATTCGAGCT‑GAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc           >  2:3022322/1‑150 (MQ=2)
                                                                          aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                <  2:2588383/145‑1 (MQ=34)
                                                                          aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCaa                >  1:2588383/1‑145 (MQ=34)
                                                                          aaCATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTAc           >  1:2016521/1‑150 (MQ=34)
                                                                             aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGtt                                                                        <  1:2779173/86‑1 (MQ=255)
                                                                             aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGtt                                                                        <  1:2779043/86‑1 (MQ=255)
                                                                             aTTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGtt                                                                        >  2:2779173/1‑86 (MQ=255)
                                                                               tC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGACCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCActgct                                                                >  2:1021392/1‑92 (MQ=255)
                                                                                  aGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAtttt    >  2:3218235/1‑150 (MQ=17)
                                                                                    cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg  <  2:383716/150‑1 (MQ=11)
                                                                                    cTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCg  <  2:4416336/150‑1 (MQ=11)
                                                                                                                                                      |                                                                                   
GGAGGTCAGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTC‑AGCTCGAATGTTGATTTCTGGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTTTCG  >  NC_002947/699576‑699808

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