Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,047 C→T 5.7% intergenic (+521/+19) 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,287,0470CT 5.7% 213.6 / 2.3 72intergenic (+521/+19)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base C (35/32);  new base T (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAG  >  NC_000913/1286900‑1287194
                                                                                                                                                   |                                                                                                                                                   
aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                                                                    <  3:105013/149‑1 (MQ=255)
       aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                                                                   >  1:269599/1‑143 (MQ=255)
       aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                                                                   <  2:269599/143‑1 (MQ=255)
         tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                           >  8:387386/1‑149 (MQ=255)
           cgcgcTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                           <  5:106920/146‑1 (MQ=255)
             cgatGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGTCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCATCCTTCCCCGCATGTGCAGAATATTTGATTGCGGATTCTCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                       <  5:386409/146‑1 (MQ=17)
                   gcgcTTCGCCCTTCGGGTCGTTGCCTGCGCCAACGCTCTCTCTCTGAATCTATCTCCGTACGTTAGTCGAACCTTCTCATCCTTAAGAGCATGGGCAGAATATTTTATTTCTGATTAGCTTGAGATTTCAGGGTCTTTTGAAAGTGAtg                                                                                                                                 <  2:224839/149‑1 (MQ=9)
                     gcTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTATCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                                                                                                                               >  8:39830/1‑149 (MQ=255)
                      cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtg                                                                                                                              >  4:165540/1‑149 (MQ=255)
                      cTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtg                                                                                                                              >  2:48633/1‑149 (MQ=255)
                          gcgcTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGGAGTATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                                                                                                                          <  7:387386/146‑1 (MQ=255)
                                   gTCGATGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTATCATCCTTCCCCGCATGGGGAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                                                                                                 >  2:109218/1‑149 (MQ=255)
                                    tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                       <  4:311171/126‑1 (MQ=255)
                                    tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                       >  3:311171/1‑126 (MQ=255)
                                       ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                                                             <  1:48633/149‑1 (MQ=255)
                                              cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTTAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                                                                                                      >  5:229668/1‑149 (MQ=255)
                                               ggCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCACCGCATGGGCAGAATAATTGATTGCGGATTCGCGTGAGAGTTCAGTGAATTTTGAAAGTGATGGTGGTGGGGGATGGATTACTCAGCGCTg                                                                                                     >  4:314803/1‑149 (MQ=17)
                                                 cAACGCTCTCTCGCTGACACTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCAGCATGGGCAGAAAATTTGATTGGGGATTCTCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                                                                                                   >  4:356566/1‑149 (MQ=17)
                                                     gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                                                                               >  3:351170/1‑149 (MQ=255)
                                                     gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACGTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt                                                                                               >  4:216692/1‑149 (MQ=255)
                                                         tctcGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTACCCCGCATGGGCAGAAAATTTGATTGCGCAGTCGCTTGAGAGTGCAGCGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc                                                                                           >  1:183803/1‑149 (MQ=17)
                                                         tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                                                             <  3:165540/147‑1 (MQ=255)
                                                             gCTGACGCTCGAGCCGCGCCTTAGTCGAGGCTTCTCATCCTTCCCCGCAGGGGCAGAATATTTGATTGCGGATTGGCGTGAGGGGTGTGGGGCTTTGGATAGTGGGGGGGGGGGGGGATGGATTACTCCGCGCTGGGCGCTTCGccggt                                                                                       >  4:50805/1‑146 (MQ=9)
                                                                   acTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTc                                                                                 >  6:82735/2‑149 (MQ=255)
                                                                    cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                                <  3:216692/149‑1 (MQ=255)
                                                                       gAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                                                              <  3:50805/148‑1 (MQ=255)
                                                                        aGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGATTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGc                                                                            <  2:183803/149‑1 (MQ=255)
                                                                        aGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGc                                                                            <  7:39830/149‑1 (MQ=255)
                                                                         gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGcc                                                                           >  3:102444/1‑149 (MQ=255)
                                                                           cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTTAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                         <  6:229668/149‑1 (MQ=255)
                                                                             aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                       <  4:351170/149‑1 (MQ=255)
                                                                                   aGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAc                                                                  <  4:102444/148‑1 (MQ=255)
                                                                                    gTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                                                                                                             >  4:261952/1‑72 (MQ=255)
                                                                                    gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                                                                                                             <  3:261952/72‑1 (MQ=255)
                                                                                     tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                  <  7:227953/98‑1 (MQ=255)
                                                                                     tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                                                                                                  >  8:227953/1‑98 (MQ=255)
                                                                                      cGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                                                                     >  6:4446/1‑94 (MQ=255)
                                                                                      cGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                                                                                                     <  5:4446/94‑1 (MQ=255)
                                                                                       gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                                                                  >  8:76477/1‑64 (MQ=255)
                                                                                       gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTAAAAGCTGAGAACATCACCCTACGGGTCGTTGCATGTGGctag                                                                  >  8:175881/1‑141 (MQ=14)
                                                                                       gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                                                                  <  7:76477/64‑1 (MQ=255)
                                                                                       gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                                                                                                          >  6:274254/1‑104 (MQ=255)
                                                                                       gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAg                                                                                                          <  5:274254/104‑1 (MQ=255)
                                                                                           cTTCTCATCCTTCCCCGCAGGGGAAGAATTTTTGATTGCGGATTCGCTTGAGAGTTTAGGGACTTTTGAAAGTGATGGTGGTGGTGGAATGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                         >  4:132774/1‑149 (MQ=14)
                                                                                                 aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                                                                                                            >  4:288914/1‑60 (MQ=255)
                                                                                                 aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACttt                                                                                                                                            <  3:288914/60‑1 (MQ=255)
                                                                                                 aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                                                     >  6:294026/1‑147 (MQ=255)
                                                                                                 aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                                                   >  6:203902/1‑149 (MQ=255)
                                                                                                 aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                                                     <  5:294026/147‑1 (MQ=255)
                                                                                                    cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                <  5:203902/149‑1 (MQ=255)
                                                                                                       ccccGCATGGGGAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                              <  1:109218/148‑1 (MQ=255)
                                                                                                       ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                              <  5:82735/148‑1 (MQ=255)
                                                                                                        cccGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg                                                                                   <  4:300819/110‑1 (MQ=255)
                                                                                                        cccACATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                            >  7:200008/1‑149 (MQ=255)
                                                                                                         ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                 >  7:293634/1‑143 (MQ=255)
                                                                                                                 gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                            >  6:242010/1‑124 (MQ=40)
                                                                                                                 gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                            <  5:242010/124‑1 (MQ=40)
                                                                                                                 gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                   <  3:314803/149‑1 (MQ=255)
                                                                                                                  cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAGTTTTGAAAGGGATGGTGGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                  <  3:356566/149‑1 (MQ=255)
                                                                                                                  cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                  >  5:196507/1‑149 (MQ=255)
                                                                                                                        atTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTTGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                            >  2:61578/1‑149 (MQ=255)
                                                                                                                            gATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                        >  6:120855/1‑149 (MQ=255)
                                                                                                                                gCGGGTGCGCTTCAGATGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                     <  1:372245/148‑1 (MQ=11)
                                                                                                                                    aTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                                                              >  2:275787/1‑87 (MQ=255)
                                                                                                                                     ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc               >  8:43828/1‑149 (MQ=255)
                                                                                                                                     ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATCCTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGCACCTTGGTCGACGCTTCTCCTCCTTc               >  1:350560/1‑149 (MQ=14)
                                                                                                                                       cGCTTGAGAGTTTAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc              <  3:132774/148‑1 (MQ=255)
                                                                                                                                        gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcccc            >  6:65965/1‑149 (MQ=255)
                                                                                                                                        gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc             <  1:61578/148‑1 (MQ=255)
                                                                                                                                        gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCCGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcccc            >  2:157445/1‑149 (MQ=255)
                                                                                                                                              gagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg      >  8:38836/1‑149 (MQ=255)
                                                                                                                                                 ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa   <  5:65965/149‑1 (MQ=255)
                                                                                                                                                   cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg  >  1:108489/1‑148 (MQ=40)
                                                                                                                                                   |                                                                                                                                                   
AAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAG  >  NC_000913/1286900‑1287194

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: