Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,282,840 G→T 53.5% E993* (GAG→TAG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,282,8400GT53.5% 15.7 / 197.6 131E993* (GAG→TAG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (26/35);  new base T (36/34);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 138 total aligned reads displayed.

GCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCC  >  NC_000913/1282741‑1282938
                                                                                                   |                                                                                                    
gCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACg                                                                                                      <  1:3125883/100‑1 (MQ=255)
gCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACg                                                                                                      >  2:37731/1‑100 (MQ=255)
  aCCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCGGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgag                                                                                                    <  1:1882185/100‑1 (MQ=255)
   ccGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACtaga                                                                                                   >  1:209064/1‑100 (MQ=255)
   ccGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACgaga                                                                                                   <  2:2478922/100‑1 (MQ=255)
    cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTTGTCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGagaa                                                                                                  <  1:3613757/100‑1 (MQ=255)
    cGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTagaa                                                                                                  <  1:669589/100‑1 (MQ=255)
     ggAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTagaag                                                                                                 >  1:2902839/1‑100 (MQ=255)
      gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGagaaga                                                                                                <  2:3852678/100‑1 (MQ=255)
      gAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGagaaga                                                                                                <  2:3455088/100‑1 (MQ=255)
         aCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATcc                                                                                             <  1:437905/100‑1 (MQ=255)
         aCCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATcc                                                                                             <  2:602867/100‑1 (MQ=255)
          ccAACGGTCAGGTAGCCGTGAACGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATCCGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCg                                                                                            >  1:4261026/1‑100 (MQ=255)
           cAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGc                                                                                           >  2:187279/1‑100 (MQ=255)
            aaCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCt                                                                                          <  2:2362129/100‑1 (MQ=255)
            aaCGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCt                                                                                          <  1:3859868/100‑1 (MQ=255)
              cGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTc                                                                                        >  2:2677603/1‑100 (MQ=255)
               ggTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTcc                                                                                       >  2:3838021/1‑100 (MQ=255)
                gTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAGTTTACCGGTCGTGACCACACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCcg                                                                                      >  1:1053234/1‑100 (MQ=255)
                  cAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCcgcg                                                                                    >  2:1343975/1‑100 (MQ=255)
                  cAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCcgcg                                                                                    <  1:2370922/100‑1 (MQ=255)
                  cAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcgcg                                                                                    <  2:1973273/100‑1 (MQ=255)
                    ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGat                                                                                  >  2:4267515/1‑100 (MQ=255)
                    ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat                                                                                  >  2:3343058/1‑100 (MQ=255)
                    ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat                                                                                  <  1:516/100‑1 (MQ=255)
                    ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat                                                                                  >  1:4110269/1‑100 (MQ=255)
                    ggTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat                                                                                  >  1:4109838/1‑100 (MQ=255)
                      tAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGatat                                                                                >  2:3276397/1‑100 (MQ=255)
                        gCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTc                                                                              >  2:2294478/1‑100 (MQ=255)
                          cGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAg                                                                            <  2:1744703/100‑1 (MQ=255)
                           gTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAgg                                                                           <  2:2931665/100‑1 (MQ=255)
                            tGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGc                                                                          <  2:1367064/100‑1 (MQ=255)
                            tGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGc                                                                          <  2:3790469/100‑1 (MQ=255)
                             gAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGca                                                                         <  2:3438644/100‑1 (MQ=255)
                                aGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAg                                                                      >  2:78037/1‑100 (MQ=255)
                                aGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAg                                                                      >  1:1802853/1‑100 (MQ=255)
                                  ccTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGcc                                                                    <  1:1005215/100‑1 (MQ=255)
                                  ccTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGcc                                                                    >  1:41424/1‑100 (MQ=255)
                                  ccTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGcc                                                                    >  1:312741/1‑100 (MQ=255)
                                    tGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgc                                                                  >  2:2749322/1‑100 (MQ=255)
                                     gggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgcg                                                                 <  2:3115686/100‑1 (MQ=255)
                                      ggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgcgc                                                                <  2:4110562/100‑1 (MQ=255)
                                      ggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgcgc                                                                <  2:4110291/100‑1 (MQ=255)
                                      ggCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCgcgc                                                                <  2:4109860/100‑1 (MQ=255)
                                       gCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCa                                                               >  1:340786/1‑100 (MQ=255)
                                        cTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaa                                                              >  2:1222632/1‑100 (MQ=255)
                                          gcgcTGAGCGGATTTACCGGTCGTGAGCTTACGCATATGGCGCTGGATAAAGAAGACTAGAAGATCCGCCTCCGCGATGTTCATGCACGGCCgcgc                                                                >  1:1443202/4‑96 (MQ=255)
                                           cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAAt                                                           <  2:2756233/100‑1 (MQ=255)
                                           cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAAt                                                           >  2:1193860/1‑100 (MQ=255)
                                            ccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAAtt                                                          <  2:3562262/100‑1 (MQ=255)
                                            ccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAAtt                                                          >  1:1608478/1‑100 (MQ=255)
                                                 gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtctc                                                     >  2:2559227/1‑100 (MQ=255)
                                                 gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtctc                                                     <  2:560906/100‑1 (MQ=255)
                                                   gAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTcc                                                    <  1:2853474/99‑1 (MQ=255)
                                                        tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg                                            >  2:2141903/1‑100 (MQ=255)
                                                        tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg                                            >  2:3649631/1‑100 (MQ=255)
                                                        tACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCg                                            >  1:3185779/1‑100 (MQ=255)
                                                         aCCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa                                           >  1:363248/1‑100 (MQ=255)
                                                          ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGAc                                          >  1:522640/1‑100 (MQ=255)
                                                           cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGAcc                                         <  2:2621665/100‑1 (MQ=255)
                                                           cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGAcc                                         <  2:609818/100‑1 (MQ=255)
                                                            ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACct                                        >  1:889251/1‑100 (MQ=255)
                                                            ggTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACct                                        >  2:2090239/1‑100 (MQ=255)
                                                             gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctg                                       >  2:3463277/1‑100 (MQ=255)
                                                             gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctg                                       >  2:4255714/1‑100 (MQ=255)
                                                             gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCGTCCGCGATATTCAGGCACAGCCGCGCAAAATTATATCCAGCCCGACctg                                       <  1:3944091/100‑1 (MQ=255)
                                                              tCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctgg                                      <  1:3163960/100‑1 (MQ=255)
                                                                gTGACCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtc                                    >  1:72977/1‑100 (MQ=255)
                                                                 tGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtct                                   <  1:75615/100‑1 (MQ=255)
                                                                   aCCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg                                 <  1:3652120/100‑1 (MQ=255)
                                                                   aCCATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg                                 >  1:3459773/1‑100 (MQ=255)
                                                                     cATACGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc                               <  1:3200312/100‑1 (MQ=255)
                                                                       tACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctggtctg                             <  1:420804/100‑1 (MQ=255)
                                                                        aCGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGtctggtctgg                            <  1:1276732/100‑1 (MQ=255)
                                                                         cGCATCTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGa                           >  1:3310009/1‑100 (MQ=255)
                                                                          gCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGACCAGACcaggtcgggctggag                           <  2:2691029/100‑10 (MQ=38)
                                                                            aTCTGGCTCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAg                         <  1:2387710/99‑1 (MQ=255)
                                                                              cTGGCGCTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCACAAGCCCGACCTGGTCTGGTCTGGAAGAt                       >  1:1370955/1‑99 (MQ=255)
                                                                              cTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATg                      <  2:1823693/100‑1 (MQ=255)
                                                                              cTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCCGGTCTGGAAGATg                      <  2:2873262/100‑1 (MQ=255)
                                                                               tGGCGCTGAATAAAGAAGACTAGGAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGa                     <  1:3104532/100‑1 (MQ=255)
                                                                               tGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGa                     <  1:2749186/100‑1 (MQ=255)
                                                                                ggCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGACCAGACCAGGTCGGGCTGGagataat                     <  1:419554/99‑7 (MQ=18)
                                                                                ggCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                    >  2:2573659/1‑100 (MQ=255)
                                                                                ggCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                    >  2:715309/1‑100 (MQ=255)
                                                                                   gcTGAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCCGGTCTGGAAGATGAacac                 >  1:3852157/1‑100 (MQ=255)
                                                                                   gcTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac                 >  1:1755281/1‑100 (MQ=255)
                                                                                    cTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACg                <  2:2130430/100‑1 (MQ=255)
                                                                                     tGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt               >  1:3485577/1‑100 (MQ=255)
                                                                                      gAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGtt              >  1:1717984/1‑100 (MQ=255)
                                                                                      gAATAAAGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGtt              <  2:3852174/100‑1 (MQ=255)
                                                                                        aTAAAGAAGACGAGAGGATCCGCTTCCGCTATATTCAGGCACAGCCGCGCAAAATTATCGCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc            <  1:959091/100‑1 (MQ=255)
                                                                                          aaaGAAGACTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCtt          <  2:1331947/100‑1 (MQ=255)
                                                                                            agaagaCTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc        <  1:3077986/100‑1 (MQ=255)
                                                                                            agaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc        <  1:4231982/100‑1 (MQ=255)
                                                                                            agaagaCGAGAAGATCCGCTTCCGCGATAGTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTAc        <  1:2282196/100‑1 (MQ=255)
                                                                                             gaagaCTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa       >  2:1683637/1‑100 (MQ=255)
                                                                                             gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa       <  2:3602297/100‑1 (MQ=255)
                                                                                              aagaCTAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACaa      <  2:1201632/100‑1 (MQ=255)
                                                                                                  cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc  >  1:2591584/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GCACCGGAAACCAACGGTCAGGTAGCCGTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCC  >  NC_000913/1282741‑1282938

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