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

TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  NC_000913/3177603‑3177896
                                                                                                                                                  |                                                                                                                                                     
tCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCATGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGGATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCaaaa                                                                                                                                                   >  2:768819/1‑150 (MQ=255)
 cACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAAt                                                                                                                                                 >  1:892433/1‑151 (MQ=255)
   cccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAAttt                                                                                                                                               >  2:314255/1‑151 (MQ=255)
    ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑c                                                                                                                                             >  2:324462/1‑151 (MQ=255)
    ccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑c                                                                                                                                             <  2:1141375/151‑1 (MQ=255)
     ctctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑cc                                                                                                                                            >  2:1183957/1‑151 (MQ=255)
       ctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑Cccg                                                                                                                                          >  2:736140/1‑151 (MQ=255)
        ttCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑Cccgc                                                                                                                                         >  1:320121/1‑151 (MQ=255)
          cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑Cccgccg                                                                                                                                       >  2:846323/1‑151 (MQ=255)
            atcatcCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGGGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑Cccgcc                                                                                                                                        <  1:757079/148‑1 (MQ=255)
                  ccGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCAt                                                                                                                               >  1:1080708/1‑151 (MQ=255)
                        aCCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTc                                                                                                                          >  1:1010879/1‑150 (MQ=255)
                         ccATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCAt                                                                                                                        >  2:484404/1‑151 (MQ=255)
                               ccAGTAGCCAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGcc                                                                                                                   <  2:850238/150‑1 (MQ=255)
                               ccAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCATTCCGAATCTGCTCAATCAGCACAACTTCATCACGCCCTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAATTTT‑CCCGCCGCACCTCATGACTCAttt                                                                                                                      <  1:325788/147‑1 (MQ=255)
                                cAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCg                                                                                                                 >  1:238510/1‑151 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                                                                                         <  2:1046656/151‑1 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                           >  2:681111/1‑149 (MQ=255)
                                        aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCAAAACTTCATCACGCACGGGGTAAAAGGGAAGAAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg                                                                                                         <  2:896179/151‑1 (MQ=255)
                                          ggggTTTCGCTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACACCTCCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc                                                                                                       <  2:843078/151‑1 (MQ=255)
                                           gggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  1:51784/1‑151 (MQ=255)
                                            ggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  2:1020152/1‑150 (MQ=255)
                                            ggTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACg                                                                                                      >  2:542307/1‑150 (MQ=255)
                                              tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  2:289325/1‑151 (MQ=255)
                                              tttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  2:646855/1‑151 (MQ=255)
                                               ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                   >  1:453547/1‑150 (MQ=255)
                                                tCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                                 >  1:562699/1‑151 (MQ=255)
                                                  gctgtggccGTCCGCGCAAATCCGAATCTGCTCAATCAGCACACCTTCACCCCGCACGGGGTAAAAGGGTAGCAGGCCTGCGGG‑GTGCCCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                                  <  1:615155/140‑1 (MQ=255)
                                                  gcgggTGTCGTCCGCGGCAATCCGAATCTGCCCAATCAGCACACCTTCATCACGCACTGGGTCAAGGGGTAGCAAAACTGCGGC‑GTGACCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                <  1:496989/147‑1 (MQ=255)
                                                  gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                >  1:600381/1‑150 (MQ=255)
                                                  gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                >  1:425198/1‑150 (MQ=255)
                                                     ggTGTCGTACGCGGCAATCCGAATCTGTTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAAt                                                                                            <  2:1161264/151‑1 (MQ=255)
                                                      gtgtCGTACGCGGCAATCCGAATCTGCTCAACCAGCACACCTTCATCCCCCCCTGGGTAAAGGGGTAGAAAGACGGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATc                                                                                           <  1:970235/151‑1 (MQ=255)
                                                       tgtCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                          >  1:90575/1‑151 (MQ=255)
                                                        gtCGTTCGCGGCAATCCGAATCTCCTCAATCAGCACACCTTCATCACGCACTGGGCCAAAGGGTAGCAAGACGGCGGC‑GGGCCCGCGCCAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                          <  2:190000/150‑1 (MQ=255)
                                                           gcaCGGGGCAATCCGAATCTGCTAAATCAGCACAACTTCATCCCGACGTGGGCAAAAGGGTAGCAAGACGGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                       <  1:599459/148‑1 (MQ=255)
                                                           gTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata                                                                                       <  2:918080/150‑1 (MQ=255)
                                                             aCGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaaa                                                                                     <  2:957973/150‑1 (MQ=255)
                                                               gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAg                                                                                    <  2:677195/149‑1 (MQ=255)
                                                               gcgGCAATCCAAATCTGCTAAATCAGCACAACTCCATCACGCACTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGACCGTGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                                   <  1:60909/150‑1 (MQ=255)
                                                                 ggaaaTCCGAATCTGCTCAACCAGCACACCTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGGGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  1:436579/147‑1 (MQ=255)
                                                                 ggCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc                                                                                 <  2:656062/150‑1 (MQ=255)
                                                                     aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt                                                                             >  1:1030531/1‑150 (MQ=255)
                                                                      tCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg                                                                            >  2:653837/1‑150 (MQ=255)
                                                                       ccGAATCTGCTAAATCAGCCCACCTTCACCCCCCCCGGGGTAAAAGGGTAGCATTTCTGCGGC‑GTGCCCGCTCAAAAAATTTG‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt                                                                             <  1:516135/148‑1 (MQ=255)
                                                                          aaTCTGCTCAATCAGCCCAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                        <  2:303590/150‑1 (MQ=255)
                                                                          aaTCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                        >  1:15270/1‑150 (MQ=255)
                                                                               ccTCAATCAGCACACTTCCACCCCCCCCTGGTTCAAAGGGAAGAAATCCTGCGGC‑GTGACCGGGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa                                                                      <  2:1065681/146‑1 (MQ=255)
                                                                                ctaaaccAGCAAACCTTAACCACGCCCTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaaa                                                                     <  2:352666/141‑1 (MQ=255)
                                                                                   aaccAGCAAAACTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGa                                                              <  1:928176/148‑1 (MQ=255)
                                                                                   aaaagccacaACTTCATCCCTTTCTGGGCCAAAGGGTACCAAATCTGCGCC‑GGGACGGCG‑TAAAAAATTTGCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                               <  2:682800/144‑1 (MQ=255)
                                                                                   aaTCAGCACACCTTCATCACGCACTGGTTCAAGGGGTAGCAAGACTGCGGC‑GTGACCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgc                                                                <  1:434108/149‑1 (MQ=255)
                                                                                   aaTCAGCACAACTTCATCACGCACTGGGTCAAGGGGTAGAAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGa                                                              <  2:284469/151‑1 (MQ=255)
                                                                                    accAGCACAACTCCATCCCCCACGGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                             <  1:1014489/149‑1 (MQ=255)
                                                                                    aTCACCACAACTTCATCACGCACGGGGCCAAGGGGTAGCAAGCCTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                               <  2:652971/149‑1 (MQ=255)
                                                                                     tCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                             >  2:320760/1‑150 (MQ=255)
                                                                                       aGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                           >  1:1045996/1‑150 (MQ=255)
                                                                                            acctccctcgcGCTCTGGTTAAAAGGAAAACAAAACTGCGGT‑GACACCGCTCAAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGTGATAAAGTGtt                                                     <  2:1105182/141‑1 (MQ=255)
                                                                                            accTTCATCACCCACTGGGTAAAAGGGAAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaaa                                                          <  1:813537/144‑1 (MQ=255)
                                                                                            aaCTTCACCACGCCCTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                      <  1:393692/150‑1 (MQ=255)
                                                                                            aaCTTCACCACGCACTGGGTCAAGGGGTAGCAAGCCTGCGGC‑GTGCCCGCGCAAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt                                                        <  1:314254/148‑1 (MQ=255)
                                                                                             aCTTCATCCCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGtt                                                     <  2:167298/150‑1 (MQ=255)
                                                                                             aCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGttt                                                    >  1:1189085/1‑151 (MQ=255)
                                                                                               tccatcaCGCACTGGGTAAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt                                                      <  1:881157/145‑1 (MQ=255)
                                                                                                    cacgcacgGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                             <  1:191509/151‑1 (MQ=255)
                                                                                                    caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtc                                                   <  2:653242/145‑1 (MQ=255)
                                                                                                    caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                              <  2:320122/150‑1 (MQ=255)
                                                                                                     ccGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCAAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCg                                                <  1:467468/146‑1 (MQ=255)
                                                                                                      cGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtctc                                                 <  1:793093/145‑1 (MQ=255)
                                                                                                        cACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                          >  2:684006/1‑150 (MQ=255)
                                                                                                         attGGGTAAAAGGGTGCCAAGACT‑CGGCGGTTACCGCTTTAAAAATTTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                         <  2:297849/148‑1 (MQ=255)
                                                                                                                cAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTa                                  >  2:309653/1‑150 (MQ=255)
                                                                                                                 aaaGGGTAGCAAGTCTGCGGC‑GTGACCGCGCTAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                                 <  1:263820/150‑1 (MQ=255)
                                                                                                                 aaaGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                                 >  1:938196/1‑150 (MQ=255)
                                                                                                                    gggTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCg                             <  2:806423/151‑1 (MQ=255)
                                                                                                                    gggTAGCAAGACTGCGGC‑GTGACCGCGCTAAAAAATTC‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACAt                               <  1:540551/149‑1 (MQ=255)
                                                                                                                     ggTAGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGt                            <  2:29389/151‑1 (MQ=255)
                                                                                                                        aGCAAGACTGCGGC‑GTGCCCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtttt                         <  2:962874/151‑1 (MQ=255)
                                                                                                                           aaGACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTg                       <  2:956352/150‑1 (MQ=255)
                                                                                                                             gACTGCGGC‑GTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa                    <  2:453549/151‑1 (MQ=255)
                                                                                                                             gACTGCGGC‑GTGACCCCACTAAAAATTTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc                     <  1:732003/150‑1 (MQ=255)
                                                                                                                                 gcggc‑gTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCa                    <  1:1176742/147‑1 (MQ=255)
                                                                                                                                   ggc‑gTGCCCGCACTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa              <  2:1080715/151‑1 (MQ=255)
                                                                                                                                   ggc‑gTGACCGCGCTAAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATg                <  2:123784/149‑1 (MQ=255)
                                                                                                                                    gc‑gTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTa              >  2:427196/1‑150 (MQ=255)
                                                                                                                                          aCCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc        >  2:429390/1‑151 (MQ=255)
                                                                                                                                            cgcgCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGc        >  2:896998/1‑149 (MQ=255)
                                                                                                                                                cTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGtt   >  1:912002/1‑150 (MQ=255)
                                                                                                                                                cTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTg  >  1:405158/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGC‑GTGACCGCGCTCAAAAATTT‑CCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTG  >  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: ‑