Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,231 A→C 12.3% M277R (ATG→AGG) ‡ narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2310AC12.3% 154.2 / 11.4 63M277R (ATG→AGG) ‡narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (31/24);  new base C (5/3);  total (36/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 = 9.77e-01

AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  NC_000913/1543092‑1543379
                                                                                                                                           |                                                                                                                                                     
aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaac                                                                                                                                            <  1:490750/150‑3 (MQ=255)
aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGAGATGGCACCGATAAATGGGCAAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaac                                                                                                                                            <  2:823939/150‑3 (MQ=255)
 aaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑AAAAacc                                                                                                                                           <  1:1082674/150‑1 (MQ=255)
            cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAcc                                                                                                                                <  1:781309/150‑1 (MQ=255)
            cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGCGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAcc                                                                                                                                <  1:848774/150‑1 (MQ=255)
                 cccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCCCCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATAGTCCCATCCGGGAACTGGGTTTTTGCCAGCATGGC‑‑ACCCCCGCACAAAAACCGATaaa                                                                                                                          >  1:134841/1‑150 (MQ=255)
                    gCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACg                                                                                                                        <  1:1314671/150‑1 (MQ=255)
                          cgccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑ACAACCCGCAGAAAAACCGATAACCGAGCCGa                                                                                                                  >  1:365085/1‑150 (MQ=255)
                           gccgAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGaa                                                                                                                 >  1:142352/1‑150 (MQ=255)
                                 cTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGc                                                                                                           >  2:347327/1‑150 (MQ=255)
                                      cGGAATTAGCACCCCCAACCGACCGCGCGATGGCACCGATTAGTGGGCCATTGAACGGCAGGCCCAGACTCTTCACATCCGGGAACTGGGTTTTTCCCTGCAGGGC‑AATCCCCTCATAAAATCCGATACACGAACCGAAGGAGGCAAggg                                                                                                      >  2:851621/1‑149 (MQ=255)
                                      cGGAAATAGCACCACCAACCGACCGCGCTAAGGCAACGATAAATTGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGt                                                                                                      <  1:243452/150‑1 (MQ=255)
                                       ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGATCTGGGTTTTTGCCAGCATGGC‑AAAAGCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                     >  1:83854/1‑150 (MQ=255)
                                       ggAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGtt                                                                                                     >  2:342691/1‑149 (MQ=255)
                                         aaaTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGAAAAACAAACCGAAGGTGGCAAGGTaaa                                                                                                   <  2:393903/150‑1 (MQ=255)
                                            tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATACACGAACCGAAGGTGGCAAGGTAAAGcc                                                                                                >  1:1146159/1‑149 (MQ=255)
                                            tAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                                >  1:717194/1‑150 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AACACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               >  1:1267119/1‑150 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  2:690124/150‑1 (MQ=255)
                                             aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                               <  1:1299836/150‑1 (MQ=255)
                                                   accaACCGACAGCGCGATGGCACCGATAAATGGGCCAAACACCGCCCGGCGCAGAATAGTCACATCCGGGAACTGGGTTTATGCCAGCATGGA‑AAAACACGCAGAAAAACCGAAAAACGCACCCAAGGTGGCCAGGTCAAGCAGGcgccg                                                                                         >  2:810583/1‑146 (MQ=255)
                                                     caACCGACCGCGCGATGGCACCGATAAATGGGCCACAGAACGCCAGGCGCAGAATATTCACATCCGGGAACGGGCTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAATGTGGCAAGGTAAAGCAGGCTcagca                                                                                       <  1:794201/150‑1 (MQ=255)
                                                     caACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGACCTGGGTTTATGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTccgcc                                                                                       >  2:811115/1‑149 (MQ=255)
                                                            ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                                <  1:904631/150‑1 (MQ=255)
                                                                 cgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                           >  2:548894/1‑150 (MQ=255)
                                                                  gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAg                                                                          >  2:1036242/1‑150 (MQ=255)
                                                                   aTGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAgg                                                                         >  1:513020/1‑150 (MQ=255)
                                                                       cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGAAAGGTAAAGCAGGATCAGAAGCCAGAGATGCAGGCGtt                                                                     >  2:477482/1‑150 (MQ=255)
                                                                                aTGGGCCAAAGACCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGA‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  2:1135611/150‑1 (MQ=255)
                                                                                aTGGGCCAAAGAACGCCAGTCGCAGATTATTCACATCCGGGAACTGGGTTTTTGCCTGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAg                                                            <  1:995087/150‑1 (MQ=255)
                                                                                 tGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                           <  2:1003134/150‑1 (MQ=255)
                                                                                   ggCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                         >  2:1246214/1‑150 (MQ=255)
                                                                                      cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGGTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                      <  2:1076878/150‑1 (MQ=255)
                                                                                          gAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTc                                                  >  1:1091023/1‑150 (MQ=255)
                                                                                            aCGCCAGGCGCAGAATATTCACATCCGGGAACTGTGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCgg                                                <  1:848695/150‑1 (MQ=255)
                                                                                               ccAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCCGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                             <  2:807116/150‑1 (MQ=255)
                                                                                                cAGGCGCAGAATATTCACATCCTGGATCTTGGTTTTTGCCAGCATTGC‑AAAACCCGCAGAAAAACCTATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCATCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAAt                                            <  1:1003710/150‑1 (MQ=255)
                                                                                                cAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACGGGGAGGTGGTCGGCatt                                            >  2:47720/1‑148 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:991960/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  1:986042/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                               >  2:174882/1‑150 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:365075/150‑1 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                              <  1:347327/150‑1 (MQ=255)
                                                                                                               cacaTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACt                             >  1:365902/1‑150 (MQ=255)
                                                                                                                  ctCCGGGAACTGGGTTTTTGCCAGCCTGGC‑AAACCCCGCAGAAAAACCGATAACCGAACCGAAGGGGGCAAGGGCAAGCAGGCACAGCAGCCAGAGATGCAGGCGTTGTAAGACAGCGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          >  2:930396/2‑150 (MQ=255)
                                                                                                                  aTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                          <  1:586029/150‑1 (MQ=255)
                                                                                                                     cGGCAACGGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAACCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGat                       <  1:954228/150‑1 (MQ=255)
                                                                                                                      gggAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGata                      <  1:1193973/150‑1 (MQ=255)
                                                                                                                           tgcgggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat                  >  2:621510/7‑150 (MQ=255)
                                                                                                                            tGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattc                <  1:477503/150‑1 (MQ=255)
                                                                                                                              gggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  2:1118953/4‑150 (MQ=255)
                                                                                                                              ggTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat              >  1:671490/1‑150 (MQ=255)
                                                                                                                               ggttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc             >  1:636112/3‑150 (MQ=255)
                                                                                                                                tttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            <  2:642606/150‑1 (MQ=255)
                                                                                                                                gttttGCCATGCTGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc            >  2:1116637/2‑150 (MQ=255)
                                                                                                                                     gCCAGCATGGC‑AAAACCCGCAGAACAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc       >  1:749867/1‑150 (MQ=255)
                                                                                                                                      ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGCACCGAAGGGGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGAGTTGTATGACAGGGAGCTGGGCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAca      >  2:1218421/1‑149 (MQ=255)
                                                                                                                                      ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCCGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGTTATCATTCCTCCCTGaac      >  2:200322/1‑148 (MQ=255)
                                                                                                                                      ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTATGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      >  2:1213583/1‑150 (MQ=255)
                                                                                                                                      ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAcc      <  1:716108/150‑1 (MQ=255)
                                                                                                                                       cAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAACGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTCGTCGGCAATTGCGGCGCGTGAACTGGCGAAATCATTCATCCCAGACCa     >  1:217151/1‑150 (MQ=255)
                                                                                                                                         gCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATg   >  2:808100/1‑150 (MQ=255)
                                                                                                                                          cATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGCGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATgc  >  1:362345/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATGC  >  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: ‑