Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,096,826 A→C 16.3% N49T (AAC→ACC)  rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,096,8260AC16.3% 163.3 / 3.0 79N49T (AAC→ACC) rdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (5/10);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATC  >  NC_000913/3096677‑3096975
                                                                                                                                                     |                                                                                                                                                      
tCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGaaaa                                                                                                                                                        >  1:387195/1‑150 (MQ=255)
tCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGACTTGATATCGTGGCCCAAACCGCCCTCGGCGTTGCTTCCGCTGAAGTACCCGGCCTGACCGTTATCGaaaa                                                                                                                                                        >  1:1219133/1‑150 (MQ=255)
 ccTCTGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTGGATTCCGCTGAAGTAACCGGCCTGACCTTTATCGAAAAc                                                                                                                                                       <  1:808416/150‑1 (MQ=255)
     tGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGa                                                                                                                                                   <  1:616385/150‑1 (MQ=255)
      gCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGCGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAACCCGGCCTGACCTTTATCGAAACCGCGAt                                                                                                                                                  >  1:695084/1‑150 (MQ=255)
      gCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGTTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAt                                                                                                                                                  <  1:343692/150‑1 (MQ=255)
      gCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAt                                                                                                                                                  >  2:138820/1‑150 (MQ=255)
       cAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGTGACTTCGGTCTTGATATCGTGGCCCAAACAGACATCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAtt                                                                                                                                                 >  2:1118370/1‑150 (MQ=255)
         aaaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTcg                                                                                                                                               >  2:695572/1‑149 (MQ=255)
           aaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGCGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAACCCGGCCTGCCCTTTATCGAAAACGCGATTCTGa                                                                                                                                             >  1:521887/1‑150 (MQ=255)
           aaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGa                                                                                                                                             <  1:1280991/150‑1 (MQ=255)
            aGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGaa                                                                                                                                            >  2:374808/1‑150 (MQ=255)
              ttGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAg                                                                                                                                          <  1:908949/150‑1 (MQ=255)
                gTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAgcg                                                                                                                                        <  2:633998/150‑1 (MQ=255)
                  ccTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCGTTATCGAAAACGCGATTCTGAAAgcgcg                                                                                                                                      >  1:1105753/1‑150 (MQ=255)
                     cGCAACCGGCAATGTCGGTAAAGTGAGTGAGCTGGCGTCGCTGCTTAGCGACGTCGGTCTTGATATCGTGGCCCAAACAGACATCCGCGTTGATTCCGCTGAAGAACCCGCCCTGCCCTTTATCTAAACCGCGATTCTGAAAGCGCGCCa                                                                                                                                   >  2:1086764/1‑150 (MQ=255)
                         aCCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCg                                                                                                                               >  1:623526/1‑150 (MQ=255)
                          ccGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCgg                                                                                                                              <  2:375595/150‑1 (MQ=255)
                           ccgCAATGTCGGTAAAGTGCGGGACCTGGCGGCGCTGCTTAGCGAGGTCGGTCATGATATCGTGGCCCAAACAGACCTCGGGGTTGATTCCGATGAAGAAACCGGCCTTACCTTTATCTAAAACGCGATTCTGAAAGCGCGCCATGCGGc                                                                                                                             <  2:210501/148‑1 (MQ=255)
                           cGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGc                                                                                                                             <  1:711401/150‑1 (MQ=255)
                           cGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGc                                                                                                                             <  2:257272/150‑1 (MQ=255)
                           cGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGc                                                                                                                             <  2:626100/150‑1 (MQ=255)
                            ggCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCa                                                                                                                            >  1:1183746/1‑150 (MQ=255)
                            ggCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCa                                                                                                                            >  1:346984/1‑150 (MQ=255)
                                  gTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑g                                                                                                                     >  1:416676/1‑150 (MQ=255)
                                       tAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGc                                                                                                                <  1:890141/150‑1 (MQ=255)
                                        aaaGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGCCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGCAGAAACCGGCCTGACCTGTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCt                                                                                                               >  1:57639/1‑150 (MQ=255)
                                                tGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATGTCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGc                                                                                                       <  1:169375/150‑1 (MQ=255)
                                                   gCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAGTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAAt                                                                                                    >  1:589340/1‑150 (MQ=255)
                                                    cTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAAtt                                                                                                   >  2:367401/1‑150 (MQ=255)
                                                      ggCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGc                                                                                                 >  1:572081/1‑150 (MQ=255)
                                                       gCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGcc                                                                                                <  2:1342716/150‑1 (MQ=255)
                                                        cGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGccc                                                                                               >  1:569808/1‑149 (MQ=255)
                                                          gcGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCcgac                                                                                             <  1:1273270/149‑1 (MQ=255)
                                                            gctgctTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCcgacga                                                                                           >  1:349705/1‑150 (MQ=255)
                                                             ctgctTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCgacgac                                                                                          >  1:275817/1‑150 (MQ=255)
                                                                  tAGCGACTTCGGTCTTGATATCGTGCCCCAACCAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACtctgg                                                                                     <  2:939840/150‑1 (MQ=255)
                                                                  tAGCGACTTCGGTCGTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGATGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGATATCTCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGACGACGACtctgg                                                                                     >  2:1027406/1‑150 (MQ=255)
                                                                    gCGACTTCGGTCTTGATATCGGGGCCCAAACTGACCTCGGCGTTTATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCAATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACtctggtc                                                                                   <  2:689821/150‑1 (MQ=255)
                                                                        ctgcGGGCTGGATCTCGGGGCCCAAACCGCCCGCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGc                                                                               <  2:695084/147‑1 (MQ=255)
                                                                        cTTCGGTCTTGATATCGTGGCCCAAACACACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTCACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATCCGGCAAAAGT‑GACAGCTTAACCGCCTATTGCCGACGCCTCTGGTCTGGc                                                                               >  2:1213044/1‑150 (MQ=255)
                                                                               cTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAt                                                                        >  1:995208/1‑150 (MQ=255)
                                                                                tggATAACGTGGCCCAAACAGACCTCGGCGTTGATTCGGCTGATGAATCCGGCGTGACCTTTATCGATACGGCGATTATGAAAGCGCGCCATGCGGCAAAAGT‑GACGGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtg                                                                       <  2:94166/148‑1 (MQ=255)
                                                                                  tatatCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtgtg                                                                     <  2:698202/149‑1 (MQ=255)
                                                                                    tatCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCt                                                                   <  2:349705/150‑1 (MQ=255)
                                                                                    tatCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTATACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCt                                                                   >  1:497538/1‑150 (MQ=255)
                                                                                    aatCGTGCCCCAAACAGACCTCGGGGTTGATACCCCTGAATAAACGGGCCTGACCTTTATCGAAACCGCTATTCTGAACGCGCGCCATGCGGAAAATGT‑CACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCt                                                                   <  2:439505/149‑1 (MQ=255)
                                                                                     cgggTGGCCCAATCAGCCCCGGGCGTTGATCCCGCTGCACAACCCGCCCTGCCCTTTATCCAAACCGCCATTCTGACAGCGCGCCATGCGGCAAACGT‑GCCCTCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCtt                                                                  <  2:466968/147‑1 (MQ=255)
                                                                                     atCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGGAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCtt                                                                  >  1:548271/1‑150 (MQ=255)
                                                                                      tCGGGGCCCAACCCGACCTCGGCGTTGATTCCGCTGACTAAACCGGCCTGACCTTTATCGAAACCGCGATTCTGAACGCGCGCCATGCGGCAAAAGT‑GCCCGCTTTACCGGAATTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTg                                                                 <  2:623526/150‑1 (MQ=255)
                                                                                        gTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAATAAACCGGCCTGACCTTTATCTAAAACGCGATTCTTAAAGCGCGCCATGCGGAAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggc                                                               <  2:462411/150‑1 (MQ=255)
                                                                                          ggCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggcgg                                                             <  2:890043/150‑1 (MQ=255)
                                                                                            cccAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGgcggcg                                                           >  1:485098/1‑150 (MQ=255)
                                                                                            cccAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCGGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGgcggcg                                                           <  2:920419/150‑1 (MQ=255)
                                                                                              cAAACAGACCTCGGCGTTGATCCCGCTGACCAACCCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCATAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGgcgcg                                                         <  1:956098/150‑1 (MQ=255)
                                                                                                   cgcccgcggcggTGAGCCCGCGGACGAACCCGCCCGGACCTTTATCCCAACCGCGATTCTGACCGCGCGCCCGGCGGCAAACGT‑GCCCGCTTTACCGGCAATGGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCggg                                                    <  1:309653/144‑1 (MQ=255)
                                                                                                      ccTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAACCGCGATTCTGAAAGCGCGCCATGCGGCTAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGAAGATGTGCTTGGCGGCGCGCCGGGctt                                                 >  1:250874/1‑148 (MQ=255)
                                                                                                       cTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGCCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGAGCCGGGCAtt                                                >  1:717623/1‑150 (MQ=255)
                                                                                                        gcgccGTTGATCCCGCTGAAGAACCCGGCCTGACCTTTATCTAAACCGCGATTCTCAACGCGCGCCAGGCGGCAAAAGG‑CACCGCTTAACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAttt                                               <  1:306600/146‑1 (MQ=255)
                                                                                                        gcGGCGTTTATTCCGCTGAAGAACCCGGCATGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGAAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAttt                                               <  2:346984/149‑1 (MQ=255)
                                                                                                        ccGGCGTTGATCCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGTTTCTGACAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTCCCGGCAATTGCCGACGACTCTGGGCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAttt                                               <  1:563798/149‑1 (MQ=255)
                                                                                                          ggCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGCCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTAc                                             >  2:458540/1‑150 (MQ=255)
                                                                                                          ggCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGAAGATGTGCTTGGCGGCGCGCCGGGCATTTAc                                             >  2:777000/1‑150 (MQ=255)
                                                                                                             gTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACAGCTTTACCGCAAATTGCCGAAGACACTGTACTGGCGGTACATTTGCTTGGCGGCGCGCCGGGCAtttcaacc                                          >  2:1027192/1‑145 (MQ=255)
                                                                                                              ttGATTCCGCTTAAGAAACCGGCCGGACGTTTATCTAAAACGCTATTCTGAAAGCGCGCCTTGCGGAAAAAGG‑GCCGGCTTTACCGGAATTTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcg                                         <  2:1070968/150‑1 (MQ=255)
                                                                                                               tcaTTCCGCTGAACAAACCGGCCTTACCTTTATCTAAACCGCTATTCTGAACGCGCGCCATGCGGCAAAAGT‑CCCCGCTTTACCGCCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgc                                        <  2:275817/148‑1 (MQ=255)
                                                                                                                      ccttaacaaaCCGGCCTGACCTAAATCCAAAACGCGATTCTCAACGCGCGCCAGGCGGCAAAATT‑GACCGCTTTACCGGAATTTCCCGACGACTCTGGTCTGGCGGGAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGGTAt                                 <  1:1324659/143‑1 (MQ=255)
                                                                                                                         gaagagACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GCCCGCTTTACCGGCAATTGCCGACGCCTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                              >  1:1274882/1‑150 (MQ=255)
                                                                                                                         gaagaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                              >  2:937654/1‑150 (MQ=255)
                                                                                                                           agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGGTTTACGGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                            <  1:695572/150‑1 (MQ=255)
                                                                                                                           agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                            >  2:944237/1‑150 (MQ=255)
                                                                                                                           agaaACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgg                            >  2:750610/1‑150 (MQ=255)
                                                                                                                             aaaCCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGCCGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTg                          >  2:1234188/1‑150 (MQ=255)
                                                                                                                               aCCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGaa                        >  1:1230093/1‑150 (MQ=255)
                                                                                                                                     cTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                  >  2:961165/1‑150 (MQ=255)
                                                                                                                                       gacaattATGGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAATGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTTGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                <  2:56416/145‑1 (MQ=255)
                                                                                                                                       gACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                <  1:9095/150‑1 (MQ=255)
                                                                                                                                        acatttATCAAAACCGCTATTCTCAACCCGCGCCATCCGCAAAACGT‑CACGGCTTAACCGTAATTTCCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAcc               <  2:995208/147‑1 (MQ=255)
                                                                                                                                            ttATCTAAACCGCTATTCTGAACGCGCGCCATGCGGCAAAAGT‑CCCCGCTTAACCGCCAATGGCCCACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATc           <  1:466943/150‑1 (MQ=255)
                                                                                                                                            taccccaaacCGCTATTCTCAAAGC‑CCCCAAGCCGGAAAACTGGACCGATTAACCGTAATTGGCCGACGACTCTGGTCTGGCGGTAGATGTGGTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTAATCCGGTGAAGACGCGACCGATc           <  2:613810/144‑1 (MQ=255)
                                                                                                                                                cGAAAACGCGATTCGGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCaaaa       >  2:1027459/1‑150 (MQ=255)
                                                                                                                                                 taaaaCGCGATTCTGCAAGCGCGCCATGCGGCAAATGT‑GCCCGCTTTCCCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAg      <  2:99505/149‑1 (MQ=255)
                                                                                                                                                 gAAAACGCGATTCTGAAAGCGCGCCATGTGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACGCTGGTCTGGGGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTTTTCCGGTGATGACGCGACCGATCaaaac      >  2:49325/1‑149 (MQ=255)
                                                                                                                                                  aaaaCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCTACGACTCTGGTCTGGCGGAAGATTTGCTTGGCGGCGAGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAAGCTACCGATCaaaaca     >  2:1013061/1‑148 (MQ=255)
                                                                                                                                                   aaaCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGaa    >  1:724788/1‑150 (MQ=255)
                                                                                                                                                     aCGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATc  >  2:418476/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCGGCAAAAGT‑GACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATC  >  NC_000913/3096677‑3096975

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