Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,793 A→G 8.2% 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.2% 140.2 / 3.0 60intergenic (+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 (16/38);  new base G (2/3);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTT  >  NC_000913/1286645‑1286940
                                                                                                                                                    |                                                                                                                                                   
gTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCa                                                                                                                                                     >  1:529743/1‑149 (MQ=255)
  acaCGTGGTGCCATCGGGCAAAATTAGCTGGGAGTAACAGGTTTTGATTATTTCAGGAAGCGATGGTGGTGGGGGAAGCATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAg                                                                                                                                                   <  1:1022630/147‑1 (MQ=255)
       gggtgCCATCGGGCAAAATTAGCGCGGAGTACCAGGTTTTGATTATTTAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGaa                                                                                                                                              <  1:196339/148‑1 (MQ=255)
             cctCGGGCAAAATTAGTTTGGCGTCCCCGGTTTTGCTTATGTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTgg                                                                                                                                        <  2:1571286/147‑1 (MQ=17)
                   gCAAAATTAGCTCGGAGTAACAGGTTTTGATTCTTTAAGGCAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCATGGTCGAAg                                                                                                                                  <  2:939747/149‑1 (MQ=255)
                     aaaaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATCGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCCGCAACGCTCTCTCCCTGGCGCACAAGTCGAACATTCGTCGAAGca                                                                                                                                >  1:1214942/1‑148 (MQ=18)
                         ttAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                            <  2:1046075/149‑1 (MQ=255)
                         ttAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                            <  1:941382/149‑1 (MQ=255)
                          tAGCTCGGAGTAACAGGTTTTGATTCTTTACGGCCGCGATGGTGGTGGGGGAAGCATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCa                                                                                                                           <  2:94534/149‑1 (MQ=17)
                             cTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCCTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTAGTCGAAGCTTCTCATcc                                                                                                                        >  2:168362/1‑149 (MQ=255)
                              tCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                            <  2:903827/144‑1 (MQ=255)
                              tCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                            >  1:903827/1‑144 (MQ=255)
                                gtaGTAACAGCTTTTGATTTTTTAAGGAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                                     <  2:529743/147‑1 (MQ=11)
                                 gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                    <  2:155620/149‑1 (MQ=255)
                                 gAGTAACAGGTTTGGCTTATTTAAGGACGCGATGGTGGTGGGGGAAGGATTACTGAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCGTGCGGCAACGCGCTCTCGCAGGCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCATCCTTcc                                                                                                                    <  1:176679/149‑1 (MQ=17)
                                  aGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                   >  2:46886/1‑149 (MQ=255)
                                      aCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTCCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCt                                                                                                               <  2:968980/149‑1 (MQ=255)
                                        atgTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCtgg                                                                                                             <  2:899126/147‑3 (MQ=255)
                                         ggTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTgg                                                                                                            <  2:594259/149‑1 (MQ=255)
                                            tttGATTATTTACGGCAGAGATGGTGGTGTGGGAAGGATAACTCAGCGCTGGGGGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGATCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGTAGCTTCTCATCCTTCCCCGCTTGTGCa                                                                                                         <  2:331791/149‑1 (MQ=14)
                                            tttGATTATTTAAGGAAGCGCTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCa                                                                                                         <  1:1017133/149‑1 (MQ=255)
                                             ttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAg                                                                                                        <  1:1118690/149‑1 (MQ=255)
                                                    tttAAGGAAGCGCTGGGGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAttt                                                                                                 <  1:1056341/149‑1 (MQ=255)
                                                     ttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTTACGCTCAAGACAACCCTTAGTCGAAGCTACTCATCCAACCCCCCAAGGGAAGAAAAtttt                                                                                                >  2:338027/1‑148 (MQ=11)
                                                      tAAGGCAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGa                                                                                               <  1:653928/149‑1 (MQ=255)
                                                         ggAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                           <  1:1244296/118‑1 (MQ=14)
                                                         ggAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                           >  2:1244296/1‑118 (MQ=14)
                                                          taaGAGATGGTGGTGGGGGATGGATTAGTCAGCGCTGCGCGCTTCGCCCTTCTGGTCGTTGCATGCGGCAACGCTCTCTCGCTGTCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGc                                                                                           <  2:896916/148‑1 (MQ=11)
                                                          gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGc                                                                                           <  1:989057/149‑1 (MQ=40)
                                                          gAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGc                                                                                           >  2:1121337/1‑149 (MQ=40)
                                                           aaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                                          <  1:1128171/149‑1 (MQ=39)
                                                            aGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCgg                                                                                         <  1:1277131/149‑1 (MQ=255)
                                                                ttggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCGTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTc                                                                                     <  2:1475229/148‑1 (MQ=255)
                                                                   ggtgtggGGGATGGAGTACTGAGTGCTGCTCGCTTCGGCCTTCGGGTCGTTGCTGGCGGAAACTCTCTATATCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGtt                                                                                  <  1:1010953/146‑1 (MQ=11)
                                                                    tgggggGGGAAGGATGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCGTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGttt                                                                                 <  1:908482/149‑1 (MQ=255)
                                                                     ggggggggTCGGATTACTCATCGCTGCGCGCGTCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTGTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTg                                                                                <  1:194593/146‑1 (MQ=11)
                                                                      gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTga                                                                               <  1:46886/149‑1 (MQ=39)
                                                                      gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTga                                                                               <  2:1359112/149‑1 (MQ=39)
                                                                      gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTga                                                                               <  2:1049831/149‑1 (MQ=39)
                                                                                     cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCttt                                                                >  2:949349/1‑149 (MQ=255)
                                                                                      ctaGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCtttt                                                               <  1:656247/147‑1 (MQ=255)
                                                                                        aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATGATTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTACGGGGCTTTTGa                                                             >  2:496971/1‑149 (MQ=255)
                                                                                         gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAAACTTGGTCGAAGCTTCTCATCCTTACCCGATTTGGCAGAATATTTGAATGCGGATTCGTTTGAGAATTCCGCGGCATTTGat                                                            >  2:656247/1‑148 (MQ=17)
                                                                                           gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAg                                                          <  1:1380419/149‑1 (MQ=255)
                                                                                              gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGa                                                       <  1:1121337/149‑1 (MQ=255)
                                                                                               cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAt                                                      <  1:686497/149‑1 (MQ=255)
                                                                                                gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCAGGGGCTTTTGAAAGTGAtg                                                     >  2:1254080/1‑149 (MQ=255)
                                                                                                gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTGGCGCTCAACTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTGGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtg                                                     <  1:989451/149‑1 (MQ=255)
                                                                                                 cgcTTCGCGCTGCGGGGCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGGCGATCCTTGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGTTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtgg                                                    <  2:1214942/149‑1 (MQ=14)
                                                                                                 cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAAGGCTATCTCGCTGGCGCTCAAGTCGATCCTTGGTCGAAGCTTCTCATCCTTGCGCGCTTGGGCAGTATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtgg                                                    <  1:369770/149‑1 (MQ=17)
                                                                                                 cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAtgg                                                    <  2:85922/149‑1 (MQ=255)
                                                                                                    ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGGTAGTTCCGGGACTTTTGAAAGTGAtggtgg                                                 >  2:501260/1‑149 (MQ=255)
                                                                                                        cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                             >  1:809237/1‑149 (MQ=255)
                                                                                                        cccTTCGGGGCGTTGCCTGCGGAAACGCGCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTggg                                             <  1:302420/149‑1 (MQ=255)
                                                                                                               ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                      >  1:1094824/1‑149 (MQ=255)
                                                                                                                 tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAATTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGAtt                                    >  2:887959/1‑149 (MQ=11)
                                                                                                                    ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAca                                 >  2:201813/1‑148 (MQ=255)
                                                                                                                                 cGCTCTCTCGCTGGCGCTCAAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                     <  1:88813/84‑1 (MQ=11)
                                                                                                                                 cGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                     >  2:88813/1‑84 (MQ=18)
                                                                                                                                      tctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGc               >  1:512365/1‑149 (MQ=255)
                                                                                                                                                   cAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt  >  2:402876/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTT  >  NC_000913/1286645‑1286940

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: