Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA42.1% 120.6 / 36.8 76E384D (GAA→GATyhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (6/38);  new base A (32/0);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 2.05e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508514‑3508810
                                                                                                                                                     |                                                                                                                                                    
aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTc                                                                                                                                                    <  2:1049455/151‑1 (MQ=255)
aCGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTc                                                                                                                                                    <  2:362951/151‑1 (MQ=255)
       ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCacacac                                                                                                                                              <  1:1109749/150‑1 (MQ=255)
        tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACg                                                                                                                                             <  2:102826/150‑1 (MQ=255)
          tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAAttaaaaacgg                                                                                                                                             >  1:1130957/1‑141 (MQ=255)
          tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTacacaccg                                                                                                                                            >  2:331240/1‑147 (MQ=255)
           aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                          <  2:767221/150‑1 (MQ=255)
            aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGGAAAATCacacacaggag                                                                                                                                         >  1:393084/1‑145 (MQ=255)
             aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATtacaaaagggaacc                                                                                                                                       >  1:537067/1‑141 (MQ=255)
                 ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCca                                                                                                                                   <  1:972670/151‑1 (MQ=255)
                    gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                                <  2:475956/151‑1 (MQ=255)
                      gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACGGGAGcacccccctc                                                                                                                               >  2:306118/1‑145 (MQ=255)
                        tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACCGGAAccccccc                                                                                                                                  >  1:239068/1‑143 (MQ=255)
                        tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                                <  1:569819/147‑1 (MQ=255)
                            tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGaaa                                                                                                                          >  2:1075660/1‑149 (MQ=255)
                             atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                                       <  2:584656/151‑1 (MQ=255)
                                cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATTACACACCGTTGCCCCCCCCCCACAAAGCggg                                                                                                                     >  1:577352/1‑148 (MQ=255)
                                 aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACCGGAACCCCCCCCTCCGAAAGCgggg                                                                                                                    >  2:113648/1‑150 (MQ=255)
                                  gATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCgggggg                                                                                                                  >  2:117279/1‑149 (MQ=255)
                                   aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATAACAAAAGGGAGCCCCCACATCAGAAAGGgtg                                                                                                                     >  1:822560/1‑147 (MQ=255)
                                     gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGGCCCCACACCAGAAAACGGGGGgcg                                                                                                                >  2:402695/1‑150 (MQ=255)
                                     gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAAGGGAGCCCCCACCTCAGAAAgagggggggg                                                                                                                >  1:1078383/1‑141 (MQ=255)
                                             gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                       <  2:753606/151‑1 (MQ=255)
                                               gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACCGGAGCCCCCCCATCAGAAAGCGGGGGGGGCACAcc                                                                                                          >  1:998705/1‑146 (MQ=255)
                                               gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                                                                     <  1:274387/151‑1 (MQ=255)
                                                 gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGCAAAAATCAAAAACGGGAGCCCCCACAACAGAAAGCGGGGGGCGGACACCGAGccg                                                                                                    >  2:1129749/1‑148 (MQ=255)
                                                     gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                                <  2:178675/150‑1 (MQ=255)
                                                       aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                              <  1:452410/150‑1 (MQ=255)
                                                         cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCAAAAACGGTAGCCCCCACATCCGAAAACGGGGGGGGGCCACCCCGCTGCGGAGCggt                                                                                           >  1:614033/1‑149 (MQ=255)
                                                          ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGCGCGGCCGGGAAAAATGAAACACGGGAGCCCCCCCACACGAAAAAGCGGTTCGCACACCAACCTCCGttgagcgg                                                                                          >  1:872668/1‑144 (MQ=255)
                                                           cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                                                                          <  1:220164/150‑1 (MQ=255)
                                                            gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                                                                        <  2:1256650/151‑1 (MQ=255)
                                                                 aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                   <  1:455012/151‑1 (MQ=255)
                                                                    aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                                <  2:613962/151‑1 (MQ=255)
                                                                     aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGGGGGGGGCACACCGAGCTGCGTAGCGGTCATCCCTTTCTg                                                                               >  2:352616/1‑151 (MQ=255)
                                                                     aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                               <  1:1036101/151‑1 (MQ=255)
                                                                       aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                             <  2:30442/151‑1 (MQ=255)
                                                                          cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                            <  1:184831/149‑1 (MQ=255)
                                                                           aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACAGTAGACCCCCCATCAGAAAACGGGGTGGGCACACCGAGCTGCGGAGCGATCATCCCCTTCTGGTaaa                                                                          >  1:1147756/1‑150 (MQ=255)
                                                                                  gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                                   <  2:97104/150‑1 (MQ=255)
                                                                                     tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACCGGAGCCCCCACATCAGAAAACGTGGGGCGGACACCGAGCTGCGGAGCGGTCATCCCTTGCTGATAAAACAACACTgg                                                                >  1:1034997/1‑150 (MQ=255)
                                                                                      ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGGGGGGACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTgg                                                                >  1:1252030/1‑149 (MQ=255)
                                                                                      ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                                <  2:930826/149‑1 (MQ=255)
                                                                                       cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCCCCTCAGAAAACGGGGGGGGCACAACGAGCTGCGGAGCGATCATCCCCTTCTGGTAAAACAGCACCgg                                                                >  1:987222/1‑148 (MQ=255)
                                                                                       cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACCGGAGCCCCCCCAACAGAAAAAGTGGGGGGCACACCGCGCGGCGGAGCGATCACACCCTGCTGATAAAACACCACCGGgac                                                             >  2:1027285/1‑151 (MQ=255)
                                                                                        tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCAAAACCCGCAACCCCCACATCAAAAAGCGGGGGGGGGACACCTAGCTGCGGGGGGACCATCCCCTTCTGGTAAAACAacacccgcccc                                                            >  2:32959/1‑144 (MQ=255)
                                                                                         gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac                                                             <  2:922745/149‑1 (MQ=255)
                                                                                           aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAAATCACACACGGTAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAACACTGGcaccac                                                          >  2:276108/1‑150 (MQ=255)
                                                                                           aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACAACAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAACCAACACTGGcaccac                                                          >  2:1189993/1‑150 (MQ=255)
                                                                                             aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                       <  2:543235/151‑1 (MQ=255)
                                                                                              ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg                                                      <  1:680176/151‑1 (MQ=255)
                                                                                                ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGTGCCCGGGAAAAATCACACACGGGAACCCCCACAACACAAAACGGGGGGCGCCCACCCAGCTGCGCAGCGGTCATCCCCTGCTGATAAAACAACACCGGCCCCACCTcgcg                                                    >  2:1232371/1‑151 (MQ=255)
                                                                                                  agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGCAGCCCCCACATCAGAAAGGGGGGTGCGCACACCAGGCTGCGTTGCGGTCATCCCCTTCTGCTAAAACAGCCCCGGCACCCCTTCGCggg                                                  >  1:472119/1‑151 (MQ=255)
                                                                                                     gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGGAGCCCCCCCATCCAAAAACGGGGTGCGCACACCCAGCCGCGTAGCGATCATCCCCTTCCGATAAACCAGCACTGGGCCCACCTCCCggg                                                  >  2:1026897/1‑148 (MQ=255)
                                                                                                      aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACAACGGAGCCCCCCCATCAGAAAACGTGGGGCGGACACCCAGCTGCGTAGCGATCATCCCCTTCTGATAAAACAGCACTGGGCCCACCTCGCGGGGAt                                               >  2:510462/1‑150 (MQ=255)
                                                                                                      aggaAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCGCGGAAAAAATAAAACACCGGAGACCCCCCCACAAAAAAAGGGGGGCGCACACCCAGCTGCGTAGGGATCATCCCCTTCTTATAAACCAGCACCGGCACCACCTCGCcggg                                                 >  2:198601/1‑147 (MQ=255)
                                                                                                             aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAAACCGTAGCCCCCCCCTCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGGGATCATCCCTTTCTGGTAAACAAACACCGGCACCACTTCGCGGGGATGGGGGGtg                                       >  1:275055/1‑151 (MQ=255)
                                                                                                              aaagaGCGACCATTTTGTGGCGGTGGGGGGGGGGAAAATTTACACACCGGAGCCCCCCCATCCAAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAAAAAGCACTGGCCCCCCTTCGCGGGGAtggtgg                                          >  1:311559/1‑147 (MQ=255)
                                                                                                               aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                                     <  1:379642/151‑1 (MQ=255)
                                                                                                                agagCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACCTCGCGGGGAtggtgg                                          >  2:1029380/1‑145 (MQ=255)
                                                                                                                 gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                                   <  2:393090/151‑1 (MQ=255)
                                                                                                                  agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCAACAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTTCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGTGGGTGTGGcc                                  >  2:471885/1‑151 (MQ=255)
                                                                                                                   gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                 <  1:1007853/151‑1 (MQ=255)
                                                                                                                         aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc                           <  1:402231/151‑1 (MQ=255)
                                                                                                                              gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                       <  2:750283/150‑1 (MQ=255)
                                                                                                                              gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                       <  2:1119380/150‑1 (MQ=255)
                                                                                                                              gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                      <  1:39309/151‑1 (MQ=255)
                                                                                                                                ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                    <  1:289274/151‑1 (MQ=255)
                                                                                                                                ggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                    <  1:597648/151‑1 (MQ=255)
                                                                                                                                  cggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                  <  1:700206/151‑1 (MQ=255)
                                                                                                                                          cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc          >  1:411692/1‑151 (MQ=255)
                                                                                                                                                aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc      >  1:1112375/1‑150 (MQ=255)
                                                                                                                                                 aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt    >  1:479741/1‑151 (MQ=255)
                                                                                                                                                 aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt    >  1:848355/1‑151 (MQ=255)
                                                                                                                                                   aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa  <  1:14621/151‑1 (MQ=255)
                                                                                                                                                   aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa  <  2:714172/151‑1 (MQ=255)
                                                                                                                                                   aaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCa  <  2:977637/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCA  >  NC_000913/3508514‑3508810

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