Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3440CT46.0% 56.4 / 38.5 63S62F (TCT→TTT) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (3/31);  new base T (29/0);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 7.14e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604195‑4604492
                                                                                                                                                     |                                                                                                                                                    
aCTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                                                                                     <  2:744891/151‑1 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCCTGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                  <  2:713928/151‑1 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                                                                                  >  2:282278/1‑151 (MQ=255)
   cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttt                                                                                                                                                  <  2:1288442/151‑1 (MQ=255)
    aTAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                                                                                 <  1:1066748/151‑1 (MQ=255)
      aGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTttttctttttct                                                                                                                                               <  2:1089775/151‑1 (MQ=255)
        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttt                                                                                                                                               >  2:1037426/1‑147 (MQ=255)
        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGc                                                                                                                                             <  2:795193/151‑1 (MQ=255)
            aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                                                                                         <  2:1559011/151‑1 (MQ=255)
               cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGcc                                                                                                                                            >  1:1169687/1‑145 (MQ=255)
               cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                      <  1:741681/151‑1 (MQ=255)
                gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATga                                                                                                                                        >  2:1042466/1‑148 (MQ=255)
                gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                                                                                      <  1:1030560/150‑1 (MQ=255)
                     gAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa                                                                                                                                 <  2:810693/150‑1 (MQ=255)
                      agtagtATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAGCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                                                                               <  2:1463521/151‑1 (MQ=255)
                       gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                                                                              <  1:860863/151‑1 (MQ=255)
                       gtagtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAATTTTTTTTTTTTTTCCCTGTGAAGTGGgcg                                                                                                                              >  2:128894/1‑151 (MQ=255)
                         agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATTAGAGGTGa                                                                                                                                 >  2:450181/1‑146 (MQ=255)
                          gtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAggg                                                                                                                              >  1:424197/1‑146 (MQ=255)
                            aTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTAATTTTTTTTTTTTCTACCATGGGAAGTGAggg                                                                                                                              >  2:224391/1‑144 (MQ=255)
                               ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCagggggg                                                                                                                      >  1:80050/1‑146 (MQ=255)
                                 tATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGGCGAGAggg                                                                                                                     >  1:1097437/1‑150 (MQ=255)
                                  aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGGAGTGAGGGCaga                                                                                                                          >  1:1527404/1‑144 (MQ=255)
                                  aTTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCCTGAGAAGAGGGCGCagggaggg                                                                                                                     >  1:1345173/1‑149 (MQ=255)
                                    tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAggg                                                                                                                     >  1:564457/1‑147 (MQ=255)
                                     gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGAAATTTTTTTTTTTTCTGTCGTTTGAAGTGGGCGCAGAGAGGGCtt                                                                                                                  >  1:410365/1‑149 (MQ=255)
                                       gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATGAGAAGTGAGCGCAGAGAggggtgat                                                                                                                >  2:119034/1‑144 (MQ=255)
                                         gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTg                                                                                                            <  2:1180651/151‑1 (MQ=255)
                                          tttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTACTTTTTTTTTTTTCTTCCCTGGGAAGAGGGAGGAGAGggggggtatttt                                                                                                             >  1:1200874/1‑141 (MQ=255)
                                           ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                           <  1:827402/150‑1 (MQ=255)
                                                cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                                                                                      <  1:1046011/150‑1 (MQ=255)
                                                 aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgactg                                                                                                    <  1:953093/151‑1 (MQ=255)
                                                    aaGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGGCGAGAGGGATTATTTTTCctgactgta                                                                                                  >  2:1070470/1‑148 (MQ=255)
                                                         tGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTATGGCCTGGGTAGTGAGCGGAGAGCGGGGTTATCTTGTctgccttgtctagcg                                                                                            >  2:851860/1‑142 (MQ=255)
                                                           cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAt                                                                                          <  1:128894/151‑1 (MQ=255)
                                                             ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                                                                                        <  1:554679/151‑1 (MQ=255)
                                                              gTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTaa                                                                                       <  1:191541/151‑1 (MQ=255)
                                                               tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGAtt                                                                                         >  1:987707/1‑148 (MQ=255)
                                                                aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAg                                                                                      <  1:527223/150‑1 (MQ=255)
                                                                 aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGGGAGCGCAGAGAGGGATTTTCTTGCCTGACTGAACTGGCGATTAAGtt                                                                                    >  2:134639/1‑151 (MQ=255)
                                                                 aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                                                                                    <  1:144233/151‑1 (MQ=255)
                                                                    cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTCTGCCTTTCGGAGGGGGGGCCGAGGGGGAATATTTTTGCTGACTGAACCggggcattaattt                                                                                   >  1:640747/1‑138 (MQ=255)
                                                                     ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGGGCAGAGAGGGATTATTTTGCCCGACCGAAATGGGGGTTAAGttttc                                                                                 >  1:1558681/1‑148 (MQ=255)
                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                 <  2:1528370/149‑1 (MQ=255)
                                                                      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTcc                                                                                 <  1:1151648/149‑1 (MQ=255)
                                                                       ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                                                                               <  2:823708/150‑1 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCATTTTTTTTTTTTTTCCATTGGGAGGGGGCGCCTGGAGGGGATTTCTTGTCTGGCTGGACTGGCGGGTAAGATTCcccg                                                                              >  2:866267/1‑149 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGGGAGGGGTTATCTTTCCTGACTGGACTGGGGGTTAAGGTTTcgcg                                                                              >  1:1299267/1‑149 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGGGAGCGCAGAGAGGGGTTATTTTGGCTGACTGAAATGGGGATTAAGTTTTCGCGt                                                                             >  1:470458/1‑150 (MQ=255)
                                                                         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTTTTTCATTAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTGCCTGAACTTGCGATTTAGTTTTcgc                                                                               >  1:1198634/1‑148 (MQ=255)
                                                                                atcccTATACAAGTGAGGAGTCATTTCAGGATGCCATGTCGCGGATCGCGGTTGCGGAGGTCATTTTTTCTATTTCTGCTATGAGAAGGGAGCGCAGATAGGGATTATCTTGCCTGCCTGAACGGGCGATTAAGTTTCCGCGTACCCggg                                                                      <  2:1447207/147‑2 (MQ=255)
                                                                                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg                                                                     <  2:913496/151‑1 (MQ=255)
                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa                                                                 <  2:535085/149‑1 (MQ=255)
                                                                                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAg                                                                <  1:1135179/150‑1 (MQ=255)
                                                                                       taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGTATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCGTTTAg                                                                <  2:1006016/149‑1 (MQ=255)
                                                                                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCCGAGAGGGGTTATCTTTCCCGATTGGACTGGCGATTTAGTTTTCGCGGACCCGGCGGTTTGTTAtt                                                           >  2:1153503/1‑151 (MQ=255)
                                                                                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTAtt                                                           <  2:1308131/151‑1 (MQ=255)
                                                                                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACCGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTa                                                             <  2:69964/149‑1 (MQ=255)
                                                                                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACTGGCGATTAAGTTTCCCCGTACCCCGCGGTTAGTTAtt                                                           >  1:616452/1‑150 (MQ=255)
                                                                                                ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCTTTTTTTTTTTTCTTTCCTGGGAAGTGAGGGGAGAGGGGGGTTATCTTTCCTGGCTTAACTTGCGATTAAGATTTCGCGTACCCGGCGGTTTGTTATTgtg                                                        >  2:848376/1‑146 (MQ=255)
                                                                                                  aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga                                                   <  1:1417975/151‑1 (MQ=255)
                                                                                                    tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGGCATGAGAAATGGGCGCACAGGGGGATTTTCTTGCCTGACTGAACTTGCGATTAAGTTTCCGCGTAACCGGGGTTTTATTATTGTGgatgat                                                  >  1:447257/1‑150 (MQ=255)
                                                                                                        ttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGGGCGCCGAGAGGGATTATTTTTCCCGGCTGAACTGGCGGTTAAGTTTCCGCGGACCCCGCGTTTAGTTTTTTCGGATGATGatat                                             >  1:1297069/1‑151 (MQ=255)
                                                                                                           aGGATGCCATGTCTCGGATCGCGTTTGCGGGGGTCATTTTTTCTTTTTCTGTCATGATAATTGAGCGCAGAGAGGGATTATCTTGCCTGACTCAACTTGCGATTAAGTTTCCGCGTACCGGGCGTGTAGTTATTGCGGATGATGatat                                             <  2:682862/148‑1 (MQ=255)
                                                                                                                 ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGGGGTTAAGTTTTCGCGTACCCCGCGGTTAGTTTTTTCGGATGATGATATTTAAGCTc                                     >  1:1487889/1‑150 (MQ=255)
                                                                                                                 ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCg                                    <  2:532362/151‑1 (MQ=255)
                                                                                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCCGAATGAACGAGCGATTAAGTTTCCCCGTCACCGGCGTTTTTTTATTTCGGATGGTGATATTTAAGc                                       >  1:1176767/1‑147 (MQ=255)
                                                                                                                                     gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg                 <  2:46203/150‑1 (MQ=255)
                                                                                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCAt             >  2:1019727/1‑151 (MQ=255)
                                                                                                                                               ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGTTTCATTGTCGCCATCACCGc       <  2:13804/150‑1 (MQ=255)
                                                                                                                                                  tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAGGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg    >  2:1303679/1‑150 (MQ=255)
                                                                                                                                                    tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc  <  1:615293/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC  >  NC_000913/4604195‑4604492

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