Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,163,578 T→C 6.7% 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,5780TC 6.7% 178.9 / 2.2 67intergenic (‑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 (30/32);  new base C (3/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

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

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: