Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,230 C→G 11.3% M277I (ATG→ATC) ‡ narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2300CG11.3% 167.6 / 14.9 62M277I (ATG→ATC) ‡narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base G (4/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543092‑1543379
                                                                                                                                          |                                                                                                                                                      
aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaaaaac                                                                                                                                            <  1:490750/150‑3 (MQ=255)
aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGAGATGGCACCGATAAATGGGCAAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCaaaaac                                                                                                                                            <  2:823939/150‑3 (MQ=255)
 aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGCAAAAacc                                                                                                                                           <  1:1082674/150‑1 (MQ=255)
            cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                <  1:781309/150‑1 (MQ=255)
            cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGCGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAAcc                                                                                                                                <  1:848774/150‑1 (MQ=255)
                 cccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCCCCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATAGTCCCATCCGGGAACTGGGTTTTTGCCAGCATGGCACCCCCGCACAAAAACCGATaaa                                                                                                                          >  1:134841/1‑150 (MQ=255)
                    gCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACg                                                                                                                        <  1:1314671/150‑1 (MQ=255)
                          cgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCACAACCCGCAGAAAAACCGATAACCGAGCCGa                                                                                                                  >  1:365085/1‑150 (MQ=255)
                           gccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGaa                                                                                                                 >  1:142352/1‑150 (MQ=255)
                                 cTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGc                                                                                                           >  2:347327/1‑150 (MQ=255)
                                      cGGAATTAGCACCCCCAACCGACCGCGCGATGGCACCGATTAGTGGGCCATTGAACGGCAGGCCCAGACTCTTCACATCCGGGAACTGGGTTTTTCCCTGCAGGGCAATCCCCTCATAAAATCCGATACACGAACCGAAGGAGGCAAggg                                                                                                      >  2:851621/1‑149 (MQ=255)
                                      cGGAAATAGCACCACCAACCGACCGCGCTAAGGCAACGATAAATTGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGt                                                                                                      <  1:243452/150‑1 (MQ=255)
                                       ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGATCTGGGTTTTTGCCAGCATGGCAAAAGCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  1:83854/1‑150 (MQ=255)
                                       ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGtt                                                                                                     >  2:342691/1‑149 (MQ=255)
                                         aaaTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGAAAAACAAACCGAAGGTGGCAAGGTaaa                                                                                                   <  2:393903/150‑1 (MQ=255)
                                            tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATACACGAACCGAAGGTGGCAAGGTAAAGcc                                                                                                >  1:1146159/1‑149 (MQ=255)
                                            tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                                >  1:717194/1‑150 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAACACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               >  1:1267119/1‑150 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  2:690124/150‑1 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  1:1299836/150‑1 (MQ=255)
                                                   accaACCGACAGCGCGATGGCACCGATAAATGGGCCAAACACCGCCCGGCGCAGAATAGTCACATCCGGGAACTGGGTTTATGCCAGCATGGAAAAACACGCAGAAAAACCGAAAAACGCACCCAAGGTGGCCAGGTCAAGCAGGcgccg                                                                                         >  2:810583/1‑146 (MQ=255)
                                                     caACCGACCGCGCGATGGCACCGATAAATGGGCCACAGAACGCCAGGCGCAGAATATTCACATCCGGGAACGGGCTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAATGTGGCAAGGTAAAGCAGGCTcagca                                                                                       <  1:794201/150‑1 (MQ=255)
                                                     caACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGACCTGGGTTTATGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTccgcc                                                                                       >  2:811115/1‑149 (MQ=255)
                                                            ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                <  1:904631/150‑1 (MQ=255)
                                                                 cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                           >  2:548894/1‑150 (MQ=255)
                                                                  gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAg                                                                          >  2:1036242/1‑150 (MQ=255)
                                                                   aTGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAgg                                                                         >  1:513020/1‑150 (MQ=255)
                                                                       cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGAAAGGTAAAGCAGGATCAGAAGCCAGAGATGCAGGCGtt                                                                     >  2:477482/1‑150 (MQ=255)
                                                                                aTGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  2:1135611/150‑1 (MQ=255)
                                                                                aTGGGCCAAAGAACGCCAGTCGCAGATTATTCACATCCGGGAACTGGGTTTTTGCCTGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  1:995087/150‑1 (MQ=255)
                                                                                 tGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                           <  2:1003134/150‑1 (MQ=255)
                                                                                   ggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                         >  2:1246214/1‑150 (MQ=255)
                                                                                      cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                      <  2:1076878/150‑1 (MQ=255)
                                                                                          gAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTc                                                  >  1:1091023/1‑150 (MQ=255)
                                                                                            aCGCCAGGCGCAGAATATTCACATCCGGGAACTGTGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCgg                                                <  1:848695/150‑1 (MQ=255)
                                                                                               ccAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCCGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             <  2:807116/150‑1 (MQ=255)
                                                                                                cAGGCGCAGAATATTCACATCCTGGATCTTGGTTTTTGCCAGCATTGCAAAACCCGCAGAAAAACCTATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCATCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAAt                                            <  1:1003710/150‑1 (MQ=255)
                                                                                                cAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACGGGGAGGTGGTCGGCatt                                            >  2:47720/1‑148 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:991960/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:986042/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  2:174882/1‑150 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:365075/150‑1 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:347327/150‑1 (MQ=255)
                                                                                                               cacaTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                             >  1:365902/1‑150 (MQ=255)
                                                                                                                  ctCCGGGAACTGGGTTTTTGCCAGCCTGGCAAACCCCGCAGAAAAACCGATAACCGAACCGAAGGGGGCAAGGGCAAGCAGGCACAGCAGCCAGAGATGCAGGCGTTGTAAGACAGCGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          >  2:930396/2‑150 (MQ=255)
                                                                                                                  aTCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          <  1:586029/150‑1 (MQ=255)
                                                                                                                     cGGCAACGGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAACCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGat                       <  1:954228/150‑1 (MQ=255)
                                                                                                                      gggAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGata                      <  1:1193973/150‑1 (MQ=255)
                                                                                                                           tgcgggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                  >  2:621510/7‑150 (MQ=255)
                                                                                                                            tGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc                <  1:477503/150‑1 (MQ=255)
                                                                                                                              gggttttGCCATGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  2:1118953/4‑150 (MQ=255)
                                                                                                                              ggTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  1:671490/1‑150 (MQ=255)
                                                                                                                               ggttttGCCATGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             >  1:636112/3‑150 (MQ=255)
                                                                                                                                tttttGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            <  2:642606/150‑1 (MQ=255)
                                                                                                                                gttttGCCATGCTGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            >  2:1116637/2‑150 (MQ=255)
                                                                                                                                     gCCAGCATGGCAAAACCCGCAGAACAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc       >  1:749867/1‑150 (MQ=255)
                                                                                                                                      ccAGCATGGCAAAACCCGCAGAAAAACCGATAAACGCACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGAGTTGTATGACAGGGAGCTGGGCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAca      >  2:1218421/1‑149 (MQ=255)
                                                                                                                                      ccAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCCGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGTTATCATTCCTCCCTGaac      >  2:200322/1‑148 (MQ=255)
                                                                                                                                      ccAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTATGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      >  2:1213583/1‑150 (MQ=255)
                                                                                                                                      ccAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      <  1:716108/150‑1 (MQ=255)
                                                                                                                                       cAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAACGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTCGTCGGCAATTGCGGCGCGTGAACTGGCGAAATCATTCATCCCAGACCa     >  1:217151/1‑150 (MQ=255)
                                                                                                                                         gCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATg   >  2:808100/1‑150 (MQ=255)
                                                                                                                                          cATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGCGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgc  >  1:362345/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                      
AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGCAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543092‑1543379

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