Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,537,972 G→T 7.7% R294R (CGC→CGA narY ← nitrate reductase 2 (NRZ), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,537,9720GT 7.7% 310.6 / 3.5 120R294R (CGC→CGAnarYnitrate reductase 2 (NRZ), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (45/63);  new base T (10/2);  total (55/65)
Fisher's exact test for biased strand distribution p-value = 1.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GGTGCAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAA  >  NC_000913/1537824‑1538121
                                                                                                                                                    |                                                                                                                                                     
ggTGCAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGc                                                                                                                                                      <  1:55098/150‑1 (MQ=255)
    cAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGTAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTc                                                                                                                                                  <  2:203201/150‑1 (MQ=255)
    cAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTc                                                                                                                                                  <  2:1272599/150‑1 (MQ=255)
    cAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTc                                                                                                                                                  <  2:496451/150‑1 (MQ=255)
    cAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCATCGATCACTGAGTGATAGTGTGGATCGAGGAACACTTCGCACTGGCGtta                                                                                                                                                  >  1:1115966/1‑149 (MQ=255)
      aCGGTAGTGCCAGTTTCCAGTCCATCGCGATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCAt                                                                                                                                                >  1:1902029/1‑150 (MQ=255)
         gTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATaga                                                                                                                                             <  1:437094/150‑1 (MQ=255)
               ccAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCaa                                                                                                                                       >  1:1600146/1‑150 (MQ=255)
               ccAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTAAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAAAACTTCGCACTGGCGGTCATAGAGGTCaa                                                                                                                                       >  1:1201344/1‑150 (MQ=255)
                 aGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACt                                                                                                                                     <  2:599411/150‑1 (MQ=255)
                     tCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTcgc                                                                                                                                 <  1:880520/150‑1 (MQ=255)
                     tCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTcgc                                                                                                                                 >  2:693403/1‑150 (MQ=255)
                      ccAGTCCATCGCCATTATGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATATTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTcgcg                                                                                                                                <  2:1459256/150‑1 (MQ=255)
                       cAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTAAATAGAGGTCAACTTCgcgc                                                                                                                               >  2:960599/1‑150 (MQ=255)
                        aGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCt                                                                                                                              >  1:1619790/1‑150 (MQ=255)
                          tCCATCGCCATTTTGTAGACTGGCGAACGCTGGTCAGCGTCAATCACATTTTGTGGAATACATTTTTTCAGGTCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCg                                                                                                                            <  2:1218495/150‑1 (MQ=255)
                           ccATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCgg                                                                                                                           >  2:661537/1‑150 (MQ=255)
                            cATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGTGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGt                                                                                                                          <  2:497674/150‑1 (MQ=255)
                            cATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGt                                                                                                                          <  2:1887565/150‑1 (MQ=255)
                                gCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTc                                                                                                                      >  1:436544/1‑150 (MQ=255)
                                      ttGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCt                                                                                                                <  1:1336628/150‑1 (MQ=255)
                                      ttGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTCTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGTACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCt                                                                                                                >  1:1426069/1‑150 (MQ=255)
                                       tGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGACCACTTCGCACTGGCGTTAATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCtt                                                                                                               >  2:1697550/1‑150 (MQ=255)
                                        gTATACTGGCGAACGCGGGGCAGCGTCAATCACATTTTGGGGAATACCTTGTTTCAGGGCGTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGATCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTc                                                                                                              <  1:140905/150‑1 (MQ=255)
                                           gACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTAGATCGAGGAACACTTCGCACTTGCGTTCAAAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTc                                                                                                           >  2:1814582/1‑150 (MQ=255)
                                             cTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGATTACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGa                                                                                                         >  1:1917414/1‑150 (MQ=255)
                                                gCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTAGCGCTCGGTGCTCGCCGCTTCCTCGATGc                                                                                                      >  1:336207/1‑150 (MQ=255)
                                                 cGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCg                                                                                                     <  1:660238/150‑1 (MQ=255)
                                                  gAACGCGGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCATCCTCGATCACTGAGGGATCGTGTGGATCGTGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCTGTGCTCGCCGCTTCCTCGATGCgg                                                                                                    <  1:773739/150‑1 (MQ=255)
                                                    aCGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCTAGGAACACTACGCAATGTATTTCATATAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTc                                                                                                  >  2:1978601/1‑150 (MQ=255)
                                                    aCGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTAGATCTAGGAACACTTCGCACAGGATTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTc                                                                                                  >  2:1923725/1‑150 (MQ=255)
                                                     cGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCg                                                                                                 <  2:1834576/150‑1 (MQ=255)
                                                       cTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCTTAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGc                                                                                               <  1:1453565/150‑1 (MQ=255)
                                                         gggCAGCGTCAATCACTTTTTGTGGAATACCTTGTTTCAGGGCATCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGcgt                                                                                             <  1:1762103/150‑1 (MQ=255)
                                                           gCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGTGCTTCCTCGATCACTGAGTGATCGTGTGGATCGAGGAACACTTCGCACTGGCTTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGcgtcg                                                                                           <  2:1533712/150‑1 (MQ=255)
                                                           gCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGcgtcg                                                                                           <  2:549157/150‑1 (MQ=255)
                                                           gCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGcgtcg                                                                                           <  2:1686497/150‑1 (MQ=255)
                                                            cAGCGTCAATCTCATTTTGTGCACTACCTTGTTTCAGGGCTTCCTCGATCACTGAGGTATCGTGTGTATCGAGGAACACTTCGCACTGGCGTTAATAGACGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGcgtagt                                                                                          <  2:1319097/150‑1 (MQ=255)
                                                              gCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTaa                                                                                        >  1:1810300/1‑150 (MQ=255)
                                                                 tCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAagc                                                                                     <  1:1830219/150‑1 (MQ=255)
                                                                 tCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAagc                                                                                     <  1:568362/150‑1 (MQ=255)
                                                                   aaTCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTTGATCGAGGAACACTTCGAAATGTATTTAATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAagcag                                                                                   >  2:1359396/1‑150 (MQ=255)
                                                                    aTCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAagcagc                                                                                  <  1:1153658/150‑1 (MQ=255)
                                                                    aTCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAagcagc                                                                                  <  1:1773693/150‑1 (MQ=255)
                                                                        caTTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTAGATCTAGGAACACTACGAAAAATAAAAAAAAAATGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGc                                                                              >  2:479765/1‑150 (MQ=255)
                                                                         aTTTTGTGGAAAACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGGCGTAAAGCAGCACGcc                                                                             <  1:1076342/150‑1 (MQ=255)
                                                                           tttGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCa                                                                           <  2:6315/150‑1 (MQ=255)
                                                                              gtgGAATACCTTGTTTCAGGGCTTCCTCGCTCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGAt                                                                        <  1:1012170/150‑1 (MQ=255)
                                                                              gtgGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTTGGCGTCGTAAAGCAGCACGCCCAGAt                                                                        <  1:1602754/150‑1 (MQ=255)
                                                                              gtgGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGAt                                                                        >  1:424442/1‑150 (MQ=255)
                                                                              gtgGAATACCTTGTTTCAGGGCTTCCTCGAGCTCTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCGCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGAt                                                                        <  1:1005313/150‑1 (MQ=255)
                                                                               tgGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATa                                                                       >  1:852321/1‑150 (MQ=255)
                                                                                ggacttccTTGGTTCAGGGCTGCCTCGATCACTGAGGGTTCGTGGGGAGCGAGGATCACTTCGCACTGGCGTGCATAGAGGTCAACTGAGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTATAGCAGCACGCCCAGATAc                                                                      <  1:85483/144‑1 (MQ=255)
                                                                                ggAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGATCACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCCATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATAc                                                                      >  1:262168/1‑150 (MQ=255)
                                                                                  aaTACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCAATGGAATTCATATAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCg                                                                    >  2:1573107/1‑150 (MQ=255)
                                                                                   aTACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCgg                                                                   <  2:434934/150‑1 (MQ=255)
                                                                                      ccTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTTGATCTAGGAACACTACGAAAACTAATTAATATAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAACGCCGCACGCGCAGATACCGGATg                                                                >  2:702798/1‑150 (MQ=255)
                                                                                      ccTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGGCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATg                                                                <  2:614348/150‑1 (MQ=255)
                                                                                      ccTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATg                                                                >  1:1323790/1‑150 (MQ=255)
                                                                                      ccTTGTTTCAGGGCTTCCTCGAGCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATg                                                                <  2:1842124/150‑1 (MQ=255)
                                                                                       cTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGc                                                               <  2:399976/150‑1 (MQ=255)
                                                                                        ttGTTTCAGGGCTTCCTCCATCACAGAGGGATCGTGTTGATCGAGTAACACTTCGCACTGTCGTTCATAGAGGTCAACTTCGCGCTCTGTGCTCGCCGATTCTTTGATGAGGTCTGCGTCGTATAGCAGCACGCCCAGATACCGGATGCg                                                              <  2:132877/150‑1 (MQ=255)
                                                                                         tGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGa                                                             <  1:639045/150‑1 (MQ=255)
                                                                                          gTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTTGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCTCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGAc                                                            <  1:1190442/150‑1 (MQ=255)
                                                                                             tCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCa                                                         >  2:265465/1‑150 (MQ=255)
                                                                                                gggCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGc                                                      >  1:1534585/1‑150 (MQ=255)
                                                                                                  gCTTCCTCGATCACTGAGGGATCGTGTTGATCGAGGAACACTACGAAATATAAATCATATAGTTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAg                                                    >  2:1998644/1‑150 (MQ=255)
                                                                                                  gCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAg                                                    >  2:1529232/1‑150 (MQ=255)
                                                                                                  gCTTCCTCGATCACTGAGGGATCGTGTAGATCTAGGAACACTACGAACTATAGATCATAGAGGTCAAATTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCTGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGAGACCCACGCAg                                                    >  2:1529235/1‑150 (MQ=255)
                                                                                                   cTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAgg                                                   <  1:273292/150‑1 (MQ=255)
                                                                                                       gtCGATCACTGAGGGATCGTGTGGATCTAGGAACACTTCGCACGGTCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTc                                               <  1:381570/149‑1 (MQ=255)
                                                                                                       cTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAATGCAGCACGCCCAGATACCGGTAGCGACCCACGCAGGTTTc                                               >  1:198594/1‑150 (MQ=255)
                                                                                                       cTCGATCACTGAGGGATCGTGTAGATCTAGGAACACTACGCAAAATATTAAATATATGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTc                                               >  2:907120/1‑150 (MQ=255)
                                                                                                         cGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTg                                             <  2:768290/150‑1 (MQ=255)
                                                                                                         cGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTg                                             <  2:436544/150‑1 (MQ=255)
                                                                                                          gATCACTGAGGGATCGTGTAGATCTAGGAACACTACGAAAAATAAAAAAAAAAAATCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTATCTGa                                            >  2:807495/1‑150 (MQ=255)
                                                                                                                 gAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACAcg                                     <  1:825366/150‑1 (MQ=255)
                                                                                                                   gggATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACAcggt                                   >  1:1614374/1‑150 (MQ=255)
                                                                                                                         gtgtGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTg                             <  1:1600313/150‑1 (MQ=255)
                                                                                                                          tgtgGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGa                            <  1:265465/150‑1 (MQ=255)
                                                                                                                          tgtgGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGa                            >  2:519886/1‑150 (MQ=255)
                                                                                                                          tgtgGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGACCACACGGACGGTTgc                            >  2:197165/1‑149 (MQ=255)
                                                                                                                          tgtgGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGACGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGCCCCACGCAGGTTTCTGAGCACACGGGCGGTAGa                            >  2:91375/1‑150 (MQ=255)
                                                                                                                            tgGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGAcc                          <  2:1425467/150‑1 (MQ=255)
                                                                                                                            cggATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGAcc                          <  1:1682251/149‑1 (MQ=255)
                                                                                                                             ggATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCATGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCg                         >  2:670675/1‑150 (MQ=255)
                                                                                                                             ggATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCg                         <  1:693403/150‑1 (MQ=255)
                                                                                                                               aTCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGCGCACACGGTCGGTTGACCGGa                       >  1:1303/1‑150 (MQ=255)
                                                                                                                               aTCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCggc                       >  1:469243/1‑149 (MQ=255)
                                                                                                                                tCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGCCCGGAc                      >  2:112757/1‑150 (MQ=255)
                                                                                                                                tCGAGGAACACTTCGCACTGGCGTTCATAGAGGGCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGAc                      <  1:469705/150‑1 (MQ=255)
                                                                                                                                 cGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGTGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACt                     <  2:534092/150‑1 (MQ=255)
                                                                                                                                  gAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTc                    <  2:167505/150‑1 (MQ=255)
                                                                                                                                          cTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTg            >  2:1161973/1‑150 (MQ=255)
                                                                                                                                          cTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCAGACTCAATTCGTg            <  2:2007036/150‑1 (MQ=255)
                                                                                                                                          cTTCGCACGGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGAGGCGGTCGGCGTCGTCAAGCAGCACGCCCAGATACGGGATGCGCCGCCCGCTGGTTTCGGAGCACACGGTCGCTTGACCGGGCTCAATTCGTg            >  1:283098/1‑150 (MQ=255)
                                                                                                                                            tctcACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTggg          <  2:1426069/147‑1 (MQ=255)
                                                                                                                                              gCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCGTCCTCGATGCGGTCGGCGTCGTCCAGCAGCTCGCCCAGATACCGGATGCGCCCCACGCTGGTTTCTGCGCCCCCGGTCGGTGGCCCGGTCGCAGTACGTgcttt        >  1:176390/1‑146 (MQ=255)
                                                                                                                                                   ggCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACaga   >  1:641932/1‑150 (MQ=255)
                                                                                                                                                    gCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACagaa  <  2:1810300/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTGCAACGGTAGTGCCAGTTTCCAGTCCATCGCCATTTTGTAGACTGGCGAACGCTGGGCAGCGTCAATCACATTTTGTGGAATACCTTGTTTCAGGGCTTCCTCGATCACTGAGGGATCGTGTGGATCGAGGAACACTTCGCACTGGCGTTCATAGAGGTCAACTTCGCGCTCGGTGCTCGCCGCTTCCTCGATGCGGTCGGCGTCGTAAAGCAGCACGCCCAGATACCGGATGCGACCCACGCAGGTTTCTGAGCACACGGTCGGTTGACCGGACTCAATTCGTGGGTAACAGAA  >  NC_000913/1537824‑1538121

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: