Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,793 A→G 8.8% intergenic (+267/+273) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,7930AG 8.8% 162.8 / 4.3 69intergenic (+267/+273)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (13/46);  new base G (4/2);  total (18/51)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTG  >  NC_000913/1286647‑1286941
                                                                                                                                                  |                                                                                                                                                    
aaaCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGTTGCCTGAGGCAACGCTATCTCGCTGGCGCTCAAgg                                                                                                                                                   >  1:28290/1‑149 (MQ=255)
  aCGTGGTGCCATCGGGCAAAATTAGCTCGTAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCg                                                                                                                                                 <  2:677253/150‑1 (MQ=255)
     tggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAc                                                                                                                                              <  1:1443305/150‑1 (MQ=255)
        tgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCtt                                                                                                                                           >  2:1135322/1‑150 (MQ=255)
         gCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCttt                                                                                                                                          >  1:512424/1‑149 (MQ=255)
          ccATCGTGCAAAATTAGCTCGGCGTAACAGGTTTTGATTATTTAAGTAAGCGATGGGGGTGGGGGAAGGATTACTCAGCTCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCATgg                                                                                                                                         <  2:1287013/150‑1 (MQ=18)
            aTCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTAGTc                                                                                                                                       >  1:1157619/1‑150 (MQ=255)
            aTCGGGCAAAATTAGCGCGGAGTAACAGGTTTTGATTTTTTAAGGACGCGATGGTGGGGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTc                                                                                                                                       <  2:1377347/150‑1 (MQ=18)
              cGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGa                                                                                                                                     <  2:1513849/150‑1 (MQ=255)
                 gCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                                                  <  1:1964323/150‑1 (MQ=255)
                 ccAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGCAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                                                  >  1:351735/2‑150 (MQ=255)
                  cAAAATTAGCTCGGAGTAAGATGTTTTGTTTATTTATGGACGCTATGGTGGTGGGGTAAGGATTACTCCGGGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGACCCTTGGTCGAAggt                                                                                                                                 <  2:1921729/150‑3 (MQ=17)
                    aaaTTAGCTCGGAGTAACAGGTTTTGATTTTTTAAGGAAGCGATGGTGGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtc                                                                                                                               <  1:355222/150‑1 (MQ=255)
                    aaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtc                                                                                                                               <  2:1314116/150‑1 (MQ=255)
                     aaTTAGCTCGGAGTAACAGGTTTTGATTATTTACCGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtct                                                                                                                              <  2:996832/150‑1 (MQ=255)
                      aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                             <  2:1268632/150‑1 (MQ=255)
                       ttAGCGCGGTGTATCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGCAGGATTCCTCAGCTCGGCGCGCTTCGCCCTGCGGGTCGTGGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctct                                                                                                                            <  1:96778/150‑2 (MQ=17)
                        tAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                                           <  2:1434381/150‑1 (MQ=255)
                          gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCCAGACGAACCTTTGTCGAAGCTTCTCATcc                                                                                                                         >  1:393243/1‑150 (MQ=255)
                          gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCAAGACGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                         >  2:1245388/1‑150 (MQ=255)
                          gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCa                                                                                                                            >  1:70872/1‑147 (MQ=255)
                          gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTAGTCGAAGCTTCTCa                                                                                                                            <  2:70872/147‑1 (MQ=255)
                           cTCGGAGTCACAGGTTTTGATTATTTAAGCACGCGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCt                                                                                                                        <  1:299718/150‑1 (MQ=255)
                           cTCGGAGTAACAGGTTTTGATTATTTACGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCt                                                                                                                        <  2:2096829/150‑1 (MQ=255)
                              ggAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                     <  1:310809/150‑1 (MQ=255)
                                   aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACCCTATCTCGCTGGCGCTCAACTCGCACCATTGGAGCAGCATCACATCCTTCCCCGcc                                                                                                                >  2:1871742/1‑149 (MQ=17)
                                        gTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGGCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATggag                                                                                                           >  2:909184/1‑148 (MQ=255)
                                        gTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCAAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGc                                                                                                           >  2:2093608/1‑150 (MQ=255)
                                             gATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                                            <  2:366988/112‑1 (MQ=255)
                                             aattattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                                            >  1:366988/2‑112 (MQ=255)
                                               ttattTAAGGAAGCGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAatat                                                                                                    <  2:1369199/150‑1 (MQ=255)
                                                tattTAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAtt                                                                                                   <  1:691235/150‑1 (MQ=255)
                                                 ttttCAGGAAGCGATGGTGTTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAttt                                                                                                  <  2:1473815/149‑1 (MQ=255)
                                                 tattAAGGAAGCGATGGTGGTGGGGGAATGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCTCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAttt                                                                                                  <  2:351735/148‑1 (MQ=255)
                                                  tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                                 <  1:425822/150‑1 (MQ=255)
                                                  tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                                 <  2:691555/150‑1 (MQ=255)
                                                  tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                                 <  1:535796/150‑1 (MQ=255)
                                                    tAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                               >  1:123306/1‑150 (MQ=9)
                                                       ggAAGCGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTGCGGTTCGTGGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGc                                                                                            <  1:1032864/150‑1 (MQ=14)
                                                        gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAGTATTTGATTGCg                                                                                           <  2:974406/150‑1 (MQ=255)
                                                        gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                                           <  2:852244/150‑1 (MQ=40)
                                                        gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                                           <  2:844837/150‑1 (MQ=40)
                                                        gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                                           <  1:866110/150‑1 (MQ=255)
                                                         aaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCgg                                                                                          <  2:1888898/150‑1 (MQ=39)
                                                            cGATGGTGGAGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGAtt                                                                                       <  2:114195/150‑1 (MQ=38)
                                                               tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGt                                                                                    <  1:994944/150‑1 (MQ=37)
                                                                ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGtt                                                                                   <  2:1431654/150‑1 (MQ=37)
                                                                  tggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTg                                                                                 <  2:1999502/150‑1 (MQ=38)
                                                                  gggggggggAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCGCGCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTGGGGCAGAATATTTGATTGCGGATTCGTTTg                                                                                 <  2:339066/146‑1 (MQ=11)
                                                                   ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTga                                                                                <  2:1682916/150‑1 (MQ=255)
                                                                    gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCATTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgag                                                                               <  2:1950472/150‑1 (MQ=255)
                                                                     ggggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgaga                                                                              <  2:1659340/149‑1 (MQ=255)
                                                                             ggATTACTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGGGAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                                      <  1:94881/150‑1 (MQ=255)
                                                                                   cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCtttt                                                                <  1:1135322/150‑1 (MQ=255)
                                                                                     cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCATTGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCAGGGGCTTTTGa                                                              >  2:1708598/1‑150 (MQ=11)
                                                                                        cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAg                                                           <  2:142960/150‑1 (MQ=255)
                                                                                          cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                         >  2:405511/1‑86 (MQ=33)
                                                                                             ggcgcTTCGCCCTGCGGGTCGTTGCCTGCGGAATCGCTCTCTCGCGGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCTTCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtg                                                      <  1:1030318/149‑1 (MQ=17)
                                                                                                     gCCGTTCGGGTGGTTTCGTGCTGCAACGTTCTCTCGCTGTCGCTCAAGACGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAGTTCTGTTACTTTTGAAAGTGATGGTGGTggg                                              <  2:28290/150‑1 (MQ=11)
                                                                                                      cccTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCCGGCGCTCCAGTCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATATGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTgggg                                             <  2:243976/150‑1 (MQ=14)
                                                                                                          tCGGGTCGTTGCCTGCGGAAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAg                                         <  2:1533228/150‑1 (MQ=255)
                                                                                                          gcGGGTCGTTGCCTGCGGCATCGGTCGCTCGCGGGCGCTCAAGCCGAACGTTGGTCGAAGGTTCTCATCCTTCCCCGCGTGGGCAGAATATATGATTGCGGATGCGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAg                                         <  2:203766/149‑1 (MQ=14)
                                                                                                           cGGGGCGGTGCCTGCGTCAACGCTCTCGCGCTGGAGCTCCAGGCGACCCGTGTGCTAATCTTCTCCTCCTTCCCCGCGGGGGCAGATTCCTTGATTGCGGAGTCGTATGAGTATTCCGGTTCTTTTGAAAGTGATGGTGGTGGGGGAAgg                                        <  2:245577/150‑1 (MQ=9)
                                                                                                                   tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCTAAGCTTCTCATCCTTCCCCGCTGGGGCAGAATTTTTGATTTCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                <  1:1889779/150‑1 (MQ=255)
                                                                                                                     gcTGAGGCAACGCTCGCGCGCCGGCGCTCCATCCGACCCGTGTTCCAATAGTAGCAGCGTTCCCCGCTTGGGCAGAATCTTATATTCCGTCTTCGTTTCAGTATTCGTGTTCTTAAGAAAGCGATGGTGGTGGGGGAAGGATTACTCAgc                              <  2:1718058/149‑1 (MQ=9)
                                                                                                                                     ctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc              >  1:147959/1‑150 (MQ=255)
                                                                                                                                       cGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCtt            >  1:1999538/1‑150 (MQ=255)
                                                                                                                                             cgcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc      >  1:1664498/1‑150 (MQ=255)
                                                                                                                                                 cAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg  <  2:147959/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTG  >  NC_000913/1286647‑1286941

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: