Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,537,548 T→C 5.7% S436G (AGC→GGC)  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,5480TC 5.7% 293.8 / 4.9 88S436G (AGC→GGC) narYnitrate reductase 2 (NRZ), beta subunit
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base C (2/3);  total (39/49)
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.92e-01

CAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATTGCACC  >  NC_000913/1537403‑1537696
                                                                                                                                                 |                                                                                                                                                    
cAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAg                                                                                                                                                   <  2:54715/149‑1 (MQ=255)
      tCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTAc                                                                                                                                             <  2:244306/149‑1 (MQ=255)
        gCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaaga                                                                                                                                           <  1:224502/149‑8 (MQ=255)
         ccAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtc                                                                                                                                          >  1:502282/1‑149 (MQ=255)
           aaaGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtcgt                                                                                                                                        <  2:614215/149‑1 (MQ=255)
            aaGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCAAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtcgtc                                                                                                                                       <  3:242919/149‑1 (MQ=255)
              ggTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGa                                                                                                                                     <  2:83106/149‑1 (MQ=255)
               gTAAAACTGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAt                                                                                                                                    <  3:594143/149‑1 (MQ=255)
               gTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAt                                                                                                                                    <  4:28978/149‑1 (MQ=255)
                 aaaaCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCGTAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATgg                                                                                                                                  <  1:135261/149‑1 (MQ=255)
                  aaaCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGATCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGc                                                                                                                                 >  1:445591/1‑149 (MQ=255)
                   aaCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAt                                                                                                                                    >  3:573883/1‑145 (MQ=255)
                   aaCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAt                                                                                                                                    <  4:573883/145‑1 (MQ=255)
                    aCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaa                                                                                                                                             >  4:654547/1‑130 (MQ=255)
                    aCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaa                                                                                                                                             <  3:654547/135‑6 (MQ=255)
                      cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                             <  3:490293/149‑1 (MQ=255)
                      cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                             <  3:53187/149‑1 (MQ=255)
                       gCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtc                                                                                                                                                >  2:550790/1‑127 (MQ=255)
                       gCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtc                                                                                                                                                <  1:550790/129‑3 (MQ=255)
                         aGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTa                                                                                                                          >  4:265048/1‑149 (MQ=255)
                           ccGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATc                                                                                                                        <  4:221935/149‑1 (MQ=255)
                           ccGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATc                                                                                                                        <  1:53549/149‑1 (MQ=255)
                                gCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTaa                                                                                                                   <  2:315002/149‑1 (MQ=255)
                                       gCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTc                                                                                                            <  4:339191/149‑1 (MQ=255)
                                        cGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTcc                                                                                                           <  3:200177/149‑1 (MQ=255)
                                               aTCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGttt                                                                                                    >  2:5134/1‑149 (MQ=255)
                                                  gcccgccATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTtgtg                                                                                                 <  3:424130/149‑1 (MQ=255)
                                                     cgccATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAAc                                                                                              >  4:346212/1‑149 (MQ=255)
                                                      gccATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACg                                                                                             >  1:330094/1‑149 (MQ=255)
                                                        cATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCa                                                                                           <  2:288976/149‑1 (MQ=255)
                                                         aTTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCat                                                                                          <  1:292249/149‑1 (MQ=255)
                                                         aTTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCat                                                                                          <  3:425925/149‑1 (MQ=255)
                                                          tttCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCata                                                                                         <  2:499372/149‑1 (MQ=255)
                                                          tttCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCata                                                                                         <  2:330094/149‑1 (MQ=255)
                                                              ccGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGt                                                                                     <  3:119285/149‑1 (MQ=255)
                                                              ccGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGt                                                                                     >  3:297041/1‑149 (MQ=255)
                                                               cGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTg                                                                                    <  1:458909/149‑1 (MQ=255)
                                                                         tCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGcc                                                                          <  2:505179/149‑1 (MQ=255)
                                                                         tCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGca                                                                          <  2:143811/149‑2 (MQ=255)
                                                                          cGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCca                                                                         <  2:47932/149‑1 (MQ=255)
                                                                             gATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTTTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatca                                                                      >  4:464444/1‑149 (MQ=255)
                                                                             gATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtcgt                                                                                                                                        >  2:417729/1‑83 (MQ=255)
                                                                             gATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtcgt                                                                                                                                        <  1:417729/83‑1 (MQ=255)
                                                                             gATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatca                                                                      >  3:546109/1‑149 (MQ=255)
                                                                                aaCAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATcc                                                                   >  3:632919/1‑149 (MQ=255)
                                                                                aaCAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATcc                                                                   >  1:156916/1‑149 (MQ=255)
                                                                                    aaaCGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCGGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGttt                                                               >  2:591756/1‑149 (MQ=255)
                                                                                     aaCGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATc                                                                             >  3:295927/1‑134 (MQ=255)
                                                                                     aaCGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATc                                                                             <  4:295927/134‑1 (MQ=255)
                                                                                      aCGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTGCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCa                                                             <  4:496941/149‑1 (MQ=255)
                                                                                         aTCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAgcg                                                          >  4:65952/1‑149 (MQ=255)
                                                                                         aTCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAgcg                                                          >  2:406632/1‑149 (MQ=255)
                                                                                          tCTTCATAGTTGGCAATGGCGAGAAAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAGGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGAGGCGCATCGCCATCATCCGTTTCAgcgc                                                         >  3:433425/1‑149 (MQ=255)
                                                                                            ttCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGccc                                                       >  3:237453/1‑149 (MQ=255)
                                                                                            ttCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGccc                                                       >  2:562895/1‑149 (MQ=255)
                                                                                             tCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCa                                                                                                          <  1:4990/97‑1 (MQ=255)
                                                                                             tCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCa                                                                                                          >  2:4990/1‑97 (MQ=255)
                                                                                                   ttGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTac                                                >  3:642396/1‑149 (MQ=255)
                                                                                                   ttGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCATACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGAGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTac                                                <  4:376966/149‑1 (MQ=255)
                                                                                                       cAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccgga                                            <  3:261106/149‑1 (MQ=255)
                                                                                                       cAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccgga                                            <  3:675798/149‑1 (MQ=255)
                                                                                                          tGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggaccg                                         <  3:511605/149‑1 (MQ=255)
                                                                                                            gCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGt                                       <  4:618807/149‑1 (MQ=255)
                                                                                                             cGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTa                                      >  3:33775/1‑149 (MQ=255)
                                                                                                                   aGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCg                                <  3:346212/149‑1 (MQ=255)
                                                                                                                         aCATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCTCCCGCAGTACCGGACCGGTATCGCCGGCACTc                          >  3:497112/1‑149 (MQ=255)
                                                                                                                            tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAAc                       <  2:420696/149‑1 (MQ=255)
                                                                                                                            tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAAc                       >  2:39906/1‑149 (MQ=255)
                                                                                                                              tctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTAt                                     <  4:284269/133‑1 (MQ=255)
                                                                                                                              tctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTAt                                     >  3:284269/1‑133 (MQ=255)
                                                                                                                              tctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACat                     >  4:553754/1‑149 (MQ=255)
                                                                                                                                ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatc                                                                       <  4:27658/97‑1 (MQ=255)
                                                                                                                                ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatc                                                                       >  3:27658/1‑97 (MQ=255)
                                                                                                                                ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACatat                   >  2:557404/1‑149 (MQ=255)
                                                                                                                                ttCGACCGGGGCGACGCTCAGCCCTCCTTCGACGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCAGATAGTGGCGCATACCCATCATACGTTTCAGCGTCCGCAATACCGGAGCGGTATCGCCCGCAGTCAACaaag                   >  2:197099/1‑146 (MQ=255)
                                                                                                                                 tCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATAtt                  <  4:237453/149‑1 (MQ=255)
                                                                                                                                 tCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATAtt                  <  2:156916/149‑1 (MQ=255)
                                                                                                                                  cGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTg                 >  1:288040/1‑149 (MQ=255)
                                                                                                                                    aCCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGc               >  2:416935/1‑149 (MQ=255)
                                                                                                                                     ccTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCg              >  3:688746/1‑149 (MQ=255)
                                                                                                                                      cTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGTGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCga             <  3:265048/149‑1 (MQ=255)
                                                                                                                                      cTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGc               <  2:398393/147‑1 (MQ=255)
                                                                                                                                      cTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGc               >  1:398393/1‑147 (MQ=255)
                                                                                                                                      cTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCga             <  1:557404/149‑1 (MQ=255)
                                                                                                                                          gCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGata         >  4:521468/1‑149 (MQ=255)
                                                                                                                                            gACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACat                     <  2:361505/135‑1 (MQ=255)
                                                                                                                                            gACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACat                     >  1:361505/1‑135 (MQ=255)
                                                                                                                                            gACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATAtt       >  3:83108/1‑149 (MQ=255)
                                                                                                                                                 tCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATTGCAcc  >  3:200/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATTGCACC  >  NC_000913/1537403‑1537696

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