Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,302,582 C→T 6.4% G58R (GGA→AGA)  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,5820CT 6.4% 294.3 / 5.8 93G58R (GGA→AGA) napAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base C (44/43);  new base T (4/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

ACCGTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAA  >  NC_000913/2302433‑2302731
                                                                                                                                                         |                                                                                                                                                     
aCCGTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTcc                                                                                                                                                       >  1:211592/1‑150 (MQ=255)
 ccGTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCg                                                                                                                                                      >  2:342100/1‑150 (MQ=255)
  cGTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGTGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCgc                                                                                                                                                     >  2:633350/1‑149 (MQ=255)
   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:606464/150‑1 (MQ=255)
                                                                           ccTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGAc                                                                                >  2:680435/1‑150 (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: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          >  2:495165/1‑150 (MQ=255)
                                                                                                                                                 ctgcGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGc          >  1:16355/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)
                                                                                                                                                         |                                                                                                                                                     
ACCGTTTTTCATACGCAGCAGCGGCTGCGT‑‑‑‑CAAACGGTCTTTACCGTACATGATTTTGGGCAGGAAATAGCCCTTAATGCAGTTCAGGCCACGGTTAACCGGTGCGTCCGGGTCGCCCTGACAGGCCACCACACGTCCCTGCTGCGTTCCGACCAGAACGCCGCAACCAGTACCGCAGAAACGGCACGGCGCTTTATCCCATTTGATGGCTTCCTGCTGACCAACAACGGCGCGGGCAACGCCCGGCACGCTGAGACCGGCAGCCGCCGCAGCGGCCGCAACGGCGTTAGCTTTCATAA  >  NC_000913/2302433‑2302731

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: