Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,793 A→G 9.6% 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 9.6% 151.7 / 4.0 61intergenic (+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/39);  new base G (3/3);  total (19/42)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTG  >  NC_000913/1286646‑1286941
                                                                                                                                                   |                                                                                                                                                    
tAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAGCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAg                                                                                                                                                    <  2:1145223/150‑1 (MQ=255)
  aaCGTGGTGCCATCGGGCAATATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTc                                                                                                                                                  <  1:718049/150‑1 (MQ=255)
   aCGTGGTGCCATCGGGCAAAATTAGCTCGGAGTATCAGGTTTTGATTATTAAAGGAAGCGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCAAGTCg                                                                                                                                                 <  1:317946/150‑1 (MQ=255)
     gtggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                                               >  1:1928538/1‑150 (MQ=255)
     gtggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGaa                                                                                                                                               <  2:1251040/150‑1 (MQ=255)
      tggtgCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCACCGCTCTCACGCTGGCGCACAAGTCGAAc                                                                                                                                              >  2:1772946/1‑150 (MQ=255)
         ggCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCtt                                                                                                                                           <  1:1329501/149‑1 (MQ=255)
           ccATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAACGCACACTCGCAGGCCCTCAACACGAAccaatt                                                                                                                                         >  1:1791934/1‑146 (MQ=21)
                 ggCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAg                                                                                                                                   <  1:261182/150‑1 (MQ=255)
                  ggaaaaTTAGCTCGGAGTTACAGGTTTTGATTATTTAAGGATGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGc                                                                                                                                  <  2:1587639/148‑1 (MQ=255)
                   cAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGGCGAACCTTTGTCGAAGCt                                                                                                                                 >  1:1921941/1‑150 (MQ=255)
                   cAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCCGCAACGCACTCTCGCAGCCGCTCAAGACGAACCTTAGTCGAAGCt                                                                                                                                 >  1:1913986/1‑150 (MQ=18)
                      aaTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtct                                                                                                                              <  2:125263/150‑1 (MQ=255)
                       aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                             <  2:550759/150‑1 (MQ=255)
                       aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCAtctc                                                                                                                             >  1:1458412/1‑150 (MQ=255)
                       aTTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTAGTCGAAGCTtctc                                                                                                                             >  1:434989/1‑150 (MQ=255)
                        ttAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                                            <  1:643397/150‑1 (MQ=255)
                           gCTCGGAGTAACCGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                         <  1:1298818/150‑1 (MQ=255)
                           gCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                                                         <  2:455260/150‑1 (MQ=255)
                            cgcGGAGTCACAGGTTTTGATTATTTAAGGACTCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCt                                                                                                                        <  2:951188/148‑1 (MQ=255)
                               ggAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                     <  1:676405/150‑1 (MQ=255)
                               ggAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                     <  2:1100162/150‑1 (MQ=255)
                                gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                    <  1:388393/150‑1 (MQ=255)
                                gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                    <  2:1188774/150‑1 (MQ=255)
                                    aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTGGTCGCATCATCTAATACTTCCCCACt                                                                                                                >  2:1507732/1‑150 (MQ=17)
                                    aaCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAACTTGGTCGAAGCTTCTAATCCTTcccccca                                                                                                                >  2:129083/1‑147 (MQ=255)
                                     aCAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGcat                                                                                                               <  2:1438649/150‑3 (MQ=255)
                                         gTTTTGAGTATTTACGTAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                           <  2:798172/150‑1 (MQ=11)
                                          ttttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTAAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTGCCCGCGTGGGCa                                                                                                          >  2:666501/1‑150 (MQ=255)
                                             tGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGaa                                                                                                       <  1:1208730/150‑1 (MQ=255)
                                               attattTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTTCGGCAACGCTCTCTCGCTGTCGCTCAAGACGAAACTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                                     >  2:1883162/1‑150 (MQ=11)
                                                  attTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAttt                                                                                                  <  1:1604172/150‑1 (MQ=255)
                                                   tttAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTg                                                                                                 <  1:1771254/150‑1 (MQ=255)
                                                     tAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTCGTCGATGCGTCTCATCCTTCCCCGCTTGGGCAGAATATTTGAt                                                                                               >  2:110204/1‑150 (MQ=11)
                                                         gcaGCGCTGGTGGTGGGGGAAGGCTTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCg                                                                                           <  1:152027/148‑1 (MQ=255)
                                                           cgCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGa                                                                                         <  2:2006230/149‑1 (MQ=38)
                                                            gggttggtggtggGGGAAGGATTACTCTGCGCTGCGCGCTTGGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTAGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATTTTTGATTGCGGAt                                                                                        <  1:491341/146‑1 (MQ=11)
                                                            gCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGAt                                                                                        <  2:967226/150‑1 (MQ=38)
                                                                   gggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTg                                                                                 <  2:232446/149‑1 (MQ=38)
                                                                      ggggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCATGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTgaga                                                                              <  2:1459147/149‑1 (MQ=255)
                                                                              ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCATTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCggg                                                                      <  1:741806/150‑1 (MQ=255)
                                                                               gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCgggg                                                                     <  1:1049221/150‑1 (MQ=255)
                                                                               gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGATTTCAgggg                                                                     >  2:1986456/1‑150 (MQ=11)
                                                                                      cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTGAAGTCGATCCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGa                                                              <  1:1254642/150‑1 (MQ=255)
                                                                                       cgcgcTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCAAGTCGAACCTGGGTCGAAGCTTCTCATCCTTCCCCGCGTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGaa                                                             <  1:1217029/149‑1 (MQ=255)
                                                                                         cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAg                                                           <  2:1297045/150‑1 (MQ=255)
                                                                                          gcTGCGCGCTTCGCCCTTCGGGTCGATGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGCACCTTGGGCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAAGATTTGAATGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGt                                                          <  1:981642/150‑1 (MQ=17)
                                                                                             gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAT                                                       <  2:1458412‑M1/150‑1 (MQ=255)
                                                                                               gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGG                                                     <  2:1825409‑M1/150‑1 (MQ=255)
                                                                                                 gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGACGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAATTAAGGGACTTTTGAAAGTGATGGTG                                                   >  2:1207304‑M1/1‑150 (MQ=255)
                                                                                                     cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGG                                               <  2:509077‑M1/150‑1 (MQ=255)
                                                                                                     cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTGGGGCAGAATATTTGATTGCGGATGCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGG                                               <  1:1207320‑M1/150‑1 (MQ=255)
                                                                                                       cccTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCGGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGG                                             <  2:93532‑M1/150‑1 (MQ=255)
                                                                                                                       cTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGCTTACTAAgcg                             >  1:730107/1‑150 (MQ=14)
                                                                                                                         gCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAATTCAGGGACTTTTGAAAGTGATGGTGGTGTGGGAAGGATTACTCAGCGCt                           >  2:345150/1‑150 (MQ=11)
                                                                                                                                     tctcGCTGGCGCTCAAGTCGACCCTTGGTCGAAGCTTCTCCTCCTTCCCCGCTGGGGCAGACTCTTTGATTGCGGATGCGTTTGCGCCTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc               <  1:1525485/150‑1 (MQ=17)
                                                                                                                                     tctcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGATAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTccaaccatctaagtcgatgacc               >  2:1664996‑M1/1‑128 (MQ=255)
                                                                                                                                         gCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                               >  1:231671/1‑82 (MQ=255)
                                                                                                                                         gCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                                               <  2:231671/82‑1 (MQ=255)
                                                                                                                                            ggCGCTCAAGTCGAACCTTGGGCGAAGCTTCTCCTCCTTCCCCGCTTGGGCAGAATTTTTGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg        <  1:1810427/150‑1 (MQ=255)
                                                                                                                                                  cAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg  >  2:1929143/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TAAACGTGGTGCCATCGGGCAAAATTAGCTCGGAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTG  >  NC_000913/1286646‑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: