Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,7480CA48.3% 72.1 / 68.2 89E55* (GAG→TAG) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (37/9);  new base A (0/43);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 7.45e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  NC_000913/3177603‑3177896
                                                                                                                                                  |                                                                                                                                                     
tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCATGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGGATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCaaaa                                                                                                                                                   >  2:768819/1‑150 (MQ=255)
 cACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAt                                                                                                                                                 >  1:892433/1‑151 (MQ=255)
   cccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAttt                                                                                                                                               >  2:314255/1‑151 (MQ=255)
    ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                             >  2:324462/1‑151 (MQ=255)
    ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTc                                                                                                                                             <  2:1141375/151‑1 (MQ=255)
     ctctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                                                                                            >  2:1183957/1‑151 (MQ=255)
       ctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccg                                                                                                                                          >  2:736140/1‑151 (MQ=255)
        ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgc                                                                                                                                         >  1:320121/1‑151 (MQ=255)
          cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                       >  2:846323/1‑151 (MQ=255)
            atcatcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGGGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                        <  1:757079/148‑1 (MQ=255)
                  ccGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCAt                                                                                                                               >  1:1080708/1‑151 (MQ=255)
                        aCCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTc                                                                                                                          >  1:1010879/1‑150 (MQ=255)
                         ccATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCAt                                                                                                                        >  2:484404/1‑151 (MQ=255)
                               ccAGTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc                                                                                                                   <  2:850238/150‑1 (MQ=255)
                               ccAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCATTCCGAATCTGCTCAATCAGCACAACTTCATCACGCCCTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAATTTTCCCGCCGCACCTCATGACTCAttt                                                                                                                      <  1:325788/147‑1 (MQ=255)
                                cAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCg                                                                                                                 >  1:238510/1‑151 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                                                                                         <  2:1046656/151‑1 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                           >  2:681111/1‑149 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCAAAACTTCATCACGCACGGGGTAAAAGGGAAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                                                                                         <  2:896179/151‑1 (MQ=255)
                                          ggggTTTCGCTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACACCTCCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                                                                                       <  2:843078/151‑1 (MQ=255)
                                           gggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  1:51784/1‑151 (MQ=255)
                                            ggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  2:1020152/1‑150 (MQ=255)
                                            ggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  2:542307/1‑150 (MQ=255)
                                              tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  2:289325/1‑151 (MQ=255)
                                              tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  2:646855/1‑151 (MQ=255)
                                               ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  1:453547/1‑150 (MQ=255)
                                                tCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                                 >  1:562699/1‑151 (MQ=255)
                                                  gctgtggccGTCCGCGCAAATCCGAATCTGCTCAATCAGCACACCTTCACCCCGCACGGGGTAAAAGGGTAGCAGGCCTGCGGGGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                                  <  1:615155/140‑1 (MQ=255)
                                                  gcgggTGTCGTCCGCGGCAATCCGAATCTGCCCAATCAGCACACCTTCATCACGCACTGGGTCAAGGGGTAGCAAAACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                <  1:496989/147‑1 (MQ=255)
                                                  gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                >  1:600381/1‑150 (MQ=255)
                                                  gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                >  1:425198/1‑150 (MQ=255)
                                                     ggTGTCGTACGCGGCAATCCGAATCTGTTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                                                                                            <  2:1161264/151‑1 (MQ=255)
                                                      gtgtCGTACGCGGCAATCCGAATCTGCTCAACCAGCACACCTTCATCCCCCCCTGGGTAAAGGGGTAGAAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc                                                                                           <  1:970235/151‑1 (MQ=255)
                                                       tgtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                          >  1:90575/1‑151 (MQ=255)
                                                        gtCGTTCGCGGCAATCCGAATCTCCTCAATCAGCACACCTTCATCACGCACTGGGCCAAAGGGTAGCAAGACGGCGGCGGGCCCGCGCCAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                          <  2:190000/150‑1 (MQ=255)
                                                           gcaCGGGGCAATCCGAATCTGCTAAATCAGCACAACTTCATCCCGACGTGGGCAAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                       <  1:599459/148‑1 (MQ=255)
                                                           gTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                       <  2:918080/150‑1 (MQ=255)
                                                             aCGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaaa                                                                                     <  2:957973/150‑1 (MQ=255)
                                                               gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAg                                                                                    <  2:677195/149‑1 (MQ=255)
                                                               gcgGCAATCCAAATCTGCTAAATCAGCACAACTCCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGTGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                                   <  1:60909/150‑1 (MQ=255)
                                                                 ggaaaTCCGAATCTGCTCAACCAGCACACCTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGGGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  1:436579/147‑1 (MQ=255)
                                                                 ggCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  2:656062/150‑1 (MQ=255)
                                                                     aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt                                                                             >  1:1030531/1‑150 (MQ=255)
                                                                      tCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                            >  2:653837/1‑150 (MQ=255)
                                                                       ccGAATCTGCTAAATCAGCCCACCTTCACCCCCCCCGGGGTAAAAGGGTAGCATTTCTGCGGCGTGCCCGCTCAAAAAATTTGCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt                                                                             <  1:516135/148‑1 (MQ=255)
                                                                          aaTCTGCTCAATCAGCCCAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                        <  2:303590/150‑1 (MQ=255)
                                                                          aaTCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                        >  1:15270/1‑150 (MQ=255)
                                                                               ccTCAATCAGCACACTTCCACCCCCCCCTGGTTCAAAGGGAAGAAATCCTGCGGCGTGACCGGGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                      <  2:1065681/146‑1 (MQ=255)
                                                                                ctaaaccAGCAAACCTTAACCACGCCCTGGGTAAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaaa                                                                     <  2:352666/141‑1 (MQ=255)
                                                                                   aaccAGCAAAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGa                                                              <  1:928176/148‑1 (MQ=255)
                                                                                   aaaagccacaACTTCATCCCTTTCTGGGCCAAAGGGTACCAAATCTGCGCCGGGACGGCGTAAAAAATTTGCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                               <  2:682800/144‑1 (MQ=255)
                                                                                   aaTCAGCACACCTTCATCACGCACTGGTTCAAGGGGTAGCAAGACTGCGGCGTGACCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgc                                                                <  1:434108/149‑1 (MQ=255)
                                                                                   aaTCAGCACAACTTCATCACGCACTGGGTCAAGGGGTAGAAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGa                                                              <  2:284469/151‑1 (MQ=255)
                                                                                    accAGCACAACTCCATCCCCCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                             <  1:1014489/149‑1 (MQ=255)
                                                                                    aTCACCACAACTTCATCACGCACGGGGCCAAGGGGTAGCAAGCCTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                               <  2:652971/149‑1 (MQ=255)
                                                                                     tCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                             >  2:320760/1‑150 (MQ=255)
                                                                                       aGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                           >  1:1045996/1‑150 (MQ=255)
                                                                                            acctccctcgcGCTCTGGTTAAAAGGAAAACAAAACTGCGGTGACACCGCTCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGTGATAAAGTGtt                                                     <  2:1105182/141‑1 (MQ=255)
                                                                                            accTTCATCACCCACTGGGTAAAAGGGAAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaaa                                                          <  1:813537/144‑1 (MQ=255)
                                                                                            aaCTTCACCACGCCCTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                      <  1:393692/150‑1 (MQ=255)
                                                                                            aaCTTCACCACGCACTGGGTCAAGGGGTAGCAAGCCTGCGGCGTGCCCGCGCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt                                                        <  1:314254/148‑1 (MQ=255)
                                                                                             aCTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGtt                                                     <  2:167298/150‑1 (MQ=255)
                                                                                             aCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                    >  1:1189085/1‑151 (MQ=255)
                                                                                               tccatcaCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                      <  1:881157/145‑1 (MQ=255)
                                                                                                    cacgcacgGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                             <  1:191509/151‑1 (MQ=255)
                                                                                                    caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                   <  2:653242/145‑1 (MQ=255)
                                                                                                    caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                              <  2:320122/150‑1 (MQ=255)
                                                                                                     ccGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCAAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCg                                                <  1:467468/146‑1 (MQ=255)
                                                                                                      cGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                 <  1:793093/145‑1 (MQ=255)
                                                                                                        cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                          >  2:684006/1‑150 (MQ=255)
                                                                                                         attGGGTAAAAGGGTGCCAAGACTCGGCGGTTACCGCTTTAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                         <  2:297849/148‑1 (MQ=255)
                                                                                                                cAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTa                                  >  2:309653/1‑150 (MQ=255)
                                                                                                                 aaaGGGTAGCAAGTCTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                                 <  1:263820/150‑1 (MQ=255)
                                                                                                                 aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                                 >  1:938196/1‑150 (MQ=255)
                                                                                                                    gggTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCg                             <  2:806423/151‑1 (MQ=255)
                                                                                                                    gggTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACAt                               <  1:540551/149‑1 (MQ=255)
                                                                                                                     ggTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                            <  2:29389/151‑1 (MQ=255)
                                                                                                                        aGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtttt                         <  2:962874/151‑1 (MQ=255)
                                                                                                                           aaGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTg                       <  2:956352/150‑1 (MQ=255)
                                                                                                                             gACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa                    <  2:453549/151‑1 (MQ=255)
                                                                                                                             gACTGCGGCGTGACCCCACTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc                     <  1:732003/150‑1 (MQ=255)
                                                                                                                                 gcggcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa                    <  1:1176742/147‑1 (MQ=255)
                                                                                                                                   ggcgTGCCCGCACTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa              <  2:1080715/151‑1 (MQ=255)
                                                                                                                                   ggcgTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATg                <  2:123784/149‑1 (MQ=255)
                                                                                                                                    gcgTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa              >  2:427196/1‑150 (MQ=255)
                                                                                                                                          aCCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc        >  2:429390/1‑151 (MQ=255)
                                                                                                                                            cgcgCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc        >  2:896998/1‑149 (MQ=255)
                                                                                                                                                cTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGtt   >  1:912002/1‑150 (MQ=255)
                                                                                                                                                cTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTg  >  1:405158/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  NC_000913/3177603‑3177896

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