Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,596 C→T 6.1% P665P (CCG→CCA narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,5960CT 6.1% 209.6 / 2.3 65P665P (CCG→CCAnarZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base T (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGC  >  NC_000913/1540449‑1540744
                                                                                                                                                   |                                                                                                                                                    
ggagCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCgg                                                                                                                                                    >  1:940125/2‑150 (MQ=255)
cgagCGCCAGACAAACAACTTACGCGGATGATTGCGGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCgg                                                                                                                                                    <  2:1325973/149‑1 (MQ=255)
         gACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGccc                                                                                                                                           >  1:92937/1‑150 (MQ=255)
           caaacaaaTTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGATTTAATCCCGAGCGGGTTACTCCCCa                                                                                                                                         >  1:281523/1‑150 (MQ=255)
             aacaaaTTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGc                                                                                                                                       >  1:392288/1‑150 (MQ=255)
              acaaaTTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCGGGGCGGTAAATTCTGTGGGGGATAATCCGGCATTGTCGGCTTCAGCTTTAATCCCGAGCGCTTTACGCCCCAGCt                                                                                                                                      >  2:410297/1‑150 (MQ=255)
                 aaTTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTgcg                                                                                                                                   >  2:23016/1‑150 (MQ=255)
                     aCGCGGATGATTGCTGCTGCTATCTGGGTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGgcgc                                                                                                                               <  1:224470/150‑1 (MQ=255)
                                tgctgctgcTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTagcc                                                                                                                    <  1:416717/150‑1 (MQ=255)
                                  ctgctgctATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCCGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCCCAGAAGGTagccac                                                                                                                  >  1:632666/1‑149 (MQ=255)
                                    ccggctATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTGGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTagccagcc                                                                                                                <  2:339753/147‑1 (MQ=255)
                                      tgctATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCa                                                                                                              <  2:1311990/150‑1 (MQ=255)
                                        ctATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATa                                                                                                            >  2:397448/1‑150 (MQ=255)
                                         tATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATAc                                                                                                           >  1:261429/1‑150 (MQ=255)
                                            cTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGTGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGtt                                                                                                        >  2:973195/1‑150 (MQ=255)
                                            cTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGTGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGtt                                                                                                        >  2:978732/1‑150 (MQ=255)
                                             tGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTc                                                                                                       <  2:503731/150‑1 (MQ=255)
                                               gttgttCGCCGGCCATACGTAACTCGCCCGATTTCAGCGCCGGGGCGGTAAATGCTGTGGGGGCTAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCgg                                                                                                     <  1:851218/150‑1 (MQ=255)
                                               gttgttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGATCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACTTTCgt                                                                                                     >  1:63078/1‑149 (MQ=255)
                                               gttgttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCCAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGGATCCAGCCCATACGTTCgg                                                                                                     <  1:453039/150‑1 (MQ=255)
                                                     cGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTACGGCGCGaa                                                                                               >  1:738693/1‑150 (MQ=255)
                                                           ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGa                                                                                         <  2:780134/150‑1 (MQ=255)
                                                             aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                       >  2:368460/1‑150 (MQ=255)
                                                             aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatggg                                                                                       <  2:1044878/150‑6 (MQ=255)
                                                             aTACGTAAAGCGCCCGCTTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCGGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                       <  1:664606/150‑1 (MQ=255)
                                                               aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAAGCCGGCCTTGTCGGCATCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                     <  1:642001/150‑1 (MQ=255)
                                                                 gTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaa                                                                                   <  2:940125/150‑1 (MQ=255)
                                                                  tAAATCGCCCGATTTCAGCGCCTTGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCTGGTTACGCCCCAGCTGAGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaat                                                                                  <  1:371765/150‑1 (MQ=255)
                                                                   acaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAAACCGGCCTTGGCGGCTTCAGCTGTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                 <  1:254947/148‑1 (MQ=255)
                                                                   aaaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                 <  1:817601/150‑1 (MQ=255)
                                                                   aaaTCGCCCGATATCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                 <  1:867048/150‑1 (MQ=255)
                                                                     aTCGCCCGATTTCAGCGCCTGGGCTGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAg                                                                               <  1:719372/150‑1 (MQ=255)
                                                                     aTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAg                                                                               >  1:98765/1‑150 (MQ=255)
                                                                        gCCCGATTTCAGCGCCTGTGCGGTAAAATCTGTGGGGGATAATCCGGCCTTGTAGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                            <  1:420958/150‑1 (MQ=255)
                                                                        gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGGGGGGGATAATCCGGCCTTGTCGGCTTCAGGTTTAATCTCGAGCGGGTTTCGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                            <  1:498420/150‑1 (MQ=255)
                                                                         cccGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTAGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGGGAACGTTGAGATCATTCAGATGa                                                                           >  2:90756/1‑150 (MQ=255)
                                                                                 cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACTGCTGta                                                                   <  1:1309094/150‑1 (MQ=255)
                                                                                 cAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                   <  1:286354/150‑1 (MQ=255)
                                                                                         gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATTCGTTCGGCGCGAACGTTGAAATCTATCAGATGACCGCTGTACTTACTgg                                                           >  1:39403/1‑150 (MQ=255)
                                                                                                         gggggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGa                                           <  2:1194271/150‑1 (MQ=255)
                                                                                                          ggggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGCGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaa                                          >  1:57661/1‑150 (MQ=255)
                                                                                                           gggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaaa                                         >  2:672503/1‑150 (MQ=255)
                                                                                                           gggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaaa                                         <  2:92937/150‑1 (MQ=255)
                                                                                                            ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTATCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                        <  1:466506/150‑1 (MQ=255)
                                                                                                               tAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTaa                                     >  1:93423/1‑150 (MQ=255)
                                                                                                                   ccGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTcc                                 <  1:890379/150‑1 (MQ=255)
                                                                                                                     ggCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGCACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTg                               >  1:1108454/1‑150 (MQ=255)
                                                                                                                     ggCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTg                               >  1:969457/1‑150 (MQ=255)
                                                                                                                     ggCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTg                               <  2:392288/150‑1 (MQ=255)
                                                                                                                       ccTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAg                             >  1:411995/1‑150 (MQ=255)
                                                                                                                        cTTGGCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGc                            <  2:377/150‑1 (MQ=255)
                                                                                                                         ttGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCa                           <  2:93423/150‑1 (MQ=255)
                                                                                                                          tGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTAGAACTCAATCAGATGACCGCTGTACTTACTGGCATCCGCGAGCGGTGAAAGTAACTCCTGACCag                          >  1:309076/1‑150 (MQ=255)
                                                                                                                                  tCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACtttt                  >  2:609736/1‑150 (MQ=255)
                                                                                                                                    aGCTTTAATCCCGAGCGGGTTACGCCCCAGCTTCGGCGCAGAAGGTAGCCAGCACATACTTTCGGCGCGAACGTTGAAATCAATAAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAATTAACTCCTGAGCAGAGACttttta                >  1:546598/1‑149 (MQ=255)
                                                                                                                                      cTTTAATCCCGAGTGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCAt              <  1:978732/150‑1 (MQ=255)
                                                                                                                                      cTTTAATCCCGAGTGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGAACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCAt              <  1:973195/150‑1 (MQ=255)
                                                                                                                                      cTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCAt              >  1:1222822/1‑150 (MQ=255)
                                                                                                                                         tAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGGAGAGACTGTTTCATAgc           >  2:605919/1‑150 (MQ=255)
                                                                                                                                         tAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgc           >  1:450422/1‑150 (MQ=255)
                                                                                                                                              ccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCAt      >  1:761949/1‑150 (MQ=255)
                                                                                                                                               cGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCAtt     >  1:414628/1‑150 (MQ=255)
                                                                                                                                                gAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTACGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTg    <  2:738693/150‑1 (MQ=255)
                                                                                                                                                gAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGGTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCAGCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTATTCATAGCGCCATTg    <  1:714571/150‑1 (MQ=255)
                                                                                                                                                  gCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGc  <  1:23016/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGC  >  NC_000913/1540449‑1540744

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: