Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,537,540 G→A 5.7% A438A (GCC→GCT 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,5400GA 5.7% 286.2 / 5.6 87A438A (GCC→GCTnarYnitrate reductase 2 (NRZ), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (35/47);  new base A (2/3);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

CCGAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGA  >  NC_000913/1537392‑1537685
                                                                                                                                                    |                                                                                                                                                 
ccGAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTgg                                                                                                                                                   <  3:567903/149‑1 (MQ=255)
ccGAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTgg                                                                                                                                                   <  2:665911/149‑1 (MQ=255)
  gAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATACCAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGc                                                                                                                                                 >  4:594143/1‑149 (MQ=255)
  gAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGACCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGc                                                                                                                                                 <  4:475585/149‑1 (MQ=255)
    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                                                                                                                         <  4:573883/145‑1 (MQ=255)
                              aaCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGAt                                                                                                                         >  3:573883/1‑145 (MQ=255)
                               aCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaa                                                                                                                                  <  3:654547/135‑6 (MQ=255)
                               aCCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGAGCGTCGCCCAggtcgaa                                                                                                                                  >  4:654547/1‑130 (MQ=255)
                                 cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                  <  3:53187/149‑1 (MQ=255)
                                 cGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGa                                                                                                                  <  3:490293/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                                                                                                             <  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                                                                               <  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                                                                                                                             <  1:417729/83‑1 (MQ=255)
                                                                                        gATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTtcgt                                                                                                                             >  2:417729/1‑83 (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                                               >  2:406632/1‑149 (MQ=255)
                                                                                                    aTCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAgcg                                               >  4:65952/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:39906/1‑149 (MQ=255)
                                                                                                                                       tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAAc            <  2:420696/149‑1 (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       <  2:156916/149‑1 (MQ=255)
                                                                                                                                            tCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATAtt       <  4:237453/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)
                                                                                                                                                    |                                                                                                                                                 
CCGAACCGTGGCAACCGTCGCCAAAGGTAAAACCGCAGCCGTTGCGTTCTGCGAAGGCATCGCCCGCCATTTCCCGATGGCTGGTCGGGATAACAAAACGATCTTCATAGTTGGCAATGGCGAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGA  >  NC_000913/1537392‑1537685

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: