Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,543,229 G→T 6.4% 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,2290GT 6.4% 551.5 / 20.7 157L278M (CTG→ATG) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (74/73);  new base T (2/8);  total (76/81)
Fisher's exact test for biased strand distribution p-value = 9.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 157 total aligned reads displayed.

CCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAC  >  NC_000913/1543082‑1543374
                                                                                                                                                   |                                                                                                                                                  
ccATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGc                                                                                                                                                   <  3:346594/149‑1 (MQ=255)
  aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                 >  3:1318487/1‑149 (MQ=255)
  aTAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCAt                                                                                                                                                 >  3:159067/1‑149 (MQ=255)
   tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                <  2:1441762/149‑1 (MQ=255)
   tAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATg                                                                                                                                                >  1:659338/1‑149 (MQ=255)
     aaaaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGc                                                                                                                                              >  4:1332514/1‑149 (MQ=255)
       aaTAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aa                                                                                                                                           >  3:784930/1‑149 (MQ=255)
         tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑aaaa                                                                                                                                         <  2:618973/149‑1 (MQ=255)
         tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑aaaa                                                                                                                                         <  3:1031558/149‑1 (MQ=255)
         tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑aaaa                                                                                                                                         >  2:1192196/1‑149 (MQ=255)
         tAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑aaaaa                                                                                                                                        <  2:134372/149‑2 (MQ=255)
          aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaa                                                                                                                                        <  4:784930/149‑2 (MQ=255)
          aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaa                                                                                                                                        <  1:238393/149‑2 (MQ=255)
          aaaGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCATGCTGGC‑aaaaa                                                                                                                                        <  3:264773/149‑2 (MQ=255)
             gTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑AAAAaccca                                                                                                                                    <  4:1098635/149‑3 (MQ=255)
              ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGc                                                                                                                                    >  1:1367679/1‑149 (MQ=255)
              ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGc                                                                                                                                    <  1:1097886/149‑1 (MQ=255)
              ttGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGc                                                                                                                                    >  3:307952/1‑149 (MQ=255)
                  tCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGG‑TTTTGCCATGCTGGC‑AAAAaccca                                                                                                                                    <  1:1331260/144‑3 (MQ=255)
                   cAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGaaa                                                                                                                               <  2:1300603/149‑1 (MQ=255)
                     aCGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAAccc                                                                                                                                      <  1:1363504/140‑1 (MQ=255)
                      cGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAc                                                                                                                            <  1:869565/149‑1 (MQ=255)
                       gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAcc                                                                                                                           >  4:748098/1‑149 (MQ=255)
                       gTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGCCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCTGAAAAAcc                                                                                                                           >  2:1208479/1‑149 (MQ=255)
                         ctcccGCAACCTGCCGAACGTATCGGATATAGCGCCACCAGCCGACCGCGCGAAGACACCGACAAATGGGCCAAAGAACGCCGGGCGCAGAATATTCACATCCGGCAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAAACGa                                                                                                                         <  2:948708/147‑1 (MQ=255)
                         cACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGCCCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCAAGCATGGC‑AAAACCCGCAGAAAAACCgc                                                                                                                         <  2:689654/149‑2 (MQ=255)
                         cACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGa                                                                                                                         <  3:897587/149‑1 (MQ=255)
                            ccGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATaa                                                                                                                      >  2:747261/1‑149 (MQ=255)
                                 cgccgcCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGaa                                                                                                                 <  4:610515/149‑1 (MQ=255)
                                  gccgccGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAAc                                                                                                                <  4:1024176/149‑1 (MQ=255)
                                        gAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAgg                                                                                                          <  1:1192196/149‑1 (MQ=255)
                                         aaCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGt                                                                                                         <  3:151678/149‑1 (MQ=255)
                                          aCTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTg                                                                                                        >  1:222984/1‑149 (MQ=255)
                                              aTCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCaa                                                                                                    >  1:187805/1‑149 (MQ=255)
                                               tCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                   >  1:784256/1‑149 (MQ=255)
                                               tCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAg                                                                                                   >  1:853960/1‑149 (MQ=255)
                                                cGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAgg                                                                                                  >  1:63430/1‑149 (MQ=255)
                                                  gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                >  4:739490/1‑149 (MQ=255)
                                                  gAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTa                                                                                                >  4:1165860/1‑149 (MQ=255)
                                                     aTAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAg                                                                                             <  4:1328553/149‑1 (MQ=255)
                                                       aGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCa                                                                                           <  4:1174888/149‑1 (MQ=255)
                                                        gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGAt                                                                                                                        <  1:360499/119‑1 (MQ=255)
                                                        gCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAg                                                                                          >  4:809994/1‑149 (MQ=255)
                                                         caccacCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAgg                                                                                         >  2:1123757/1‑149 (MQ=255)
                                                          accaccAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGc                                                                                        >  1:1222846/1‑149 (MQ=255)
                                                           ccaccaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCt                                                                                       >  2:685965/1‑149 (MQ=255)
                                                             accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                     <  2:486329/149‑1 (MQ=255)
                                                             accaACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTca                                                                                     >  1:493805/1‑149 (MQ=255)
                                                                 accgaccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTcagcag                                                                                 >  4:59080/1‑149 (MQ=255)
                                                                     accgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCag                                                                             <  3:174268/149‑1 (MQ=255)
                                                                      ccgCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCaga                                                                            <  2:721663/149‑1 (MQ=255)
                                                                        gcgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAgaga                                                                          <  2:437057/149‑1 (MQ=255)
                                                                         cgcgATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGAt                                                                         <  3:410264/149‑1 (MQ=255)
                                                                            gATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCa                                                                      <  1:212087/149‑1 (MQ=255)
                                                                                 cACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGt                                                                 <  2:265173/149‑1 (MQ=255)
                                                                                  aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                >  1:202359/1‑149 (MQ=255)
                                                                                  aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                <  3:670137/149‑1 (MQ=255)
                                                                                  aCCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGtt                                                                >  4:1073518/1‑149 (MQ=255)
                                                                                   ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                               >  2:656768/1‑149 (MQ=255)
                                                                                   ccGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTg                                                               <  3:172363/149‑1 (MQ=255)
                                                                                    cGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCagag                                                                           <  2:860495/136‑1 (MQ=255)
                                                                                     gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATaaa                                                                                                                     <  1:1157276/93‑1 (MQ=255)
                                                                                     gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATaaa                                                                                                                     >  2:1157276/1‑93 (MQ=255)
                                                                                     gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                             >  2:810483/1‑149 (MQ=255)
                                                                                     gATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTa                                                             >  4:339388/1‑149 (MQ=255)
                                                                                      aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            <  2:222984/149‑1 (MQ=255)
                                                                                      aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            <  2:1367679/149‑1 (MQ=255)
                                                                                      aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            <  2:659338/149‑1 (MQ=255)
                                                                                      aTAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            <  3:739490/149‑1 (MQ=255)
                                                                                              gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTg                                                                                                        <  2:205969/97‑1 (MQ=255)
                                                                                              gCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGa                                                    <  3:1073518/149‑1 (MQ=255)
                                                                                                cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                  <  3:284776/149‑1 (MQ=255)
                                                                                                cAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGc                                                  >  3:145491/1‑149 (MQ=255)
                                                                                                 aaaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCt                                                 >  3:1415113/1‑149 (MQ=255)
                                                                                                  aaGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAgg                                                      <  2:782507/143‑1 (MQ=255)
                                                                                                   aGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTgg                                               >  1:1381766/1‑149 (MQ=255)
                                                                                                     aaCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTc                                             <  1:1271434/149‑1 (MQ=255)
                                                                                                      aCGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAGCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCg                                            >  1:1066545/1‑149 (MQ=255)
                                                                                                        gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            <  4:1051734/131‑1 (MQ=255)
                                                                                                        gCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTaa                                                            >  3:1051734/1‑131 (MQ=255)
                                                                                                          cAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                        >  1:479114/1‑149 (MQ=255)
                                                                                                              cgcAGAATATTGACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAATCCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCaa                                        >  1:1207054/1‑145 (MQ=255)
                                                                                                              cgcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGa                                    <  2:936337/149‑1 (MQ=255)
                                                                                                               gcAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAg                                   <  4:196140/149‑1 (MQ=255)
                                                                                                                  gAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGgc                                 >  4:1141462/1‑148 (MQ=255)
                                                                                                                     tatTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGt                             <  1:1456042/149‑1 (MQ=255)
                                                                                                                           caTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTg                       <  1:576238/149‑1 (MQ=255)
                                                                                                                            aTCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTgg                      >  4:13006/1‑149 (MQ=255)
                                                                                                                               cGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGc                     >  2:120103/1‑147 (MQ=255)
                                                                                                                                     cTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcat             <  3:718067/149‑1 (MQ=255)
                                                                                                                                        gggttttGCCATGCTGGCAAAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAg                                                   >  4:151807/4‑109 (MQ=255)
                                                                                                                                        ggTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGC‑‑GCGTTGTAAGACAGGGAGc                                                  <  4:433684/107‑1 (MQ=255)
                                                                                                                                        ggTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGC‑‑GCGTTGTAAGACAGGGAGc                                                  >  3:433684/1‑107 (MQ=255)
                                                                                                                                         gTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATAtcattcat         >  1:147105/1‑149 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc        <  4:145491/149‑1 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc        <  3:1142828/149‑1 (MQ=255)
                                                                                                                                          tttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATc        <  2:124416/149‑1 (MQ=255)
                                                                                                                                           ttttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATcc       <  4:785578/149‑1 (MQ=255)
                                                                                                                                            tttGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATccc      <  1:517924/149‑1 (MQ=255)
                                                                                                                                                ccAGCATGGC‑AAAACCCGCAGAAAAACCGATAAGCGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAc  >  3:1242021/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCATAAAAATAAAGTTGATCAACGTCACCCGCACGCCGCCGAACTTATCGGAAATAGCACCACCAACCGACCGCGCGATGGCACCGATAAATGGGCCAAAGAACGCCAGGCGCAGAATATTCACATCCGGGAACTGGGTTTTTGCCAGCATGGC‑AAAACCCGCAGAAAAACCGATAAACGAACCGAAGGTGGCAAGGTAAAGCAGGCTCAGCAGCCAGAGATGCAGGCGTTGTAAGACAGGGAGCTGGTCGGCAATTGAGGCGCGTGAACTGGCGATATCATTCATCCCTGAC  >  NC_000913/1543082‑1543374

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑