Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,587 A→G 6.0% L56P (CTG→CCG)  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,5870AG 6.0% 291.4 / 7.2 99L56P (CTG→CCG) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base A (46/47);  new base G (4/2);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 102 total aligned reads displayed.

TCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCTA  >  NC_000913/2302441‑2302736
                                                                                                                                                  |                                                                                                                                                     
tCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGGTAACCGGGGCGTCCGGGTCTCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGaa                                                                                                                                                    >  2:253762/1‑150 (MQ=255)
 ttgaCGCAGCCGCTGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCAACCACCCGTCCCTGCTGCGTTCCGACCAGAAc                                                                                                                                                   >  2:264755/4‑150 (MQ=255)
   tACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgc                                                                                                                                                 >  2:1524313/1‑150 (MQ=255)
    aCGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgcc                                                                                                                                                <  2:178941/150‑1 (MQ=255)
     cGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGCAcgccg                                                                                                                                               >  1:266497/1‑150 (MQ=255)
     cGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccg                                                                                                                                               >  1:413752/1‑150 (MQ=255)
      gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccgc                                                                                                                                              <  2:1384386/150‑1 (MQ=255)
      gcagcaGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAAcgccgc                                                                                                                                              <  1:508595/150‑1 (MQ=255)
        agcagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagcag                                                                                                                                            <  2:726954/150‑2 (MQ=255)
        agcagcGGCTGCGTCAAACGGTCTTTACCGTAAATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTAAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCaa                                                                                                                                            >  1:1110041/1‑150 (MQ=255)
         gcagcGGCTGCGGCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCTGGTCGGAACgcagcagg                                                                                                                                           <  2:1256429/150‑3 (MQ=255)
          cagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAAcc                                                                                                                                          <  1:884105/150‑1 (MQ=255)
          cagcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACATCCCTGCTGCGTTCCGACCAGAACGCCGCAAcc                                                                                                                                          <  2:460261/150‑1 (MQ=255)
           agcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCa                                                                                                                                         <  1:791548/150‑1 (MQ=255)
            gcGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAg                                                                                                                                        >  2:203062/1‑150 (MQ=255)
              ggCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTa                                                                                                                                      >  2:1525770/1‑150 (MQ=255)
               gCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTAc                                                                                                                                     >  1:606464/1‑150 (MQ=255)
                cTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTAAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTAcc                                                                                                                                    <  2:884950/150‑1 (MQ=255)
                 tGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTTATGCAGTTCCGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGCCCAGAACGCCGCAACCAGTACCg                                                                                                                                   >  2:308333/1‑150 (MQ=255)
                   cGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGAAGTACAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCa                                                                                                                                 >  1:1074043/1‑150 (MQ=255)
                    gTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGATCGCCGCAACCAGTACCGCAg                                                                                                                                <  1:216116/150‑1 (MQ=255)
                    gTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAg                                                                                                                                <  2:238272/150‑1 (MQ=255)
                    gTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAg                                                                                                                                <  2:243137/150‑1 (MQ=255)
                     tCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGa                                                                                                                               <  2:1156858/150‑1 (MQ=255)
                       aaaCGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGaaa                                                                                                                             >  2:247776/1‑150 (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:1465796/150‑1 (MQ=255)
                                                              cccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGa                                                                                      <  1:137761/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)
                                                                                                 gTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTgag                                                   <  1:483105/150‑1 (MQ=255)
                                                                                                    ctggTCGCCCTGACAGGCCACCACACGTCCCTGCGGCGTTCCGAGCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                                <  1:285078/148‑1 (MQ=255)
                                                                                                    cGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                                <  1:391459/150‑1 (MQ=255)
                                                                                                    cGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGAcc                                                <  1:217189/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                 >  2:1123820/1‑150 (MQ=255)
                                                                                                                                   tgcggcGTTCTGGTCGGAACGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTa                 >  1:581046/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      >  2:346715/1‑150 (MQ=255)
                                                                                                                                              gACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATaaa      >  1:855153/1‑150 (MQ=255)
                                                                                                                                              gACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCCGCAGCCGCCGCAGCGGCCGCCACGCCGATAGCTGTCATaac      >  1:272785/1‑149 (MQ=255)
                                                                                                                                               aCCAGACCGCAGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAg     <  2:167792/150‑1 (MQ=255)
                                                                                                                                                ccAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAg     >  2:554886/1‑149 (MQ=255)
                                                                                                                                                ccAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAg     <  1:554886/149‑1 (MQ=255)
                                                                                                                                                ccAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGc    >  1:151947/1‑150 (MQ=255)
                                                                                                                                                 cAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCt   <  1:1096716/150‑1 (MQ=255)
                                                                                                                                                 cAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCt   <  1:1018635/150‑1 (MQ=255)
                                                                                                                                                  aGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGCCCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCGTTCAGAACGCTa  >  2:105910/1‑150 (MQ=255)
                                                                                                                                                  aGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCTa  <  1:1481763/150‑1 (MQ=255)
                                                                                                                                                  aGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCTa  >  2:187084/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCATACGCAGCAGCGGCTGCGTCAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAAAGCTA  >  NC_000913/2302441‑2302736

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: