Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,177,5310AC44.9% 133.3 / 18.9 91V127G (GTG→GGG) nudFADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (32/17);  new base C (0/40);  total (32/59)
Fisher's exact test for biased strand distribution p-value = 4.85e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATC  >  NC_000913/3177384‑3177674
                                                                                                                                                     |                                                                                                                                                
gCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACAcc                                                                                                                                                <  2:767913/150‑1 (MQ=255)
    gTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGttt                                                                                                                                           <  1:835005/151‑1 (MQ=255)
    gTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGtt                                                                                                                                            <  1:63324/150‑1 (MQ=255)
       cccGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                                                                                                        <  2:284963/151‑1 (MQ=255)
         cGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTcc                                                                                                                                      <  2:1134876/151‑1 (MQ=255)
          ggCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                                                                                                     <  2:242398/151‑1 (MQ=255)
          ggCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                                                                                                     >  2:911928/1‑151 (MQ=255)
           gCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGt                                                                                                                                    <  2:47134/151‑1 (MQ=255)
            cTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                                                                                                   >  1:478755/1‑151 (MQ=255)
             ttACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGttt                                                                                                                                  >  1:213874/1‑151 (MQ=255)
               aCCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                                                                                                                 <  2:904455/150‑1 (MQ=255)
                   caTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTa                                                                                                                             <  1:105661/150‑1 (MQ=255)
                    aTGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATc                                                                                                                           >  2:762823/1‑151 (MQ=255)
                      gAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                                                                                                         <  2:489843/150‑1 (MQ=255)
                       aaCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                                                                                                       >  1:905102/1‑151 (MQ=255)
                       aaCGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                                                                                                       >  1:862833/1‑151 (MQ=255)
                          gcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                                                                                                    <  2:5773/151‑1 (MQ=255)
                          gcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGATTGACTATCAGTcc                                                                                                                     >  2:936120/1‑150 (MQ=255)
                           cgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                                                                                                    <  1:1092640/150‑1 (MQ=255)
                           cgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCg                                                                                                                   >  2:430051/1‑151 (MQ=255)
                              atatCTTCGTTTTCATCAGCCAGCCCGTGAATACCGCTTGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCct                                                                                                                <  1:462241/151‑1 (MQ=255)
                                atCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                                                                                                              >  2:359235/1‑151 (MQ=255)
                                    tCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCa                                                                                                           >  1:796269/1‑150 (MQ=255)
                                       ttttCACCCCCCAGCCCGTAAATCCCGCTTGCGGTCGGGGCGCCCCCTTCGCCCCCCATAATTGACGAACCCTCACTGGTCCCCCCCGGGCTTGCCAGGAAACTTAACCCGGTTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAAt                                                                                                         <  1:81689/149‑1 (MQ=255)
                                       ttttCACCACCAAGCCCGTGAATCCCGTTGGGGGTGTGGGCGCCCCCTCCCCCCCCAATATTTAACAAACCCTCGCTGGTCCCCCCCGGGCTGGCCAGAAACCTTACCCCCGGTTTGGTCCTTTTGACTATCAGTCCCGCctct                                                                                                              <  1:1006313/144‑1 (MQ=255)
                                            accgccgagcccttGAATCCCGCTGGCGGGCGTGGCCTCCACCCCGCCCACCATATTTGACGACGCCTCACTGGTGCCCCCGGGGTTGGCCAGGATATTGAACCCCGTTTCGGTCCGTTTAACTATCAGGCCCGCCTCTTCAATCGCTTc                                                                                                   <  2:878507/137‑1 (MQ=255)
                                             ccgccaaaaCCGTGAACCCCGCGTGGGGTGGGGGGGTCAACTCCCCCCACAATAATTGCCGAGCCCTCTTGGGTGCCCCCCGGGCTTGCTAGGAAACTTACCCCCGGTTGGTTCGGTTTGTCGAGCCGTCCCGCCTCTTCAATCGCTTcg                                                                                                  <  1:969815/142‑1 (MQ=255)
                                              ccgCCAGCCCGTAAATCCCGTTTGGGGCGGGGGCCCCCCCTCCCCCCCCCATAATTGACGACCGCTCACTGGGCCCCCCCGGGTTGGCCAGGAAACTTAACCCCGTTTTGGTCCGTTTGACTATAAGTCCCGCCTCTTCAATCGCt                                                                                                     <  1:688802/144‑1 (MQ=255)
                                                gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa                                                                                               >  2:832026/1‑150 (MQ=255)
                                                 ccaccccGTAAATCCCGTTGGCGGTCGTGGCGCCCCCTCCGCCCCCCATAATTGACGACCGCTCACTGGTCCCCCCCGGGTTTGCCAGAAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcg                                                                                                  <  1:855816/141‑1 (MQ=255)
                                                 ccAGCCCGTGAATACCGCTTGCGGTCGTGGCGCCCCCTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACg                                                                                             <  2:827288/151‑1 (MQ=255)
                                                 ccAGCCCGTGAACCCCGCTTGCGGTGGGGGCCCCCACTCCCCCCCCCATATTGGACGACCCCTCCCTGGTCCCCCCCGGGCTTCCCAGGAAACTAAACCCCGTTTTGGTCCTTTGGACTATCAGCCCCGCCTCTTCAATCGCTTcgcg                                                                                                <  1:752668/148‑1 (MQ=255)
                                                 ccAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACg                                                                                             >  2:984497/1‑151 (MQ=255)
                                                   aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACggg                                                                                           >  2:897663/1‑151 (MQ=255)
                                                      ccGTGAATCCCGCTTGCGCTCGGGGCGCCCCCTTCCCCCCCCAAAATGGCCAAACGCTCTCTGGTGCCCCCCGGCCTCGCCAGGAAACTTAACCCCGTTTGGGTCGGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgc                                                                                                 <  1:474215/142‑1 (MQ=255)
                                                      ccGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCa                                                                                         >  1:262564/1‑150 (MQ=255)
                                                         tGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACAt                                                                                     >  1:258181/1‑151 (MQ=255)
                                                              cccGCTTGCGGTCGGGGCGCCCCCTCCCCCCCCCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGAAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACAt                                                                                     <  1:617899/145‑1 (MQ=255)
                                                               cggttggCGGCCGGGGCCTCCCCTCCCCCCCCAAAAATTGACGAACGCTCCCGGGGCCCCCCGGGGCTTCCAAGAAAACTTACCCCCGGTTTGGCCCGTTTGACATCAAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAAc                                                                                       <  2:1067239/137‑1 (MQ=255)
                                                                  cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCacac                                                                            >  1:6595/1‑151 (MQ=255)
                                                                   tggggggcggggCGCCCCCTCCGCCCCCCAAAATGGACGACCGCTCACTGGCGCCCCCCGGGCTTGCCAGAAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCcaca                                                                             <  1:415013/139‑1 (MQ=255)
                                                                     gcggcgtggccgcccACTCCCCCCCCCATAATTGACAACCCCGCAGGGGTGCCCCCGGGGTTTGCCAGAAACCTTAACCCCGTTTGGGTCGGTTTACCTTTCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCcaca                                                                             <  1:548013/134‑1 (MQ=255)
                                                                       ggtgggggCGCCCCCTCCCCCCCCCAAAATGGACGACCGCTCCCGGGTGCCCCCCGGGTTTGCCAGAAAACTTAACCCCGGTTTGGTCCGTTTGACTATCATTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCacac                                                                            <  1:269983/140‑1 (MQ=255)
                                                                       ggtgggggCGCCCACTCCGCCCCCAAAATTCGACCACCGCTCCCTGGTCCCCCCCGGGCTTGACAGAAAACTTAACTCCGGTTTGTTCTGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCa                                                                       <  2:690496/145‑1 (MQ=255)
                                                                       ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTc                                                                        >  1:136008/1‑150 (MQ=255)
                                                                       ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTc                                                                        >  2:126183/1‑150 (MQ=255)
                                                                       ggTCGGGGCCCCCCCTTCCCCCACCATAATTGACGACCGCTCACTGTTCCCCCCCGGGCTTGCCAGAAAACTTACCCCCTGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTc                                                                        <  2:191908/150‑1 (MQ=255)
                                                                        gccGTGGCCTCCTCTTCCCCCCCCATAATTGACGAACGCTCGCTGGTGCCCCCGGGGCTTGCCAGGAATCTTATCCCCGGTTTGTCCCGTTTGACTATAACTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCa                                                                       <  2:501606/148‑1 (MQ=255)
                                                                           gTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                    >  1:1202424/1‑150 (MQ=255)
                                                                               cgccccctcccccccccATAATTGACGAACGCTCCCTGGTCCCCCCCGGGCTTGCCAGAAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCct                                                                   <  2:862871/136‑1 (MQ=255)
                                                                               cGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTc                                                               >  2:656839/1‑151 (MQ=255)
                                                                                gTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTc                                                               >  1:834654/1‑150 (MQ=255)
                                                                                 cccgctccccccACAATAATTGACGAACCCTCCCTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCaa                                                             <  1:662941/142‑1 (MQ=255)
                                                                                  cccctcccccccccATAATTGACGAACCCTCCCGGGTCCCCCCCGGGCTTGCAAGGAAATTTACCCCCGGTTTGGCCCTTTTGACTACCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTc                                                               <  1:1133640/140‑1 (MQ=255)
                                                                                  cccctccccccACCATAATTGACGACCGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCtt                                                                <  2:1249249/139‑1 (MQ=255)
                                                                                  cccctccGCCCCCCATATTGGCCGAACCCGCCGTGGTGCCCCCGGGGCTGGCCAGAAAACTTAACCCCGGTTGGGTCCTTTTGCCTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCa                                                              <  1:164621/143‑1 (MQ=255)
                                                                                  ccactcccccacccaTAATTGACAAACGCTTAGTGGTGCCCCCCGCACTGCGCAGCAAAATTAACACCGGTTTGGCCCTTTTCCTTTCCCCTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACGCTTTCACCCTCtt                                                                <  1:604693/135‑1 (MQ=255)
                                                                                  ccACTCCGCCCACAATAATAGCCGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGAAATTTAAACACCGTTTTGGACCGTTCGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCctct                                                                 <  2:1006386/146‑1 (MQ=255)
                                                                                   cACTCCCCCCCCCAAAATTGCCGACCCCTCACGGTCCCCCCCCGGGTTTCCCAGGAATCTCACCCCGGGTTTGGCCCGTTTGTTTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                    <  1:242775/142‑1 (MQ=255)
                                                                                      tcccccccccAAAATGGACGACCCCTCCCGGGTGCCCCCCGGGTTTGCCAGAAAATTAACCCCCGTTTTGGTCGTTTTAATTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatca                                                          <  1:1043708/145‑1 (MQ=255)
                                                                                       ccccccACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGTTTGCCAGAAAACTTAACCCCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATcc                                                       <  2:639969/148‑1 (MQ=255)
                                                                                        ccccccccAAAATTGACAAACGCTCACTGGTCCCCCCCGGGCTTGCCAGAAAACTTAACCCCGTTTGGGCCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatcatc                                                        <  1:841679/147‑1 (MQ=255)
                                                                                         cccccccATAATGGACGACCGCTCACGGGTCCCCCCCGGTCTTGCCAGAAAACTTAACCCCGTTTTGGTCCGTTTGCCTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCg                                                     <  1:874021/150‑1 (MQ=255)
                                                                                            ccccaTAATGGCGGAACGCTCACTGGTGCCCCCCGGTTTTGCCAGGAAACTAACCCCCGGTTTGGCCCTTTTGTCTTTCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCa                                                  <  2:1200739/149‑1 (MQ=255)
                                                                                                aTAATTGACGACCCCTCACTGTTCCCCCCCGGGCTTGCCAGGAAACTTACCCCCGGTTTGGCCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAAcc                                               <  2:1263185/150‑1 (MQ=255)
                                                                                                     tGACGGAAGCTCACCGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGGTTTGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCCTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTcc                                         >  1:574975/1‑151 (MQ=255)
                                                                                                       gACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTcc                                         >  1:630712/1‑150 (MQ=255)
                                                                                                        aCGAACGCTCCCTGGTGCCCCCCGGGCTTGCCAGGAAACTTACCCCCGGTTGGGCCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCa                                        <  2:875233/150‑1 (MQ=255)
                                                                                                           aaCGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAg                                    <  2:361658/151‑1 (MQ=255)
                                                                                                            aCGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTACCCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGt                                      <  1:1181571/148‑1 (MQ=255)
                                                                                                             cGCTCCGTGGTCCCCCCCGGGCTTGCCAGGAAACTTACCCCCGTTTGGGCCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGcc                                  <  2:149033/151‑1 (MQ=255)
                                                                                                              gCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGcc                                  >  1:1013452/1‑150 (MQ=255)
                                                                                                               cTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCa                                 >  2:537336/1‑150 (MQ=255)
                                                                                                               cTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCa                                 >  2:34086/1‑150 (MQ=255)
                                                                                                               cTCACTGGTCCCCCCCGGCCTTGCCAGAAAACTTAACCCGGGTTGGGTCCTTTTGACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAg                                       <  1:207773/144‑1 (MQ=255)
                                                                                                                 cccgggTCCCCCCCGGGCTTCCCAGAAACTTAACCCCCGTTTTGGCCCTTTTGACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCaa                                <  1:1146196/145‑1 (MQ=255)
                                                                                                                 cACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAgg                              >  1:121447/1‑151 (MQ=255)
                                                                                                                 cACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAgg                              >  1:244603/1‑151 (MQ=255)
                                                                                                                     ggTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAgggg                            <  2:1145598/149‑1 (MQ=255)
                                                                                                                     ggTGCCCCCCGGGCTTGCCAGAAAACTTAACACCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAggg                             <  2:1212565/148‑1 (MQ=255)
                                                                                                                            cccgtcttggcgaGGAAACTTACCGCCGGTGTGTTCTGTTTGACTTTCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCTCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTgg                   <  1:54717/139‑1 (MQ=255)
                                                                                                                                ggCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTc               >  2:224813/1‑151 (MQ=255)
                                                                                                                                   ttGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGt             >  1:745018/1‑150 (MQ=255)
                                                                                                                                   ttGCCAGAAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGt             <  2:1178831/150‑1 (MQ=255)
                                                                                                                                    tGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTa            >  1:407989/1‑150 (MQ=255)
                                                                                                                                      ccAGGAACCTAAACCCCGTTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTAcgc         <  2:461410/151‑1 (MQ=255)
                                                                                                                                      ccAGGAAACTTACCACCGTTTGGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTAcgc         <  2:858117/151‑1 (MQ=255)
                                                                                                                                         ggAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGt             <  2:197613/144‑1 (MQ=255)
                                                                                                                                         ggAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGc      >  1:29657/1‑151 (MQ=255)
                                                                                                                                              cTTACCCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATc  <  2:499756/150‑1 (MQ=255)
                                                                                                                                                 aaCCCCGTTTGGGCCCTTTTTACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATc  <  1:359226/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATC  >  NC_000913/3177384‑3177674

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