Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,686,0280TC54.9% 75.6 / 127.2 93E1049G (GAG→GGG) bcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base T (28/13);  new base C (0/50);  total (28/65)
Fisher's exact test for biased strand distribution p-value = 7.88e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-14

CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCC  >  NC_000913/3685878‑3686170
                                                                                                                                                       |                                                                                                                                                  
cgcgCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCt                                                                                                                                                    <  2:553494/151‑1 (MQ=255)
 gcgcCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCt                                                                                                                                                    >  1:1434732/1‑150 (MQ=255)
    ccGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCa                                                                                                                                                <  2:204629/151‑1 (MQ=255)
     cGTGTAGCCGAATCCCTGGCTGCTGCCCCCGTCCTTGGATTGCCGCGCAGCTTCTTCCTGACAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGTCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTcccg                                                                                                                                                <  1:1555785/150‑2 (MQ=255)
      gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                              >  2:918896/1‑151 (MQ=255)
      gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                              >  2:555296/1‑151 (MQ=255)
      gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                              >  2:362482/1‑151 (MQ=255)
      gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACAACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                              >  1:88699/1‑151 (MQ=255)
       tgtAGCCAAAGCCCTGGCTGCTGCCGCCAACGTTGGTTTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGGCCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                              <  2:913325/150‑1 (MQ=255)
           gCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                         >  1:172478/1‑151 (MQ=255)
            ccGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                         <  2:1355204/150‑1 (MQ=255)
            ccGAAGCCCTGACTGCTGCCCCCATCGTGGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCTACCACGAGCCAGACGCACCCAGCTCCCACGACCAg                                                                                                                                         <  1:379758/150‑1 (MQ=255)
               aaGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGttt                                                                                                                                     <  2:796690/150‑1 (MQ=255)
                 gCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTc                                                                                                                                  >  1:77632/1‑150 (MQ=255)
                 gCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTc                                                                                                                                  >  1:1196476/1‑150 (MQ=255)
                  cccTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTc                                                                                                                                  <  1:406895/149‑1 (MQ=255)
                         tgctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTc                                                                                                                       <  1:280238/151‑1 (MQ=255)
                                    tCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCa                                                                                                            >  2:1365774/1‑151 (MQ=255)
                                     ctttGGTTTGCCGCGCAGTTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGTTTTTGGTGCGTGAATGCGCCCCCGAGCCAGACGCCCCCAGCCCCCACGCCCATTTTTCCGTGCGCCCCCGCCACATCAc                                                                                                           <  1:955502/149‑1 (MQ=255)
                                     cGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCa                                                                                                            >  2:1061377/1‑150 (MQ=255)
                                        tGGTTTGCCGCCCGCCTTTTCCCTGCCGGTCGGTCGGGATAAAATCCATCAGGGGAAAACGGGGCATGTTTTTGGTGCGTGAATGCGCCCACGAGCCAGACGCCCCCAGCCCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCgg                                                                                                        <  2:1280475/151‑1 (MQ=255)
                                        tGGATGGCCGCGCAGCTTCTTCCTGCCAGTCGGTGGGGATCAGATTCATCAGCGGAAAACGCGGCATGTTTTTGGGGCGTGAATGCGCCCACGAGCCAGCCGCCCCCAGCCCCCCCGACCAGTTTTCCGTGCGCTCCCGCCACATCACCgg                                                                                                        <  1:1380794/151‑1 (MQ=255)
                                           tttcccgcccggtttctccCTGCCAGTCGGTGGGGATAAGATTCATCAGCGGATAGCGCGGCATTGTTTTGGGGCGTAAATGCGACCACGAGCCAGACGCCCCCAGCCCCCCCGACCGGTTTTCCGTGCGCTCCCGCCACATCACCGGTa                                                                                                      <  2:88728/137‑1 (MQ=255)
                                           ttcgCCGCGCATTTTCTCCCGGCCGGTCGGCGGGGATCAGATTCATCAGCGAATAACGCGGCATGTTTTTGGTGCGTGATTGCGCCCACGACCCAGCCCCCCCCAGCCCCCCCGCCCAGTTTTCGGTGCGCTCCCGCCACATCACCGGTAt                                                                                                     <  1:969266/148‑1 (MQ=255)
                                            ttGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATg                                                                                                    >  1:68026/1‑151 (MQ=255)
                                             tgcgccgcattttctccCCGCCAGTCGGTAGGGATCAGATTCATCAGCGGATACCGCGGCATGGTTTTGGGGCGGGAAGGCGCCCACGACCCAGCCGCCCCCTCCCCCCACGACCATTTTTCCGGGCGCCCCCGCCACATCACCGGTATgg                                                                                                   <  1:1674530/140‑1 (MQ=255)
                                             tGCCGCCCAGCTTTTCCCTGCGGGTCGGGCGGGATCAAATTCAACAGCGGAAAGCGCGCCGTGGTTTTGGGGCGGGAATGCGACCACGGGCCAGCCCCCCCCAGCTCCCGCGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATgg                                                                                                   <  1:40166/151‑1 (MQ=255)
                                                  cgcattttctccCTCCCAGTCGGTCGGAATCAGATTCATCAGCGAATACCGCGGCATGGTTTTGGTGGGTGAATGCGCCCCCGACCCAGACCCCCCCAGCCCCCCCACCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAAc                                                                                              <  1:1548749/145‑1 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                          >  1:665905/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                          >  1:1285119/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                          >  1:1363349/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                          >  1:322222/1‑151 (MQ=255)
                                                       tttcttcCTGCCGGTCGGTGGTGACGAGATCCGGGGGGGGATAGCGAGGCATGTTTTTGGTGGGTAAAGGCAACCACGAGCCAGACGCACCCGGCACCCACAACAAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGAc                                                                                           <  1:813147/148‑1 (MQ=255)
                                                          ttcccTGCCGGTGGGTCGGGATCAAATTCACCGGGGGAAACCGCGGCAGGTTTTGGGGGCGTGAATGCGCCCACGCCCCGGCCGCCCCCGGCCCCCCCGACCTTTTTTCGGTCCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTa                                                                                      <  1:348902/148‑1 (MQ=255)
                                                           tccccccgaGTGGGTGGGGATAAGATTCACCAGGGGAAAACGCGGCATTTTTTTGGGGCGAAAATGCGACCACAAGCCAGGCGCACCCAGCTCCCCCACCCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTAt                                                                                     <  2:442512/143‑1 (MQ=255)
                                                              cggccggtggggggggATAAGATCCACCGGGGAATACCGCGGCGGGGTTTGGGGGGGGGAAGGAGCCCACGGGCAAAACGCCCCCAGCCCCCCCGCCAGTTTTTCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                                  <  2:1129137/138‑1 (MQ=255)
                                                              cgcgcAGTGGGTCGGAATCAGATTCATCAGCGAATACCGCGGCTTGTTTTTGGGGCGTGAACGCGCCCACGAGCCAGCCGCACCCAGCCCCCACAACCTTTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                                  <  1:1568246/147‑1 (MQ=255)
                                                              cTGCCAGTGGGTCGGGATCAGATTCACCAGCGAAAACCGCGGCAGGTTTTTGGGGCGTAAAGGCGCCCCCGAGCCAGACGCCCCCGGCCCCCACGCCCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTcc                                                                                  <  2:1061395/151‑1 (MQ=255)
                                                                 ccggtgggtgggGAGCAGATTCATCAGCGGAAACCGCGGCTTGTTTTGGGGGCGTGAATGCGCCCACGAGCCGAACCCCCCCAGCCCCCACGCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTg                                                                                <  1:1496024/144‑1 (MQ=255)
                                                                     tcgggggggATCAGATCCATCAGCGGATACCGCGGCATGTTTTGGGTGCGTGAATGCGCCCACGAGCCAGACGCACCCAGCCCCCACGACCAGTTTTCCGGGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGa                                                                            <  1:402986/144‑1 (MQ=255)
                                                                       gggcggGACAAGATTCAGCGCCGGAAAACGCGGCATGTTTTGGGGGGGTGAATGCGCCCCCGAGCCAGACGCCCCCACCCCCCACGCCCATTTTTCCGGGCGCGCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACt                                                                          <  1:845578/147‑1 (MQ=255)
                                                                        gtcggGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGa                                                                            >  1:11348/1‑147 (MQ=255)
                                                                           gggacaaGATACACCGGCGAATACCGCGGCGGGGTTTGGGTGCGTGAATGCGACCCCGAGCCAGACGCACCCGGCCCCCACAACCGGTTTTCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtag                                                                      <  1:1356456/144‑1 (MQ=255)
                                                                           gggATAAGATTCATCAGGGGAAACCGGGGCTTGTTTTTGGTGGGGGAAGGCGCCCACGAGCCAGACGCCCCCAGCCCCCCCGACCAGTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtag                                                                      <  2:114913/150‑1 (MQ=255)
                                                                            ggATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtagta                                                                    >  1:731225/1‑151 (MQ=255)
                                                                            ggACCAAATTCCCCAGGGGAAAACGGGGCATGGTTTTGGGGCGGGAATGCGACCACGAGCCAGACCCACCCACCTCCCCCGACCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTgtagt                                                                     <  2:1706749/150‑1 (MQ=255)
                                                                              acaaaactcataaGCGGAAACGGCGTCATGTTTTGGGGGGTTGAGTCCGACCACGGCCAAGCCCCCCCCAGCCCCCACACCCATTTTCCCGTCCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                                  <  1:108052/144‑1 (MQ=255)
                                                                              aacAGATTCACCGCGGGAAAAGGCGGCATGGTTTGGGGGGTGAAATGCGACCGCGCCCCAGACGCCCCCAGCCCCCACGACCGGTTTCCGGGCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                                  <  1:450279/149‑1 (MQ=255)
                                                                              aacAGATACACCGGCGGAAAACGCGGCTTGGTTTTGGGGCGTGAAGGCGACCACGACCCAGCCCCCCCCCGCCCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                                  <  2:836966/149‑1 (MQ=255)
                                                                              aTCAGATTCACCAGCGAAACACGCGGCGTGTTTTTGGGGCGTGACTGCGCCCACGAGCCGGCCGCCCCCGCCCCCCACGACCATTTTTCCTGGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAc                                                                  <  1:1305090/151‑1 (MQ=255)
                                                                               taaGATCCACCGGCGGAAACCGCGGCATGTTTTTGGTCCGGAAATGCCCCCACGACCCAGCCGCCCCCCGCCCCCACGACCATTTTCCGGCTCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                                 <  1:1364377/149‑1 (MQ=255)
                                                                               tCAGTTTCATCGGCGGATACCGCGGCTTGGTTTTGGTGCGGAAATGCGCCCACGGCCCAGCCGCCCCCAGCCCCCCCGCCCATTTTTCCGTGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                                 <  2:280687/151‑1 (MQ=255)
                                                                               tCAGATTAACCAGCGAATAACGCGGCGTGTTTTTGGTGTGTGAAGGCGACCCCGAGCAAGACGCACCAGGCCCCCCCGCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                                 <  1:512608/151‑1 (MQ=255)
                                                                               tCAGATCCACCAGCGGATACCGCGCCATGTTTTTGGTGCGTGAATGCGACCACGAGCAAGACCCCCCCACCTCCCCCCCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                                 <  2:1628079/151‑1 (MQ=255)
                                                                               aaaGATTCACCAGCGAAAACCGGGGCTTGTTTTGGGTGGGTGAAGGCGCCCACGGCCCAGACGCCCCCAGCCCCCACGCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAcc                                                                 <  1:1283476/149‑1 (MQ=255)
                                                                                   tttcaccgggggAAAACGCGGCATGTTTTGGGTGGGTAAAGGCGCCCCCGAGCCAGACGCACCCGGCCCCCACACCCATTTTTCCGGGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGccc                                                             <  2:1609024/141‑1 (MQ=255)
                                                                                   tttCACCAGGGGAAAACGGGGCATGTTTTTGGTGCGTGAAGGCGACCCCGAGCCAGACCCCCCCAGCTCCCCCGACCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAAccc                                                                <  2:571026/147‑2 (MQ=255)
                                                                                       ctgaGCGAAAAACGCGGCAGGTTTTTGGTGCGTAAATGCGCCCACGCCCCAGCCGCCCCCAGCCCCCCCGCCCGTTTTCCCGTCCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAc                                                         <  1:834159/148‑1 (MQ=255)
                                                                                       cccgggggAAAACGGGGCATGTTTTTGGTGGGTGAAGGCGCCCACGGGCAAGACGCCCCCAGCCCCCACGACCATTTTCCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAc                                                         <  2:1360295/145‑1 (MQ=255)
                                                                                        tcagggaaaaccgcgGCGTGGTTTTGGGGGGGAAATGCGCCCACGACCGAGAGCCCCCCAGCCCCCACGCCCATTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAcc                                                        <  2:1467437/140‑1 (MQ=255)
                                                                                        ccaGCGGAAACCGGGGCATGGTTTTGGTGCGTAAATGCGACCCCGAGCCAGACGCCCCCAGCCCCCACGACCAGTTTTCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGAcc                                                        <  2:187801/150‑1 (MQ=255)
                                                                                           gCGAATACCGCGCCTTGTTTTTGGTGGGTGGAGGCGCCCACGAGCCAGACCCCCCAACCCCCCCCGCCCATTTCCCGGGGGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgag                                                     <  1:39521/150‑1 (MQ=255)
                                                                                            cGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagt                                                    >  2:1628035/1‑151 (MQ=255)
                                                                                             ggATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                                   >  1:1728244/1‑151 (MQ=255)
                                                                                              gaaaccgcgCCTTGTTTTTGGTGCGTAAATGCGCCCCCGCCCCAGCCGCCCCCACCCCCCCCGCCCATTTTTCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                                   <  2:489204/145‑1 (MQ=255)
                                                                                              gaaaaCGCGGCATGTTTTTGGTGGGGGAAGGCGCCCACGAGCCAGCCGCCCCCAGCCCCCACGACCATTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                                   <  2:344814/147‑1 (MQ=255)
                                                                                              gATAACGGGGCGTGGTTTGGGTGCGGAAATGCGACCACGGCCCAGACGCCCCCGGCCCCCCCGACCAGTTTTCCGGGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtg                                                   <  1:1353330/150‑1 (MQ=255)
                                                                                              gATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCgagtga                                                  >  2:1253151/1‑151 (MQ=255)
                                                                                                  aCGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAg                                              >  2:864074/1‑151 (MQ=255)
                                                                                                   cgcgGTGTGTTTTTGGGGCGGAAAGGCCCCCCCGCGCCACACGCCCCCAGCCCCCCCGCCCTGTTTTCCGGGCCCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGc                                             <  1:1015189/151‑1 (MQ=255)
                                                                                                     cgGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCa                                           >  2:950698/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                                    >  2:413793/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                                    >  2:1594543/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                                    >  2:163423/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                                    >  2:1551350/1‑151 (MQ=255)
                                                                                                             ttttGGTGCGTGAATGCGCCCACGACCAAGCCGCCCCCGGCCCCCACGACCAGTTTCCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                                   <  2:914662/151‑1 (MQ=255)
                                                                                                             tttgGGTGCGTGAATGGCCCCACGCGCCAGACCCACCCGGCCCCCCGGCCATTTTTCCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGagc                                   <  1:1175468/151‑3 (MQ=255)
                                                                                                             tttgGGGGCGTAAATGCCCCCACGGGCCCGCCCCCCCCAGCCCCCGGGACCTTTTTTCCGTCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                                   <  1:680621/151‑1 (MQ=255)
                                                                                                             ttgtGGTGCGTGAATGCCCCCACGAGCCAGACCCCCCCGGCCCCCCCGCCCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                                   <  1:1067711/148‑1 (MQ=255)
                                                                                                              tttggggggggaaaTGCCCCCACGACCCAGCCCCCCCCAGCCCCCACGACCAGTTTTCCGGGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATc                                   <  2:1455827/139‑1 (MQ=255)
                                                                                                              tttGGGGCGGGAAGGCGCCCACGAGCCAGACCCCCCCGGCCCCCACGACCATTTTCCCTGCCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCt                                  <  1:1358736/151‑1 (MQ=255)
                                                                                                                  gTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGt                              >  2:1032678/1‑151 (MQ=255)
                                                                                                                    gCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCg                            >  1:1246753/1‑151 (MQ=255)
                                                                                                                              cgcccACGACCCAGACGCCCCCAGCCCCCACGACCAGTTTCCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAg                   <  1:470186/148‑1 (MQ=255)
                                                                                                                                cccccggcccAGACCCCCCCACCCCCCACGACCATTTTCCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGAt                 <  2:1153290/143‑1 (MQ=255)
                                                                                                                                    cGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtg             >  1:730588/1‑151 (MQ=255)
                                                                                                                                      agccagccCCCCCCACCCCCCCCGCCCATTTTCCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtgt            <  2:89030/150‑1 (MQ=255)
                                                                                                                                      acccagccccccccAGCCCCCCCGCCCGTTTTCCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtgt            <  2:1325791/139‑1 (MQ=255)
                                                                                                                                      acccagccccccccAGCCCCCCCGACCGTTTTCCCGTGCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtgt            <  2:1244794/139‑1 (MQ=255)
                                                                                                                                       cccagccccccccGGCCCCCACGCCCATTTTCCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCAtgttg          <  2:746217/141‑1 (MQ=255)
                                                                                                                                          cgccccccccAGCCCCCCCGACCAGTTTCCCGTGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc       <  2:1482162/144‑1 (MQ=255)
                                                                                                                                          cgACCCACCCCGCCCCCAGAACCATTTTCCCGCGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc       <  1:362366/150‑1 (MQ=255)
                                                                                                                                           gccGCCCCCGGCCCCCCCGACCGTTTTTCGGTCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc       <  2:1707002/148‑1 (MQ=255)
                                                                                                                                            ccccccccACCCCCCCCGCCCATTTTCCCGGGCCCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc       <  2:228785/144‑1 (MQ=255)
                                                                                                                                               cccccAGCCCCCACGCCCGTTTTCCGGTGCACACCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGcc  <  1:187463/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTCAGGCC  >  NC_000913/3685878‑3686170

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