Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,264,246 G→A intergenic (‑15/‑68) qorA ← / → dnaB quinone oxidoreductase, NADPH‑dependent/replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,264,2460GA99.3% 571.7 / NA 145intergenic (‑15/‑68)qorA/dnaBquinone oxidoreductase, NADPH‑dependent/replicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (74/69);  minor base T (0/1);  total (75/70)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCC  >  NC_000913/4264098‑4264393
                                                                                                                                                    |                                                                                                                                                   
gtgtCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATa                                                                                                                                                   <  4:1308362/151‑1 (MQ=255)
gtgtCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATa                                                                                                                                                   <  2:2444/151‑1 (MQ=255)
      aTAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTc                                                                                                                                              >  1:843952/1‑150 (MQ=255)
        aaaaTTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAg                                                                                                                                            <  1:292826/150‑1 (MQ=255)
        aaaaTTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGGGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTc                                                                                                                                              >  1:895507/1‑148 (MQ=255)
         aaaTTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAggg                                                                                                                                          <  1:265884/151‑1 (MQ=255)
          aaTTGATGCCGATGGCTTTATTTTCGCCCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGt                                                                                                                                         <  4:337762/151‑1 (MQ=255)
           aTTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTc                                                                                                                                        >  1:1592729/1‑151 (MQ=255)
             tGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTg                                                                                                                                      >  2:257499/1‑151 (MQ=255)
                tGCCGATGGCTTTATTTTCGTCCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAATTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCt                                                                                                                                       <  3:171522/147‑1 (MQ=255)
                tGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCtt                                                                                                                                   >  1:61329/1‑151 (MQ=255)
                  ccGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTc                                                                                                                                  >  4:663308/1‑150 (MQ=255)
                  ccGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCa                                                                                                                                 <  4:1102355/151‑1 (MQ=255)
                    gATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCata                                                                                                                               <  1:967763/151‑1 (MQ=255)
                     aTGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCatat                                                                                                                              <  4:696085/151‑1 (MQ=255)
                       ggCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATg                                                                                                                             <  1:1584330/150‑1 (MQ=255)
                       ggCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATg                                                                                                                             >  4:716704/1‑150 (MQ=255)
                        gCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGAt                                                                                                                           <  1:1664369/151‑1 (MQ=255)
                         cTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGAt                                                                                                                           >  4:1251125/1‑150 (MQ=255)
                          tttATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACCTATCGTCAGGGTCTGCTTCATATGAt                                                                                                                           >  3:299079/1‑149 (MQ=255)
                           ttATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATaa                                                                                                                         <  1:1548592/150‑1 (MQ=255)
                           ttATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATa                                                                                                                          <  4:1133400/149‑1 (MQ=255)
                            tATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATaaa                                                                                                                        <  1:239393/150‑1 (MQ=255)
                             aTTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAg                                                                                                                       >  4:263531/1‑150 (MQ=255)
                             aTTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGt                                                                                                                      >  1:810940/1‑151 (MQ=255)
                               tttCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGtt                                                                                                                     <  1:1549346/150‑1 (MQ=255)
                                 tCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCg                                                                                                                  >  4:1678111/1‑151 (MQ=255)
                                  cGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCg                                                                                                                  <  1:968592/150‑1 (MQ=255)
                                  cGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCg                                                                                                                  >  4:365339/1‑150 (MQ=255)
                                  cGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGa                                                                                                                 <  4:1732570/151‑1 (MQ=255)
                                      cTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGAccc                                                                                                              >  3:512376/1‑150 (MQ=255)
                                        ggATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCa                                                                                                             <  1:191762/149‑1 (MQ=255)
                                        ggATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCAtt                                                                                                           >  3:520004/1‑151 (MQ=255)
                                        ggATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCAtt                                                                                                           <  3:517782/151‑1 (MQ=255)
                                         gATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTc                                                                                                          <  1:466156/151‑1 (MQ=255)
                                          aTTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCt                                                                                                         >  3:385263/1‑151 (MQ=255)
                                          aTTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCt                                                                                                         >  4:152454/1‑151 (MQ=255)
                                          aTTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCt                                                                                                         >  2:223224/1‑151 (MQ=255)
                                            ttCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCttt                                                                                                       >  1:727189/1‑151 (MQ=255)
                                            ttCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCtt                                                                                                        <  1:1185680/150‑1 (MQ=255)
                                             tCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCttt                                                                                                       <  3:330114/150‑1 (MQ=255)
                                              cATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTAt                                                                                                     <  2:170044/151‑1 (MQ=255)
                                                 tctcCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTAtct                                                                                                   <  3:319598/150‑1 (MQ=255)
                                                      gccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTa                                                                                              >  4:1551822/1‑150 (MQ=255)
                                                      gccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTa                                                                                              >  1:1403817/1‑150 (MQ=255)
                                                      gccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGCGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTa                                                                                              >  3:281967/1‑150 (MQ=255)
                                                       ccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTaa                                                                                             <  1:636466/150‑1 (MQ=255)
                                                       ccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAAc                                                                                            >  1:1855292/1‑151 (MQ=255)
                                                             cGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCAt                                                                                       <  4:1104397/150‑1 (MQ=255)
                                                              ggCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTc                                                                                     <  1:1113797/151‑1 (MQ=255)
                                                               gCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCa                                                                                    <  1:1248770/151‑1 (MQ=255)
                                                                 aGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACt                                                                                  <  4:929388/151‑1 (MQ=255)
                                                                 aGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACt                                                                                  >  2:51888/1‑151 (MQ=255)
                                                                 aGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACt                                                                                  >  1:458368/1‑151 (MQ=255)
                                                                   gAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTAt                                                                                <  1:1894053/151‑1 (MQ=255)
                                                                     gTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATgg                                                                              <  1:1656510/151‑1 (MQ=255)
                                                                         aCTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAgg                                                                          >  1:665143/1‑151 (MQ=255)
                                                                            ctACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGa                                                                         <  1:71574/149‑1 (MQ=255)
                                                                              aCGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGaaa                                                                       <  3:66863/149‑1 (MQ=255)
                                                                              aCGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAAt                                                                      <  4:130085/150‑1 (MQ=255)
                                                                              aCGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATa                                                                     <  4:1681748/151‑1 (MQ=255)
                                                                               cGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATaa                                                                    <  1:1559076/151‑1 (MQ=255)
                                                                                 gCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAAc                                                                  <  4:1198016/151‑1 (MQ=255)
                                                                                    tGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAAccc                                                                <  4:1575816/150‑1 (MQ=255)
                                                                                     gAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCtt                                                              <  1:313862/151‑1 (MQ=255)
                                                                                          aCTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTcaac                                                          <  4:1732816/150‑1 (MQ=255)
                                                                                          aCTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTcaac                                                          >  2:503010/1‑150 (MQ=255)
                                                                                             tCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCaacaaa                                                       <  2:500823/150‑1 (MQ=255)
                                                                                               cGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAAcag                                                    <  4:1165768/151‑1 (MQ=255)
                                                                                                 ggCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAAcagca                                                  <  1:1386966/151‑1 (MQ=255)
                                                                                                  gCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAAcagcag                                                 <  3:314580/151‑1 (MQ=255)
                                                                                                  gCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAAcagc                                                   >  4:284654/1‑149 (MQ=255)
                                                                                                   ccACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAgg                                                <  1:1285601/151‑1 (MQ=255)
                                                                                                     aCCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGc                                               >  4:18442/1‑150 (MQ=255)
                                                                                                       cGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGa                                            <  4:387929/151‑1 (MQ=255)
                                                                                                            ttGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAAcccc                                       >  1:1011123/1‑151 (MQ=255)
                                                                                                             tgtgAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAAcccc                                       <  1:508445/150‑1 (MQ=255)
                                                                                                             tgtgAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAAcccc                                       <  4:907202/150‑1 (MQ=255)
                                                                                                               tgAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCcgcg                                    >  1:51530/1‑151 (MQ=255)
                                                                                                               tgAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCcgcg                                    >  1:1796630/1‑151 (MQ=255)
                                                                                                                 aaaTTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGaa                                  <  1:1094854/151‑1 (MQ=255)
                                                                                                                 aaaTTCAATTCGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGa                                   >  1:709205/1‑150 (MQ=255)
                                                                                                                   attcaattcGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAAcg                                >  1:1821908/1‑151 (MQ=255)
                                                                                                                      caattcGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAAcgcg                              >  1:562893/1‑150 (MQ=255)
                                                                                                                      caattcGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAAcgcg                              >  1:1899132/1‑150 (MQ=255)
                                                                                                                       aattcGTGTTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGAt                            >  4:1687304/1‑151 (MQ=255)
                                                                                                                             tgtTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCcaca                       >  1:1174304/1‑150 (MQ=255)
                                                                                                                             tgtTGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCcaca                       >  4:1782091/1‑150 (MQ=255)
                                                                                                                               ttGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGt                    >  4:1478718/1‑151 (MQ=255)
                                                                                                                                tGCCATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGtt                   >  1:326081/1‑151 (MQ=255)
                                                                                                                                  ccATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGc                 >  4:1277481/1‑151 (MQ=255)
                                                                                                                                   cATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGcc                >  1:728342/1‑151 (MQ=255)
                                                                                                                                   cATGTGTCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGcc                >  4:611173/1‑151 (MQ=255)
                                                                                                                                      gtgtCCTCCAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCgg              <  4:326187/150‑1 (MQ=255)
                                                                                                                                         tcctccAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTg          >  4:1395372/1‑151 (MQ=255)
                                                                                                                                           ctccAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGa         >  1:631283/1‑150 (MQ=255)
                                                                                                                                            tccAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGaa        >  4:779896/1‑150 (MQ=255)
                                                                                                                                            tccAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGaa        >  4:1576396/1‑150 (MQ=255)
                                                                                                                                              cAGAACATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAg      >  3:112975/1‑150 (MQ=255)
                                                                                                                                                 aaCATATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGcc  >  3:68908/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCC  >  NC_000913/4264098‑4264393

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