Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,686,0210AC45.1% 135.9 / 15.0 91G1051G (GGT→GGGbcsCcellulose synthase subunit
Reads supporting (aligned to +/- strand):  ref base A (30/20);  new base C (0/41);  total (30/61)
Fisher's exact test for biased strand distribution p-value = 6.88e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTC  >  NC_000913/3685878‑3686165
                                                                                                                                                |                                                                                                                                                   
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:362482/1‑151 (MQ=255)
      gtgtAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACg                                                                                                                                        >  2:555296/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)
                          gctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGc                                                                                                                                               <  2:116838/124‑1 (MQ=255)
                          gctgcCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGc                                                                                                                                               >  1:116804/1‑124 (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:1285119/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:322222/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:1363349/1‑151 (MQ=255)
                                                      gCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACa                                                                                     >  1:665905/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:163423/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:1594543/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:1551350/1‑151 (MQ=255)
                                                                                                            gTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGAt                               >  2:413793/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)
                                                                                                                                          cgACCCACCCCGCCCCCAGAACCATTTTCCCGCGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc  <  1:362366/150‑1 (MQ=255)
                                                                                                                                           gccGCCCCCGGCCCCCCCGACCGTTTTTCGGTCCGCCCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTc  <  2:1707002/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGCGCCGTGTAGCCGAAGCCCTGACTGCTGCCGCCATCGTTGGATTGCCGCGCAGCTTCTTCCTGCCAGTCGGTCGGGATCAGATTCATCAGCGGATAACGCGGCATGGTTTTGGTGCGTGAATGCGACCACGAGCCAGACGCACCCAGCTCCCACGACCAGTTTTCCGTGCGCTCCCGCCACATCACCGGTATGGCAAACGACAGGTATTCCTGCGGACTGTAGTAACCGCCCTGACCGAGTGAGTAGCCACTCAGATCTTTGTCGTAATGCCAGATCATGTTGTTC  >  NC_000913/3685878‑3686165

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: