Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,229 G→T 12.9% L278M (CTG→ATG)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,543,2290GT12.9% 130.9 / 11.4 61L278M (CTG→ATG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base T (5/3);  total (34/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 = 6.25e-01

AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATG  >  NC_000913/1543092‑1543378
                                                                                                                                         |                                                                                                                                                      
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                              >  2:174882/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                              >  1:991960/1‑150 (MQ=255)
                                                                                                             ttCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGaa                              >  1:986042/1‑150 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                             <  1:347327/150‑1 (MQ=255)
                                                                                                              tCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAAc                             <  1:365075/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)
                                                                                                                                         |                                                                                                                                                      
AAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGACCATG  >  NC_000913/1543092‑1543378

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: ‑