Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,584 A→G 6.4% V57A (GTC→GCC) ‡ napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,302,5840AG 6.4% 260.7 / 6.9 96V57A (GTC→GCC) ‡napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base G (4/2);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

GTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAG  >  NC_000913/2302436‑2302733
                                                                                                                                                        |                                                                                                                                                     
gTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc                                                                                                                                                      >  2:979557/1‑150 (MQ=255)
gTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc                                                                                                                                                      >  2:791548/1‑150 (MQ=255)
    ttCATACGCAGCAGCGGCTGCGTATGAAAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGAc                                                                                                                                                      >  1:1672497/1‑150 (MQ=255)
     tCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGGTAACCGGGGCGTCCGGGTCTCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa                                                                                                                                                 >  2:253762/1‑150 (MQ=255)
      ttgaCGCAGCCGCTGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCAACCACCCGTCCCTGCTGCGTTCCGACCAGAAc                                                                                                                                                >  2:264755/4‑150 (MQ=255)
        tACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgc                                                                                                                                              >  2:1524313/1‑150 (MQ=255)
         aCGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgcc                                                                                                                                             <  2:178941/150‑1 (MQ=255)
          cGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGCAcgccg                                                                                                                                            >  1:266497/1‑150 (MQ=255)
          cGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccg                                                                                                                                            >  1:413752/1‑150 (MQ=255)
           gcagcaGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccgc                                                                                                                                           <  1:508595/150‑1 (MQ=255)
           gcagcaGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccgc                                                                                                                                           <  2:1384386/150‑1 (MQ=255)
             agcagcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagcag                                                                                                                                         <  2:726954/150‑2 (MQ=255)
             agcagcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTAAATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTAAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCaa                                                                                                                                         >  1:1110041/1‑150 (MQ=255)
              gcagcGGCTGCGG‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagcagg                                                                                                                                        <  2:1256429/150‑3 (MQ=255)
               cagcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAAcc                                                                                                                                       <  1:884105/150‑1 (MQ=255)
               cagcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACATCCCTGCTGCGTTCCGACCAGAACGCCGCAAcc                                                                                                                                       <  2:460261/150‑1 (MQ=255)
                agcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCa                                                                                                                                      <  1:791548/150‑1 (MQ=255)
                 gcGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAg                                                                                                                                     >  2:203062/1‑150 (MQ=255)
                   ggCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTa                                                                                                                                   >  2:1525770/1‑150 (MQ=255)
                    gCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTAc                                                                                                                                  >  1:606464/1‑150 (MQ=255)
                     cTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTAAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTAcc                                                                                                                                 <  2:884950/150‑1 (MQ=255)
                      tGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTTATGCAGTTCCGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGCCCAGAACGCCGCAACCAGTACCg                                                                                                                                >  2:308333/1‑150 (MQ=255)
                        cGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGAAGTACAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCa                                                                                                                              >  1:1074043/1‑150 (MQ=255)
                         gT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGATCGCCGCAACCAGTACCGCAg                                                                                                                             <  1:216116/150‑1 (MQ=255)
                         gT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAg                                                                                                                             <  2:238272/150‑1 (MQ=255)
                         gT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAg                                                                                                                             <  2:243137/150‑1 (MQ=255)
                          t‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGa                                                                                                                            <  2:1156858/150‑1 (MQ=255)
                                aaaCGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGaaa                                                                                                                          >  2:247776/1‑150 (MQ=255)
                                aaaCGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGGGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGaaa                                                                                                                          <  1:983378/150‑1 (MQ=255)
                                     gTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAacggca                                                                                                                     >  2:1386785/1‑150 (MQ=255)
                                        tttACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAacggcacgg                                                                                                                  <  1:1328596/150‑1 (MQ=255)
                                            ccGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCt                                                                                                              >  1:493348/1‑150 (MQ=255)
                                             cGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCtt                                                                                                             >  2:181928/1‑150 (MQ=255)
                                              gTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCttt                                                                                                            <  2:1706263/150‑1 (MQ=255)
                                                 cATGATTTTGCGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACGGGTGCGTCCGGGTCGCCCTGACAGGCCACACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATcc                                                                                                        <  2:11521/150‑1 (MQ=255)
                                                   tGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATccc                                                                                                       <  2:531440/150‑1 (MQ=255)
                                                       tttGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACTGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCAttt                                                                                                   <  1:851845/150‑1 (MQ=255)
                                                           ggCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATg                                                                                               <  1:1344312/150‑1 (MQ=255)
                                                               ggAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCACCCAGTACCGCAGAAACGGCCCGGCGCTTTATCCCATTTGATGGCtt                                                                                           >  2:823410/1‑150 (MQ=255)
                                                                 aaaTAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTcc                                                                                         <  2:600661/150‑1 (MQ=255)
                                                                  aaTAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGACACGGCGCTTTATCCCATTAGATGGCTTCct                                                                                        >  1:1550668/1‑150 (MQ=255)
                                                                   aTAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCctg                                                                                       <  1:342100/150‑1 (MQ=255)
                                                                       cccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGa                                                                                   <  1:137761/150‑1 (MQ=255)
                                                                       cccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGa                                                                                   <  1:1465796/150‑1 (MQ=255)
                                                                        ccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGGTGAc                                                                                  >  2:104810/1‑150 (MQ=255)
                                                                        ccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGAc                                                                                  >  2:680435/1‑150 (MQ=255)
                                                                        ccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGAc                                                                                  <  2:606464/150‑1 (MQ=255)
                                                                           tAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACTGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGAc                                                                                  <  1:1100273/147‑1 (MQ=255)
                                                                           tAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGAc                                                                                  >  2:1100273/1‑147 (MQ=255)
                                                                           tAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACcaa                                                                               <  1:1524313/150‑1 (MQ=255)
                                                                              tGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACcaacaa                                                                            >  2:234029/1‑150 (MQ=255)
                                                                               gCAGTTTAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCaacaac                                                                           >  1:939116/1‑150 (MQ=255)
                                                                               gCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCaacaac                                                                           >  2:1481763/1‑150 (MQ=255)
                                                                                 aGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACgg                                                                         >  1:1420296/1‑150 (MQ=255)
                                                                                    tCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAGCGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGgcgc                                                                      >  2:1530671/1‑150 (MQ=255)
                                                                                      aGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGAACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCgg                                                                    <  2:365198/150‑1 (MQ=255)
                                                                                       ggCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCggg                                                                   >  1:448122/1‑150 (MQ=255)
                                                                                              gTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCACCAACGGCGCGGGCAACGcc                                                            >  2:1583821/1‑150 (MQ=255)
                                                                                                tAACCGTTCCGCCCGGGTGGCCCTGAAAGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCgg                                                         <  1:781788/150‑1 (MQ=255)
                                                                                                  aCCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGc                                                        >  2:407317/1‑150 (MQ=255)
                                                                                                   ccGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCa                                                       >  1:1464143/1‑150 (MQ=255)
                                                                                                        gCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTg                                                  <  2:939116/150‑1 (MQ=255)
                                                                                                          gTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTgag                                                <  1:483105/150‑1 (MQ=255)
                                                                                                             ctggTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCCGAGCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                             <  1:285078/148‑1 (MQ=255)
                                                                                                             cGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                             <  1:217189/150‑1 (MQ=255)
                                                                                                             cGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                             <  1:391459/150‑1 (MQ=255)
                                                                                                               ggTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCgg                                           <  1:994543/150‑1 (MQ=255)
                                                                                                                  cGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAg                                        >  1:642158/1‑150 (MQ=255)
                                                                                                                    cccTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAgcc                                      >  2:1286075/1‑150 (MQ=255)
                                                                                                                    cccTGACAGGCCACCACACGTCCCTGCTGCGTTACGACCAGAACGCCGCCACCAGTACCGCAGAAACGGCACGGCGCGTTATCCCATTTGAGGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCAGGCTGGGCCCGGCAgcc                                      >  1:1643540/1‑150 (MQ=255)
                                                                                                                          cAGTCCACCACACGACCCTGCTGCGCTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGTCGCGTTATCCCATTTGATGGCTTCCTGCTGACCAACTACGGCTCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCa                                <  1:264755/150‑1 (MQ=255)
                                                                                                                             gCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCg                             <  2:448122/150‑1 (MQ=255)
                                                                                                                              ccaccaCACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCACTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCgg                            >  1:1702369/1‑150 (MQ=255)
                                                                                                                               caccacACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGc                           >  2:1172103/1‑150 (MQ=255)
                                                                                                                                      cGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACACCGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACg                    >  2:65266/1‑150 (MQ=255)
                                                                                                                                      cGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACg                    <  1:979557/150‑1 (MQ=255)
                                                                                                                                        tCCCTGCTGCGTTCCGTCCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGc                  <  2:1074043/150‑1 (MQ=255)
                                                                                                                                         cccTGCTGCGTTCCGACCAGACCGCCGCAACCAGTACCGCAGAAACGGCACGGCGCGTAATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCg                 <  1:1530671/150‑1 (MQ=255)
                                                                                                                                           ctgctgCGTTCCGACCAGACCGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGtt               <  1:1598094/150‑1 (MQ=255)
                                                                                                                                            tgctgcGTTCCGACCAGAATGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTa              >  2:1115259/1‑150 (MQ=255)
                                                                                                                                            tgcggcGTTCTGGTCGGAACGCAGCACCCAGTACCGCAGAAACGGCACGGCGCTTTATCCCAGTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCCGCGGCCGCCACGGCGTTa              >  1:1512974/1‑150 (MQ=255)
                                                                                                                                            tgcggcGTTCTGGTCGGAACGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTa              >  1:581046/1‑150 (MQ=255)
                                                                                                                                            tgcggcGTTCTGGTCGGAACGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTa              >  2:1123820/1‑150 (MQ=255)
                                                                                                                                            tgcggcGTTCTGGTCGGAACGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCAATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTa              >  1:1506460/1‑150 (MQ=255)
                                                                                                                                             gctgcGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAg             >  2:900255/1‑150 (MQ=255)
                                                                                                                                              ctgcGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGc            >  1:16355/1‑150 (MQ=255)
                                                                                                                                              ctgcGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGc            >  2:495165/1‑150 (MQ=255)
                                                                                                                                                gcGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCTGCCGCAACGGCGTTAGCtt          <  1:1339559/150‑1 (MQ=255)
                                                                                                                                                gcGTCCCGACCAGAACGCCGCAACCAGGACCGCAGATCCGGCACGGCGCGTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCtt          <  1:1082000/150‑1 (MQ=255)
                                                                                                                                                 cgatCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCGTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCttt         <  1:140895/147‑1 (MQ=255)
                                                                                                                                                  gTTCCGACCAGAAAGCCGCAACCATTACCGCAGAAAGGGCACGGCGCTTTAACCCATTTGATGGCTTCCTGCTGACCAACATCGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGAGGCCGCAACGGCGTTAGCTTTc        <  2:1537500/150‑1 (MQ=255)
                                                                                                                                                     ccGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCGTTCAtt     >  2:1500519/1‑149 (MQ=255)
                                                                                                                                                      cGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGCCCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGCCCGGCAGCCGCCGCCGCGGCCGCAACGGCGTTAGCTTTCAtca    >  1:820246/1‑148 (MQ=255)
                                                                                                                                                       gACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATaaa   >  1:855153/1‑150 (MQ=255)
                                                                                                                                                       gACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATaaa   >  2:346715/1‑150 (MQ=255)
                                                                                                                                                       gACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCCGCAGCCGCCGCAGCGGCCGCCACGCCGATAGCTGTCATaa  >  1:272785/1‑149 (MQ=255)
                                                                                                                                                        aCCAGACCGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAg  <  2:167792/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAG  >  NC_000913/2302436‑2302733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: