Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,537,544 A→T 5.7% V437D (GTC→GAC)  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,5440AT 5.7% 302.2 / 7.0 88V437D (GTC→GAC) narYnitrate reductase 2 (NRZ), beta subunit
Reads supporting (aligned to +/- strand):  ref base A (37/46);  new base T (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 = 3.42e-01

ACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATT  >  NC_000913/1537396‑1537691
                                                                                                                                                    |                                                                                                                                                   
aCCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGa                                                                                                                                                     <  3:655369/149‑1 (MQ=255)
       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                                                                                                                                        <  3:654547/135‑6 (MQ=255)
                           aCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaa                                                                                                                                        >  4:654547/1‑130 (MQ=255)
                             cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                        <  3:490293/149‑1 (MQ=255)
                             cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                        <  3:53187/149‑1 (MQ=255)
                              gCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtc                                                                                                                                           <  1:550790/129‑3 (MQ=255)
                              gCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtc                                                                                                                                           >  2:550790/1‑127 (MQ=255)
                                aGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTa                                                                                                                     >  4:265048/1‑149 (MQ=255)
                                  ccGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATc                                                                                                                   <  1:53549/149‑1 (MQ=255)
                                  ccGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATc                                                                                                                   <  4:221935/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                                                                                     <  3:425925/149‑1 (MQ=255)
                                                                aTTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCat                                                                                     <  1:292249/149‑1 (MQ=255)
                                                                 tttCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCata                                                                                    <  2:330094/149‑1 (MQ=255)
                                                                 tttCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCata                                                                                    <  2:499372/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                                                                        <  4:295927/134‑1 (MQ=255)
                                                                                            aaCGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATc                                                                        >  3:295927/1‑134 (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:675798/149‑1 (MQ=255)
                                                                                                              cAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccgga                                       <  3:261106/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                                >  3:284269/1‑133 (MQ=255)
                                                                                                                                     tctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTAt                                <  4:284269/133‑1 (MQ=255)
                                                                                                                                     tctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACat                >  4:553754/1‑149 (MQ=255)
                                                                                                                                       ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatc                                                                  >  3:27658/1‑97 (MQ=255)
                                                                                                                                       ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCcatc                                                                  <  4:27658/97‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
ACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATT  >  NC_000913/1537396‑1537691

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: