Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,965 A→G 11.4% intergenic (+439/+101) 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,9650AG11.4% 144.9 / 2.7 79intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/8);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 8.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCT  >  NC_000913/1286819‑1287108
                                                                                                                                                  |                                                                                                                                               
ctCCTTTCCCCTTTGGGCAGATTATTTGATTGGGGCGTCGTTTGAGGTTTCCGGGGCTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                                               <  2:1489324/148‑1 (MQ=14)
       cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                          >  1:673881/1‑147 (MQ=255)
       cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                                          <  2:673881/147‑1 (MQ=255)
       cccGCTGGGGCAGAAGATTTGGTTGCGGATTCGTTTCGGCCGCCGAGTGCTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                        <  1:784997/149‑1 (MQ=14)
        ccGCTTGGGCAGAATATTTGATTGCGGATACGAATCAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGcc                                                                                                                                       >  1:437760/1‑149 (MQ=14)
          gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                     <  1:1167948/149‑1 (MQ=255)
                  cgAATAATTGATTGCCGATTCGTTTGAGGATTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTGCTTCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                             <  1:1053100/148‑1 (MQ=14)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                             >  1:1048752/1‑149 (MQ=255)
                  aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                                             >  1:1048758/1‑149 (MQ=255)
                       gcgtggttGCGGCTGCGTTTGCGGATTCGGGGGTTTTTGGAAGTTATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                        <  2:48957/143‑1 (MQ=14)
                       atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                                        <  1:1514572/149‑1 (MQ=255)
                       atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                        <  1:1551605/149‑1 (MQ=255)
                       atTTGATTGCGGATTCGTTTGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                        <  1:555272/149‑1 (MQ=255)
                                 ggATTCGTTTGAGGATTCCGGGTCTTGTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGTTCGGCTAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                              <  1:435931/149‑1 (MQ=14)
                                 ggAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                                              <  1:602632/149‑1 (MQ=255)
                                  gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACTGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGACGCTTCTCATCCt                                                                                                             >  2:1239347/1‑149 (MQ=255)
                                  gATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                                                                                             >  1:1524907/1‑149 (MQ=255)
                                     tCGTTTGAGAATTCCGGGGCTTTTCAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcc                                                                                                          >  1:915181/1‑149 (MQ=255)
                                       gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTACGCTCTTCGGGTGGTTGCCTGCGGCAGCGCTCTGTAGCTGTCGCTCGAGCCGAACCTTAGTCGACGCTTCTCATCCTTcccc                                                                                                        >  2:256266/1‑149 (MQ=14)
                                       gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGACCCTTCATCGAAGCTCTCCATCCTTcccc                                                                                                        >  2:1440376/1‑149 (MQ=14)
                                        tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCGTCCTTCCCCg                                                                                                       >  1:121684/1‑149 (MQ=255)
                                        tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAACCTGCGAGGCGCTCTTCCCGGTGTTGGTCGTTGCCTGCGGCAACGGTCTATCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCATTCCCCg                                                                                                       >  2:207240/1‑149 (MQ=14)
                                         ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                      >  1:619093/1‑149 (MQ=255)
                                          tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                     <  1:693723/149‑1 (MQ=255)
                                          tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGAATTACTCGGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCCGCGGCCACGCCCTCGCGCGGACGCGCGCGCCGAACCCGGGGCGCAGCTTCTCATCCTACCCCGcc                                                                                                     >  2:1573176/1‑148 (MQ=14)
                                            agaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATg                                                                                                   >  2:781914/1‑149 (MQ=255)
                                                   cGGGGTTTTAGGAAGCGGTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                                            <  2:1797899/149‑1 (MQ=14)
                                                    ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                           <  1:1561171/149‑1 (MQ=255)
                                                     gggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                          >  1:1439952/1‑149 (MQ=39)
                                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGGTTAGCCCTTCGGGTCGTTGCCTGCGGCAAGGGTCTCCCGCTTCCGCTCGAGCCGGGCCGTCGCCGAAGCTTCTCATCCTTCCCCGCATGGGCAGGAGATTTGa                                                                                     >  2:930332/1‑149 (MQ=11)
                                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGGTGACGCTCGAGCCGAACCTTCGTCGTAGCTTCTGATCCTTCGCCGCGTGGGCAGAGTGTTTGa                                                                                     >  2:1423169/1‑149 (MQ=11)
                                                          tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTACGCGCTTCGCCCTTCGGGGCGTCGCCTGCGGCAAGGCTCTCTCGCTGCCGCTCGAGCCGGACCGTAGTCGAAGCTTCTCATCATTCCCCGAATGGGCAGCATATTTGa                                                                                     >  2:930480/1‑149 (MQ=11)
                                                          ttaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCCCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                     <  2:1514182/144‑1 (MQ=255)
                                                            aggatgcgATGGTGGGGGGGGAAGGATGCCTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                   <  2:1211114/142‑1 (MQ=11)
                                                              aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                                 <  1:781914/149‑1 (MQ=18)
                                                                 gcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                              <  1:626102/147‑1 (MQ=255)
                                                                 gcgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                              <  2:24148/147‑1 (MQ=255)
                                                                  tGATGGTGGTGGGGGAAGGATTCCTCCGCGCTGCGCGCTTGGCGGTACGGGAGGGTGGTTACGCGATCGCACACGCACTCACGCTCGAGCCGAACCTTGGTCGACGTTTCACTAACTACCCACCACGGGCGGAAGCTTTGATTGCGGAt                                                                             >  2:830101/1‑149 (MQ=2)
                                                                      ggtggtggGGGAAGGATTCCTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGc                                                                         <  1:508585/149‑1 (MQ=255)
                                                                       gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTGGCCGTTGGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGGACCTTGGTCGAGGGTTCTCATCCTTCCGCGCATGGGCAGAATATTTGATTGCGGAtttgtt                                                                        >  2:178738/1‑145 (MQ=11)
                                                                          gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                     <  2:1048758/149‑1 (MQ=255)
                                                                          gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                     <  2:1048752/149‑1 (MQ=255)
                                                                              ggaaaGGATTACCCAGCGCTGCGCGCTGCGCCCTTCGGGTCATTGCCTGCGGCATCGCTCTCTCGCTGACGCCCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGGATATTTGATTGCGGATTCGCTTGAGAGt                                                                 <  2:1465037/146‑1 (MQ=14)
                                                                               ggAAGGATTACCCAGCGCCGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                <  1:1739487/149‑1 (MQ=255)
                                                                                    gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                           >  1:257978/1‑149 (MQ=255)
                                                                                     aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCTCATGGGCAGTATATTTGATTGCGGATTCGCTGGAGAGTTCAGGGa                                                          >  1:488777/1‑149 (MQ=255)
                                                                                      ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGAAGGTCAGGGAc                                                         >  2:1773337/1‑149 (MQ=255)
                                                                                         cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                      <  1:327336/149‑1 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                     <  1:735237/149‑1 (MQ=255)
                                                                                          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCCCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGCGGTCAGGGACcttt                                                     >  2:1481186/1‑149 (MQ=14)
                                                                                            aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                                           >  1:45970/1‑125 (MQ=25)
                                                                                            aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                                           >  2:1643631/1‑125 (MQ=25)
                                                                                            aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                                           <  2:45970/125‑1 (MQ=25)
                                                                                               gccgcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                <  1:579978/146‑1 (MQ=255)
                                                                                                     cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                     >  1:1459487/1‑90 (MQ=18)
                                                                                                        ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGGCTTTTGACAGTGAtggtgg                                       >  1:1751107/1‑149 (MQ=255)
                                                                                                        ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtgg                                       >  1:1559760/1‑149 (MQ=255)
                                                                                                           gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                    <  2:1261059/149‑1 (MQ=255)
                                                                                                           gCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                    <  2:1669613/149‑1 (MQ=255)
                                                                                                            cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTCGAAAGTGATGGTGGTggg                                   >  1:1815219/1‑149 (MQ=255)
                                                                                                            accgtggggTCCTCGCCTGCACTACCGCACTCTCGCCCTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTGCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                   <  2:1587969/143‑1 (MQ=14)
                                                                                                             ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                  <  1:1692492/149‑1 (MQ=255)
                                                                                                             ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCGGGGACTTTTGAAAGTGATGGTGGTgggg                                  >  2:494349/1‑149 (MQ=255)
                                                                                                               ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTTTTGAAAGTGAAGATggggggggg                                >  2:192724/1‑148 (MQ=255)
                                                                                                                  gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                             <  1:609980/149‑1 (MQ=255)
                                                                                                                      cGTTGCCTGCGTCTTCCGTCCTTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                         <  2:1523877/149‑1 (MQ=14)
                                                                                                                      cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGCGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGATATTGATGGTGGTGGGGGAAGGATTa                         <  1:1010217/149‑1 (MQ=255)
                                                                                                                      cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                         <  2:1524907/149‑1 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                       <  1:1091794/149‑1 (MQ=255)
                                                                                                                        ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACt                       >  2:902998/1‑149 (MQ=255)
                                                                                                                          gCCCGCGGCAGCGCTCTCTCGCTGACGCTCCAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCACAGCATGGGCAGACTCTTTTATTTCGTATTCGCTTGAGAGTTCCGCGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                     <  1:506296/149‑1 (MQ=11)
                                                                                                                                cgccaCGCTCTCTCGCTGGCTCTCGAGCCGAACCTTGGCCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGAATCGCATGAAATTTCGGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTg               <  1:802415/145‑1 (MQ=11)
                                                                                                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGc                                                                                 >  1:165402/1‑80 (MQ=255)
                                                                                                                                    aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCtaagcg           >  2:276912/1‑144 (MQ=255)
                                                                                                                                      gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGAAGGTGAGGGTGGTGGGGGCAGGATTACTCAGCGCTtcgcgcg         >  1:3884/1‑148 (MQ=17)
                                                                                                                                        tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTCCTCATCCTTCCCCGCATGGCCCGAATCCTTGATTGCGGATTCGCTGGAACGCTCCGGGCCTTTTGAAAGGGATGGCGGTGGGGGCAGGATTAGTCCCCTCTGCGCGCTTc       >  2:52606/1‑149 (MQ=255)
                                                                                                                                           ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGACTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc    <  1:790490/149‑1 (MQ=255)
                                                                                                                                            tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGccc   >  1:828812/1‑149 (MQ=255)
                                                                                                                                             cGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCt  >  1:917914/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCT  >  NC_000913/1286819‑1287108

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