Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913827,5070TG44.6% 164.5 / 18.7 92D492A (GAC→GCC) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (6/45);  new base G (41/0);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.54e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATG  >  NC_000913/827360‑827653
                                                                                                                                                   |                                                                                                                                                  
ccGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTGCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCaa                                                                                                                                                 <  2:642087/151‑1 (MQ=255)
ccGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCaa                                                                                                                                                 <  1:1195290/151‑1 (MQ=255)
     ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                                                                             <  1:848806/150‑1 (MQ=255)
     ggCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCg                                                                                                                                            <  2:85587/151‑1 (MQ=255)
       cGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGa                                                                                                                                          <  2:1098507/151‑1 (MQ=255)
       cGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGa                                                                                                                                          <  1:826866/151‑1 (MQ=255)
       cGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCggg                                                                                                                                          >  1:486893/1‑150 (MQ=255)
       cGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCgg                                                                                                                                           >  1:884784/1‑150 (MQ=255)
         aTTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGa                                                                                                                                          >  1:813881/1‑149 (MQ=255)
            aaTTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCg                                                                                                                                     <  1:152682/151‑1 (MQ=255)
               ttCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACTCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGGCgggt                                                                                                                                  >  1:976196/1‑149 (MQ=255)
                tCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCgggtt                                                                                                                                 >  1:488530/1‑148 (MQ=255)
                 cccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcg                                                                                                                                <  2:429464/151‑1 (MQ=255)
                 cccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                                                                                 <  2:1187456/150‑1 (MQ=255)
                  ccGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCGGGAGTCCGGtcgt                                                                                                                               >  1:92066/1‑151 (MQ=255)
                     cgGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGCGGGGGGCAACGCCGGAAGGCGGGtcgtcg                                                                                                                             >  2:627039/1‑150 (MQ=255)
                        tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTTGTCgaga                                                                                                                          >  1:532237/1‑150 (MQ=255)
                        tACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCCAAGCCGGGAAGCGGGTCGTCggga                                                                                                                          >  2:272162/1‑150 (MQ=255)
                            ccAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAg                                                                                                                     <  1:315252/151‑1 (MQ=255)
                            ccAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGGGAAAc                                                                                                                       >  2:14686/1‑149 (MQ=255)
                              aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                                                                                   <  2:1138278/151‑1 (MQ=255)
                              aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                                                                                   <  2:1013358/151‑1 (MQ=255)
                              aGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGa                                                                                                                    >  1:200912/1‑150 (MQ=255)
                                gCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATg                                                                                                                 <  2:1087050/151‑1 (MQ=255)
                                gCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGGCGGGaaaaaaaa                                                                                                                   >  1:864547/1‑144 (MQ=255)
                                 ccGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGGAACAGAATGt                                                                                                                >  1:355062/1‑151 (MQ=255)
                                   gATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                                                                              >  1:446292/1‑151 (MQ=255)
                                    aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGCGGGGGGGCAAGCCGGCAGCCGGTTAGCTGACAAACAAAAtatgcg                                                                                                             >  2:58932/1‑146 (MQ=255)
                                    aTGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCg                                                                                                             <  2:422807/151‑1 (MQ=255)
                                     tGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCCGAAGACGGTTCGTCGAGAAAAAGAAAGTcccg                                                                                                            >  2:643434/1‑149 (MQ=255)
                                      gCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGt                                                                                                           <  2:552790/151‑1 (MQ=255)
                                      gCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGTAACGCCCGAAATCGGGTCGTCGAGAAAAAAGATGTCCggg                                                                                                           >  1:1336603/1‑150 (MQ=255)
                                       cGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                                                                                          <  1:336674/151‑1 (MQ=255)
                                       cGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                                                                                          >  1:213947/1‑151 (MQ=255)
                                        ggTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAGTCGGGTTGTCGAGAAACAGAATGTCCGGTTc                                                                                                         >  1:1122338/1‑151 (MQ=255)
                                        ggTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTc                                                                                                         >  1:648019/1‑151 (MQ=255)
                                           cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                                                                                       <  2:451876/150‑1 (MQ=255)
                                           cGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTCGTTGAGGAAAAAAAAGTCCgggttatt                                                                                                      >  2:1129836/1‑145 (MQ=255)
                                                tatTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATc                                                                                                  >  1:426236/1‑150 (MQ=255)
                                                 atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCa                                                                                                 <  1:610316/150‑1 (MQ=255)
                                                 atTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCa                                                                                                 <  1:1283499/150‑1 (MQ=255)
                                                   tCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCCACGCCCGAAGGCGGTTCGTCGAGGAAAAAAAAGTCCGGGTTATGCATCAg                                                                                                >  2:635518/1‑149 (MQ=255)
                                                    ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTAACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACCCCGGAAGGCGGTTCGTCGAGAAACAGAATGTCCGGGTCATTCATCAGCg                                                                                              >  1:132601/1‑150 (MQ=255)
                                                      gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                                                                           <  2:271140/151‑1 (MQ=255)
                                                         tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACa                                                                                          <  1:871634/149‑1 (MQ=255)
                                                          cATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgc                                                                                       <  1:1174802/151‑1 (MQ=255)
                                                           aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                                                                                      <  1:499541/151‑1 (MQ=255)
                                                            tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGAAGGTTGGTCTGGAAAAAGAATATCCGGGTGATTGCTTAGCGCACAGgcc                                                                                      >  2:1049437/1‑150 (MQ=255)
                                                             ccATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcca                                                                                     <  2:1113600/150‑1 (MQ=255)
                                                              cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                   <  2:1243605/151‑1 (MQ=255)
                                                              cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                   <  1:1130067/151‑1 (MQ=255)
                                                              cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagc                                                                                   <  2:872150/151‑1 (MQ=255)
                                                                      gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc                                                                            <  1:311314/150‑1 (MQ=255)
                                                                       ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                                                                          <  1:1111923/151‑1 (MQ=255)
                                                                          ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTg                                                                        <  2:1268340/150‑1 (MQ=255)
                                                                           gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGtt                                                                      <  1:645156/151‑1 (MQ=255)
                                                                             gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGACGGGTCGTCGAAAAAAAGAAAGTTCGGGTCCTGCATCAGCGAAAAGGCCAGCGCCAGCCGCCGtttt                                                                    >  1:1168935/1‑150 (MQ=255)
                                                                                cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa                                                                  >  1:364730/1‑150 (MQ=255)
                                                                                cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaaa                                                                 <  2:393249/151‑1 (MQ=255)
                                                                                cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaaa                                                                 <  1:1065780/151‑1 (MQ=255)
                                                                                cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGACGGTTCGTCGAGAAACAGAAAGGCCGGTTCATGCATCAGCGGAAAGGGCAGCGCCACCCGCTGTTTaaa                                                                  >  1:621058/1‑150 (MQ=255)
                                                                                cATCACCGTGACGCCTTTCTCTACCATGCTGTGGATGTGCAGCCAAAATTCACGGAGGGTGAGGGGGTCAAGGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaaa                                                                 <  2:92063/151‑1 (MQ=255)
                                                                                 aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAAc                                                                <  1:354698/151‑1 (MQ=255)
                                                                                 aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAAc                                                                >  1:267177/1‑151 (MQ=255)
                                                                                    aCCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAAACt                                                              >  1:3148/1‑150 (MQ=255)
                                                                                     ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa                                                            <  1:352892/151‑1 (MQ=255)
                                                                                        tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTAAGGTCGTCGAGAAAAAAAAAGTCCCGGTCCCGCCTCAGCGAACAAGCCAGCGCCAGACGGCGTTTAAAAACTAAttg                                                         >  1:196404/1‑149 (MQ=255)
                                                                                        tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCCGCGCCAGCCGCTGTTTAAAAACTAATgg                                                         >  1:754583/1‑151 (MQ=255)
                                                                                         gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGc                                                        <  1:871961/151‑1 (MQ=255)
                                                                                         gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAACCTAATggg                                                        >  1:1099586/1‑150 (MQ=255)
                                                                                            gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGTTCATGCATCAGCGAACTGGGCAGCGCCAGCCGCCGTTTTAAAACTAATGGCAGt                                                     >  2:61218/1‑151 (MQ=255)
                                                                                            gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCCGAAGTCGGTTTGTCGAGAAACAGAATGTCCGGGTTATTCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAAACTAATGGCAgg                                                     >  1:786963/1‑150 (MQ=255)
                                                                                              cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAACCTAATGGCAGttt                                                   >  1:50315/1‑150 (MQ=255)
                                                                                                ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAAAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTTAAACCTAATGGCAGTTCAt                                                 >  1:764854/1‑151 (MQ=255)
                                                                                                   tctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTTAAACCTAATGGGAGTTTATCgg                                              >  1:1320448/1‑151 (MQ=255)
                                                                                                      aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCgg                                              <  1:852849/148‑1 (MQ=255)
                                                                                                             tGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGGAGGCGGTTCGGCGAGAAACAAAAAGTCCGGGTCATGCATCAGCGAACAAGGCAGCGCCAGCCGCCGTTTTAAAACTAATGGCAGGTCATCGGGGGGGTgggg                                    >  1:947518/1‑150 (MQ=255)
                                                                                                              gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAATGTTCGGGTCATGCATCAGCGAACAAGCCAGCGCCCGCCGCCGTTTAAAAACTAATGGCCGGTCATCGGTGGCGgggggg                                   >  2:568690/1‑149 (MQ=255)
                                                                                                                 gATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGc                                 <  2:886142/150‑1 (MQ=255)
                                                                                                                 gATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGc                                 <  1:88105/150‑1 (MQ=255)
                                                                                                                  aTGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGACGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAAACTAATGGCAGGTCATCggggggggggggggggg                               >  2:923429/1‑148 (MQ=255)
                                                                                                                   tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACCCCCGGAGGCGGGTTGTCGAGAAAAAGGATGGTCGGGTTATGCATCCGCGAAAAGGGCAGCGCCCGCCCCTGGTTAAAACATAATTGCCGGTCCTCCGGGGCGCgggggggggt                              >  2:1008292/1‑147 (MQ=255)
                                                                                                                          ccAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCCGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGGGGCGTGGGAGGGGATACttttt                       >  1:1350302/1‑150 (MQ=255)
                                                                                                                            aaaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAg                     <  2:813812/151‑1 (MQ=255)
                                                                                                                             aaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAgg                    <  2:792841/151‑1 (MQ=255)
                                                                                                                             aaaTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAgg                    <  2:320440/151‑1 (MQ=255)
                                                                                                                                        gggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTc         <  2:741750/151‑1 (MQ=255)
                                                                                                                                          gTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGc       >  2:678546/1‑151 (MQ=255)
                                                                                                                                             aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTc     >  2:865077/1‑150 (MQ=255)
                                                                                                                                             aGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTc     >  2:522098/1‑150 (MQ=255)
                                                                                                                                              gggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCAt   >  1:1318782/1‑151 (MQ=255)
                                                                                                                                               ggggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATg  <  1:938382/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGCTGGCGATTAATTTCCCGCGGTACACCAGGCCGATGCGGTCGCAATATTCCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAATGGCAGTTCATCGGTGGCGTGGGAGGCGATACTTTTCAGGCCGAACGCCTCGCTCATG  >  NC_000913/827360‑827653

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