Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,163,577 T→C 6.3% intergenic (‑308/‑61) trmA ← / → btuB tRNA m(5)U54 methyltransferase, SAM‑dependent, tmRNA m(5)U341 methyltransferase/vitamin B12/cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,163,5770TC 6.3% 191.5 / 2.1 71intergenic (‑308/‑61)trmA/btuBtRNA m(5)U54 methyltransferase, SAM‑dependent, tmRNA m(5)U341 methyltransferase/vitamin B12/cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base C (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACC  >  NC_000913/4163355‑4163787
                                                                                                                                                                                                                              |                                                                                                                                                                                                                  
gttcGTCTGGGCATTAACTTGTGTGACTCGGTAGTATCCAGTTGCCGGCCGGGGGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCGCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCAACTGGTTCtcatca                                                                                                                                                                                                               <  1:688340/225‑1 (MQ=255)
cgaCGTCTGGCAATTAATTTGGTACACCCTGTAGCTTCCCGGTGCCGGGCCTGTGAGTTAATAGGGAATCCAGTGGGAATCTGGGGCTGACGGGCCGCGGTAAGGAAAGGTGCGATGATTTGGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCtcatca                                                                                                                                                                                                               <  1:503396/226‑1 (MQ=255)
 agcGTTTGGAAAGTAACTTGTTACAACCTGTAGCCTCCACTTGCCGGTCCTGTGAGTTAATGGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCtcatca                                                                                                                                                                                                               <  2:416303/225‑1 (MQ=255)
   cGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCATGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTACGGATAGGTGCGATGATTGCGTTATGCGGCCACTGCCGTTCGGTGCGAAGTCATCATCTGTTAGTATCTTAGCTACCCCTCCAAGCCCGAAGACCTGCCGGCAAACGTAGCATCTGTTTCTCATCTTcg                                                                                                                                                                                                            >  2:238807/1‑228 (MQ=255)
      cTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACGGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGCCCAACGTCGCACCAGGTtctctcactcgcgg                                                                                                                                                                                                         >  2:681886/1‑218 (MQ=255)
          gcAGTAACTTGTTACAACCTTGAGCATCCAGTTGCCGGTCGTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAata                                                                                                                                                                                                     <  2:329550/227‑1 (MQ=255)
           cAGTAATGTGTTACAAGCTGTAGCATCCACTTGCCGGTCCTGGGAGTTAAAAGGGAATCCAGTGCGAGTCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGTACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAata                                                                                                                                                                                                     <  2:292221/227‑1 (MQ=255)
              tAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGCCCAAAGTCGCCTCTGGGTAccatcaccccgcaaattcct                                                                                                                                                                                                 >  1:416303/1‑214 (MQ=255)
                      tACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTATTAGGGAATCCAGTGCGAATCTGGAGCTGCCGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGTCATTCGTTGGGAAGTCATCATCTCTTAGTATCTGAGATACCCCTCCAAGCCCGAAGACCGGCCGGCCAACGTCGCACCAGGTACCCACCATCCCGCAAAATTGTTGACAccc                                                                                                                                                                                         >  1:329550/1‑227 (MQ=255)
                               tAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCTAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATcc                                                                                                                                                                                <  2:542891/228‑1 (MQ=255)
                               tAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATc                                                                                                                                                                                 >  2:79315/1‑227 (MQ=255)
                                aGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATcc                                                                                                                                                                                <  1:79315/227‑1 (MQ=255)
                                   aTCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCcttct                                                                                                                                                                            >  2:40542/1‑228 (MQ=255)
                                      cACTTGCCGGTTCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGGAAGGAGCGATTATTGCGTTATGCGGACATGGCCATCCGGTGGGAAGTGTTCATCTCTTAGTATTTTAGATATCACTCCAAGCCCGAAGACCTCCCGGAAAACGACCCATCTGGTCATCACCATCGCGTAATATTGATAAAAACTTCCGCACCCTtcttct                                                                                                                                                                         >  1:665902/1‑228 (MQ=255)
                                       aCTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCACCCTTCTTCTa                                                                                                                                                                        >  2:556836/1‑228 (MQ=255)
                                         ttGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCCCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTAtt                                                                                                                                                                      >  2:343292/1‑228 (MQ=255)
                                              gggccTGGGAGTTAATGGGGAATTCAGTGGGGAGGATGGGCTGACGCGCAGGTGTAAAGGAAGGGTCGATTATTGGGTTAGGCGGACACTGCCATTCGGTGGGAAGTCCTCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAAACCTGCCGGCGAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCAGCATCCTTCTTCTATTGTGGa                                                                                                                                                                 <  1:678398/225‑1 (MQ=255)
                                              ggTCCTTTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTTACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCACCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTgg                                                                                                                                                                  <  2:465219/227‑1 (MQ=255)
                                                 ccTGTGAGTTAATAGGGAATCCAGGGCGAATCTGGTGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCTGACCCTGCCCTTCGGTGGGAAGGCATCTTCTCTTAGTATCTTAGATAACCCTCCAAACCCGAAGACCTGCAGGCCAACGTCTCCTCCGGTTCTCAACAACCCGTAACATTGACGAAACCCACGGccccatttcaccatttaggatgc                                                                                                                                                              >  1:302261/1‑206 (MQ=255)
                                                       aGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACa                                                                                                                                                        <  2:276935/228‑1 (MQ=255)
                                                        gTTAATGGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGTAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACa                                                                                                                                                        <  2:92353/227‑1 (MQ=255)
                                                        gTTAATAGGGAATTCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACa                                                                                                                                                        <  2:663990/227‑1 (MQ=255)
                                                        gTTAAAGGGGAATCCAGTGGGAATCTGGAGCTGACGCGCCGCGGGAAGGAAAGGTGTGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCCTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACa                                                                                                                                                        <  1:681886/227‑1 (MQ=255)
                                                           aaaggggAATCCATGGTGAATCTGGGGCTGACGCGCAGGGGGAAGGAAAGGTGCGATGGTTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATg                                                                                                                                                     <  1:40542/223‑1 (MQ=255)
                                                               gggAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGAtt                                                                                                                                                  >  2:169564/1‑226 (MQ=255)
                                                               gggAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGGGGGAAGTCATCATCTCTTAGTATCTTTGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCCTCATCCCGTAACATTGCGGAAACCTGCGGCATCCTTCTTCTATTTTGCTTGCTTTACAATGATaaa                                                                                                                                                >  1:201538/1‑228 (MQ=255)
                                                                ggAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGAGCGATGATTGCGTGATGCGGACACTGCCATTCGGTGGAAAGTCATCATCTCTCAGTATCTTGGATACCCCTCCAAGCCCGAAGCCCTGCCGCGCAACGTCGCATCTGGTTCTCATCACCCCTGAATATTGATGAAACCTGCGGCAACCGTCTTCTATTGGCGATGCTTTACAGTGAttt                                                                                                                                                 >  2:306291/1‑225 (MQ=255)
                                                                 gAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGCTGAAACCTGCGCCCTCCTTCTTCTTTTGTGGCTGCTTTCAAATGATTaaa                                                                                                                                               >  1:499035/1‑227 (MQ=255)
                                                                   aTCCAGTGCGAGTCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGCTGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGGCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCAGGCCCACGTCGCATCTGGTTCTCCTCATCGCCGAATATCGATGAAACCTGCGGCTACCTTCTCCTATTGTGGTTGTTTAACACAGATTaaa                                                                                                                                               >  1:625311/1‑225 (MQ=255)
                                                                    tCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATtgat                                                                                                                                                                                                <  2:590174/175‑1 (MQ=255)
                                                                         tGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAAGGTCGCATCTGGTTCTCTTCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTTTGGATGCTTTACAATGATTAAAAACGCTTcc                                                                                                                                      >  1:476677/1‑227 (MQ=255)
                                                                          gCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCCCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCgc                                                                                                                                     >  2:95304/1‑228 (MQ=255)
                                                                                                           aGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGGGGGAAGTCATCATATCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCCCATCAGGTTCTCATCATCGCGTAATCTTGATGAAACCTGCGGCATCCTTCTTCTCTTTCGGATGCTTTTCCATGATATAAAAAGCTTCTTTGCTGTCGGCGGGTTTCGTCAAGGCATTTTcac                                                                                                    >  1:425803/1‑226 (MQ=255)
                                                                                                                 gtggtttGGGTTTTGCGGAACGTGCCGTTCGGTGGGAAGCCATCATCTCTTAGTATATTAGATCCCCCGCCAAGCCCGAAGACCGGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCGTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTggg                                                                                              <  1:60948/223‑1 (MQ=255)
                                                                                                                    gtttgacttgttgtgatacggCCATTCGGTGGAAAGTCTTCACCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGca                                                                                            <  2:676682/207‑1 (MQ=255)
                                                                                                                    gatTGCGGTATGCGGACCCTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGca                                                                                            <  2:499035/227‑1 (MQ=255)
                                                                                                                    gatGGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATCCCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGca                                                                                            <  2:476677/227‑1 (MQ=255)
                                                                                                                     atTGGGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTTTTAGTAACTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGca                                                                                            <  2:675562/226‑1 (MQ=255)
                                                                                                                      ttGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATTTCTTAGTATCTTAGATACCCCTTCAAGGCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGcaca                                                                                          <  2:425803/227‑1 (MQ=255)
                                                                                                                            ttttcggaacgtGCCATTCGGTGGGAAGTCCTCACCTCTTAGTAACTTGGATCCCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATAc                                                                                    <  1:5860/216‑1 (MQ=255)
                                                                                                                              tGCGGACCCTGCCATTTGGTGGGAAGTCCACATCTCTTAGTATCTTAGATTCCCCTCCAAGCCCCAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAg                                                                                 <  2:685681/228‑1 (MQ=255)
                                                                                                                                    cccgtcccttCGGTGGGAAGTCATCACCTCCTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCCCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGa                                                                           <  1:343292/220‑1 (MQ=255)
                                                                                                                                     gcTGCCATTCGGTGGGAAGTTATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCGTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGa                                                                           <  2:178542/226‑1 (MQ=255)
                                                                                                                                            ttCGGTGGGCAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCtcatc                                                                                                                                                                                                                >  1:597839/1‑87 (MQ=255)
                                                                                                                                             tCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCtcatca                                                                                                                                                                                                               <  2:597839/87‑1 (MQ=255)
                                                                                                                                                 tGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCACGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAAGATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGTTGTCGTCGCGTTCCGTGAAGGCGTTTTCCGCTTCGTCACCAAATACCAACCCGGATaaagatgtcg                                                                >  2:264981/1‑217 (MQ=255)
                                                                                                                                                          atcatcTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGAGACTCTCGACGTTACTCCTacc                                                     >  2:430589/1‑226 (MQ=255)
                                                                                                                                                                    aGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCCCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGTCGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCAAAGGCAACCAACCCTGTTACCCTCGTCGTTCCTGACAACCGTTTTGAAc                                           >  1:689338/1‑228 (MQ=255)
                                                                                                                                                                    aGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAAAGTCGCATCTGGTTCTCATCATCGCGTCATATTGATGTAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGCCGGCGTGTTCCGTCACGGCCTTTTCCGCTTGGCCATAGGTTACCAGCCCGGCTATTCTCGCCGTTCATGCTGCCCGTTTTGaa                                            >  2:471596/1‑227 (MQ=255)
                                                                                                                                                                      tagttttgATTCCCCCCAAACCCCGAAGACCTGCCGGCCAACGTAGCATCTGGTTCTCATCAGCGTGTAATATTGATGAAACCTGCGGCCTCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAg                                         <  1:235706/221‑1 (MQ=255)
                                                                                                                                                                      tATCTTAGATCCCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAg                                         <  2:355574/228‑1 (MQ=255)
                                                                                                                                                                         cTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACCCTCGCCGTTTCTGCTCACCGTTTTGAACAGCcg                                      >  2:390217/1‑228 (MQ=255)
                                                                                                                                                                           tgattcccccTCCAAGCCGCAAACCCTGCCGGCCAACGTCGCCTCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCcgcg                                    <  2:625311/222‑1 (MQ=255)
                                                                                                                                                                           tAGATCCCCGTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCAACTGGTTCTCATCATCGCGTAATAATGATGAAACCTGCGGCATCCTGCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCGCGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCcgcg                                    <  2:658672/228‑1 (MQ=255)
                                                                                                                                                                               tACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGCTACCACCCGGGATACTCTCGTCGTTACTGCTAACCGGTTTGAACAAccagccacc                                >  1:184398/1‑221 (MQ=255)
                                                                                                                                                                                aCCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCCCTGCTGACTGCGTGTTCCGTCACGGCTTTTTCCCATTGGTCACAGGTGACCACCCCGGCTACTCTCGCCTTTACTGCTAACCCTTTTGAAAACCAGCgcacaa                               >  1:339526/1‑225 (MQ=255)
                                                                                                                                                                                 ccccccccAGCCCGAAACCCTGCCGGCCAACGTCGCATTTGGTTCTCATTATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCAc                              <  1:95304/228‑1 (MQ=255)
                                                                                                                                                                                 ccccccCAAGCCCGAAGACCTGCCGGCCAACCGCGCAACTGGTTCTCATCATCGCGTAATTATGATGAAACCTGCGGCAACCTTCTTCGATTGTGGATGCTTTACAATGATTAAAAAAGCGTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCAc                              <  2:302261/228‑1 (MQ=255)
                                                                                                                                                                                 ccccTTCACTCCCCAAGACCCTGCGGGCAACGGCGCGTCTGGTTCTCAACATCGCTTAATATTGATGAAACCTGCGGCTTCCTTCATCTAGTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCAc                              <  1:438821/228‑1 (MQ=255)
                                                                                                                                                                                 ccccTCCCACCCCCAAGACCTGCCGGCCCACGTCGCATCTGGTTCTCATCGTCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCAc                              <  1:556836/228‑1 (MQ=255)
                                                                                                                                                                                  ccccccAAGCCCGAAGGCCTGCCGGCCAACGTTGCATCTGGGTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCAc                              <  2:223811/227‑1 (MQ=255)
                                                                                                                                                                                  cccTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACCGCTAACCGGTTTGAACAGCCGCGCCGCCCt                             >  2:28470/1‑228 (MQ=255)
                                                                                                                                                                                  cccTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTCCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCCtttctcttataaacatctccg                              >  2:6059/1‑209 (MQ=255)
                                                                                                                                                                                  cccTCCAAGCCCGAAGACCTGCCGGCCAAAGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCCGCATCCTTCTTCTATTGTGGCTGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGAGTTTTCCGCTTGGGCTCAGGATACCAGTCCGTATACACCCGCGGTTACTCCTACACGGTTTGAAAACCCGCgcacaacg                             >  1:225/1‑223 (MQ=255)
                                                                                                                                                                                    cTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGa                                                                           <  2:420267/180‑1 (MQ=255)
                                                                                                                                                                                    cTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGa                                                                           >  1:420267/1‑180 (MQ=255)
                                                                                                                                                                                    cTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGGGTTCCGTCACGGCATTTTCCGCTTGGGCCCAGGATACCAGCCCGGATACTACCGTCGTTCCGGCTAACCGTTTTGAACAGCCGCGCCGCAcct                            >  2:235706/1‑225 (MQ=255)
                                                                                                                                                                                    cTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTCACGGCGTGTTCCGGTACGGCATTTTCCGCTTGGGCACAGGAACACAGCCCGGTAACTTCCCTCCGTCATTCTACCCGTTTTGAAAAGCCGCccagcaacag                           >  1:660148/1‑224 (MQ=255)
                                                                                                                                                                                         aGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGCTGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATAATCTCGCCGTTACTGCTAACCGTTTTGAACAGCCGCGCCGCACTTGGCTTg                      >  2:640227/1‑228 (MQ=255)
                                                                                                                                                                                         aGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTACTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATCAAAAAAGCTTCGCTGCTGACGGCGCGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAACCCgg                                                                            >  2:589630/1‑174 (MQ=255)
                                                                                                                                                                                         aGCCCGAAGACCTGACGGCAAACGTCGCATCTGGTTCCCATCATCGCATAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTg                      <  2:689338/228‑1 (MQ=255)
                                                                                                                                                                                            ccGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGCAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAACAACCTTCGCTGCTGTCGGCTTGTTCCGTCACGGCATTTTCCGCTTTTGCACAGGATACCAGCCCCGACACTCTCGACGTTACTGCTCACCGTTTTTAACCGCCGCCGGGAACTGTGGTTgacc                   >  1:564793/1‑225 (MQ=255)
                                                                                                                                                                                            ccGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGGGCAGCACTGGCCTTGccc                   >  2:438550/1‑226 (MQ=255)
                                                                                                                                                                                            ccGAAAACCGTCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCa                    <  1:437923/227‑1 (MQ=255)
                                                                                                                                                                                              gAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAACCGCGCAGCACTAGGCTTGCacca                 >  1:84131/1‑228 (MQ=255)
                                                                                                                                                                                               aaGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGCTACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCCCTGTGCTTGc                     >  2:437923/1‑223 (MQ=255)
                                                                                                                                                                                                         ggccaacgcgccacttgtttCTCATCCTCGCGTAATATTGATGAAACCTGCCGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTtgt      <  1:169564/210‑1 (MQ=255)
                                                                                                                                                                                                               cggcgtttttGGTTTTCATCTTCGAGTAATTCTGATGAAACCTGCGGCATCCTTTTTCTATGGTGGATGCTTTACAATGATTAAAAAAGCTTCGGTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGAcc  <  1:264981/217‑1 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                  
CTACGTCTGGACAGTAACTTGTTACAACCTGTAGCATCCACTTGCCGGTCCTGTGAGTTAATAGGGAATCCAGTGCGAATCTGGAGCTGACGCGCAGCGGTAAGGAAAGGTGCGATGATTGCGTTATGCGGACACTGCCATTCGGTGGGAAGTCATCATCTCTTAGTATCTTAGATACCCCTCCAAGCCCGAAGACCTGCCGGCCAACGTCGCATCTGGTTCTCATCATCGCGTAATATTGATGAAACCTGCGGCATCCTTCTTCTATTGTGGATGCTTTACAATGATTAAAAAAGCTTCGCTGCTGACGGCGTGTTCCGTCACGGCATTTTCCGCTTGGGCACAGGATACCAGCCCGGATACTCTCGTCGTTACTGCTAACCGTTTTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACC  >  NC_000913/4163355‑4163787

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