Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG42.9% 108.5 / 19.4 71F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (5/35);  new base G (30/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.17e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508509‑3508806
                                                                                                                                                      |                                                                                                                                                     
aGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGaaa                                                                                                                                                     <  2:913262/151‑1 (MQ=255)
   ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGgaaaa                                                                                                                                                   >  2:67868/1‑150 (MQ=255)
     aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTc                                                                                                                                                <  2:1049455/151‑1 (MQ=255)
     aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTc                                                                                                                                                <  2:362951/151‑1 (MQ=255)
       gttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGCGG‑GGaaaaa                                                                                                                                                   >  1:750274/1‑147 (MQ=255)
            ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCacacac                                                                                                                                          <  1:1109749/150‑1 (MQ=255)
             tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACg                                                                                                                                         <  2:102826/150‑1 (MQ=255)
               tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAAAttaaaaacgg                                                                                                                                         >  1:1130957/1‑141 (MQ=255)
               tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATTacacaccg                                                                                                                                        >  2:331240/1‑147 (MQ=255)
                aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTa                                                                                                                                      <  2:767221/150‑1 (MQ=255)
                 aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GGGC‑GGGGAAAATCacacacaggag                                                                                                                                     >  1:393084/1‑145 (MQ=255)
                  aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATtacaaaagggaacc                                                                                                                                   >  1:537067/1‑141 (MQ=255)
                      ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCca                                                                                                                               <  1:972670/151‑1 (MQ=255)
                         gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                            <  2:475956/151‑1 (MQ=255)
                           gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATAACACACGGGAGcacccccctc                                                                                                                           >  2:306118/1‑145 (MQ=255)
                             tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTACAAACCGGAAccccccc                                                                                                                              >  1:239068/1‑143 (MQ=255)
                             tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                            <  1:569819/147‑1 (MQ=255)
                                 tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATCACACACGGGAGCCCCCACATCAGaaa                                                                                                                      >  2:1075660/1‑149 (MQ=255)
                                  atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                   <  2:584656/151‑1 (MQ=255)
                                     cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGG‑GCGC‑GGGAAAAATTACACACCGTTGCCCCCCCCCCACAAAGCggg                                                                                                                 >  1:577352/1‑148 (MQ=255)
                                      aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTAAACACCGGAACCCCCCCCTCCGAAAGCgggg                                                                                                                >  2:113648/1‑150 (MQ=255)
                                       gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCgggggg                                                                                                              >  2:117279/1‑149 (MQ=255)
                                        aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GGGC‑GGGAAAAATAACAAAAGGGAGCCCCCACATCAGAAAGGgtg                                                                                                                 >  1:822560/1‑147 (MQ=255)
                                          gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTACACACGGTAGGCCCCACACCAGAAAACGGGGGgcg                                                                                                            >  2:402695/1‑150 (MQ=255)
                                          gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTACACAAGGGAGCCCCCACCTCAGAAAgagggggggg                                                                                                            >  1:1078383/1‑141 (MQ=255)
                                                  gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                   <  2:753606/151‑1 (MQ=255)
                                                    gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACAAACCGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACAcc                                                                                                      >  1:998705/1‑146 (MQ=255)
                                                    gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                 <  1:274387/151‑1 (MQ=255)
                                                      gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGCAAAAATCAAAAACGGGAGCCCCCACAACAGAAAGCGGGGGGCGGACACCGAGccg                                                                                                >  2:1129749/1‑148 (MQ=255)
                                                          gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                            <  2:178675/150‑1 (MQ=255)
                                                            aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                          <  1:452410/150‑1 (MQ=255)
                                                              cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATCAAAAACGGTAGCCCCCACATCCGAAAACGGGGGGGGGCCACCCCGCTGCGGAGCggt                                                                                       >  1:614033/1‑149 (MQ=255)
                                                               ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGCGC‑GGCC‑GGGAAAAATGAAACACGGGAGCCCCCCCACACGAAAAAGCGGTTCGCACACCAACCTCCGttgagcgg                                                                                      >  1:872668/1‑144 (MQ=255)
                                                                cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                      <  1:220164/150‑1 (MQ=255)
                                                                 gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                    <  2:1256650/151‑1 (MQ=255)
                                                                      aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                               <  1:455012/151‑1 (MQ=255)
                                                                         aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                            <  2:613962/151‑1 (MQ=255)
                                                                          aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGGGGGGGGCACACCGAGCTGCGTAGCGGTCATCCCTTTCTg                                                                           >  2:352616/1‑151 (MQ=255)
                                                                          aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                           <  1:1036101/151‑1 (MQ=255)
                                                                            aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                         <  2:30442/151‑1 (MQ=255)
                                                                               cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                        <  1:184831/149‑1 (MQ=255)
                                                                                aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACACACAGTAGACCCCCCATCAGAAAACGGGGTGGGCACACCGAGCTGCGGAGCGATCATCCCCTTCTGGTaaa                                                                      >  1:1147756/1‑150 (MQ=255)
                                                                                       gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                               <  2:97104/150‑1 (MQ=255)
                                                                                          tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATAACACACCGGAGCCCCCACATCAGAAAACGTGGGGCGGACACCGAGCTGCGGAGCGGTCATCCCTTGCTGATAAAACAACACTgg                                                            >  1:1034997/1‑150 (MQ=255)
                                                                                           ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGGGGGGACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTgg                                                            >  1:1252030/1‑149 (MQ=255)
                                                                                           ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                            <  2:930826/149‑1 (MQ=255)
                                                                                            cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTACACACCGTAGCCCCCCCCTCAGAAAACGGGGGGGGCACAACGAGCTGCGGAGCGATCATCCCCTTCTGGTAAAACAGCACCgg                                                            >  1:987222/1‑148 (MQ=255)
                                                                                            cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAAATCACACACCGGAGCCCCCCCAACAGAAAAAGTGGGGGGCACACCGCGCGGCGGAGCGATCACACCCTGCTGATAAAACACCACCGGgac                                                         >  2:1027285/1‑151 (MQ=255)
                                                                                             tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGG‑GCGC‑GGGAAAAATCAAAACCCGCAACCCCCACATCAAAAAGCGGGGGGGGGACACCTAGCTGCGGGGGGACCATCCCCTTCTGGTAAAACAacacccgcccc                                                        >  2:32959/1‑144 (MQ=255)
                                                                                              gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                                                         <  2:922745/149‑1 (MQ=255)
                                                                                                aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGAAAAAATCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAACACTGGcaccac                                                      >  2:276108/1‑150 (MQ=255)
                                                                                                aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACACACGGGAGCCCCCACAACAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAACCAACACTGGcaccac                                                      >  2:1189993/1‑150 (MQ=255)
                                                                                                  aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                   <  2:543235/151‑1 (MQ=255)
                                                                                                   ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg                                                  <  1:680176/151‑1 (MQ=255)
                                                                                                     ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGT‑GCCC‑GGGAAAAATCACACACGGGAACCCCCACAACACAAAACGGGGGGCGCCCACCCAGCTGCGCAGCGGTCATCCCCTGCTGATAAAACAACACCGGCCCCACCTcgcg                                                >  2:1232371/1‑151 (MQ=255)
                                                                                                       agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATCACACACGGCAGCCCCCACATCAGAAAGGGGGGTGCGCACACCAGGCTGCGTTGCGGTCATCCCCTTCTGCTAAAACAGCCCCGGCACCCCTTCGCggg                                              >  1:472119/1‑151 (MQ=255)
                                                                                                          gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCAAACACGGGAGCCCCCCCATCCAAAAACGGGGTGCGCACACCCAGCCGCGTAGCGATCATCCCCTTCCGATAAACCAGCACTGGGCCCACCTCCCggg                                              >  2:1026897/1‑148 (MQ=255)
                                                                                                           aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGC‑GGGAAAAATTACACAACGGAGCCCCCCCATCAGAAAACGTGGGGCGGACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTGGGCCCACCTCGCGGGGAt                                           >  2:510462/1‑150 (MQ=255)
                                                                                                           aggaAAGAAAAGAGCGACCATTTTGTGGCGGCGGG‑GCGC‑GGAAAAAATAAAACACCGGAGACCCCCCCACAAAAAAAGGGGGGCGCACACCCAGCTGCGTAGGGATCATCCCCTTCTTATAAACCAGCACCGGCACCACCTCGCcggg                                             >  2:198601/1‑147 (MQ=255)
                                                                                                                  aaaagaGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATTAAAAACCGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGGGATCATCCCTTTCTGGTAAACAAACACCGGCACCACTTCGCGGGGATGGGGGGtg                                   >  1:275055/1‑151 (MQ=255)
                                                                                                                   aaagaGCGACCATTTTGTGGCGGTGGG‑GGGG‑GGGAAAATTTACACACCGGAGCCCCCCCATCCAAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAAAAGCACTGGCCCCCCTTCGCGGGGAtggtgg                                      >  1:311559/1‑147 (MQ=255)
                                                                                                                    aagaGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                 <  1:379642/151‑1 (MQ=255)
                                                                                                                     agagCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACCTCGCGGGGAtggtgg                                      >  2:1029380/1‑145 (MQ=255)
                                                                                                                      gagCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                               <  2:393090/151‑1 (MQ=255)
                                                                                                                       agCGACCATTTTGTGGCGGTGGC‑GCGC‑GGGAAAAATCACACACGGTAGCCCCCCCAACAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGTGGGTGTGGcc                              >  2:471885/1‑151 (MQ=255)
                                                                                                                        gCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  1:1007853/151‑1 (MQ=255)
                                                                                                                              aTTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                       <  1:402231/151‑1 (MQ=255)
                                                                                                                                   gtggcggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                   <  2:750283/150‑1 (MQ=255)
                                                                                                                                   gtggcggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                   <  2:1119380/150‑1 (MQ=255)
                                                                                                                                   gtggcggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                  <  1:39309/151‑1 (MQ=255)
                                                                                                                                     ggcggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                <  1:289274/151‑1 (MQ=255)
                                                                                                                                     ggcggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                <  1:597648/151‑1 (MQ=255)
                                                                                                                                       cggtggc‑gCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg              <  1:700206/151‑1 (MQ=255)
                                                                                                                                                cgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc      >  1:411692/1‑151 (MQ=255)
                                                                                                                                                      aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  1:1112375/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508509‑3508806

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