Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,195 A→C 8.0% intergenic (+593/+21) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,1950AC 8.0% 192.9 / 7.5 73intergenic (+593/+21)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base C (3/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270063‑3270342
                                                                                                                                    |                                                                                                                                                   
cAACAGCAAAACCCCGCCGAACCGGGTTTTTACGTAAACCAGCTGAAACTGACCGATCAGCCGCTTTCTTTTGCGTCTAGCGGCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt                                                                                                                                     <  1:140297/149‑6 (MQ=11)
 aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                                >  6:374099/1‑137 (MQ=37)
 aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                                <  5:374099/137‑1 (MQ=37)
       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCGTTGTGCGCCGGCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                                                              >  7:313147/1‑149 (MQ=11)
                       gggTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCGTTCTTTTGGGTATAGCGTCGTGGCCAGGCATTCATCTTTGTGCGCCGGCAAACGCAAAAACCCGCCGAAGCGGGTTTTTACCTAAATCACGGGCAACTGACCGaa                                                                                                                >  1:60077/1‑146 (MQ=11)
                                 gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt                                                                                                                                          <  1:320078/111‑1 (MQ=11)
                                 gTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCATCGGCGGGttttt                                                                                                                                          >  2:320078/1‑111 (MQ=11)
                                   aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                                                              >  2:373242/1‑121 (MQ=17)
                                   aaaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                                                                              <  1:373242/121‑1 (MQ=17)
                                     aacAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCCTCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAc                                                                                                                                        <  5:302958/107‑1 (MQ=11)
                                              aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                                >  4:159799/1‑92 (MQ=33)
                                              aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                                <  3:159799/92‑1 (MQ=33)
                                               aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa                                                                                          <  8:182618/145‑1 (MQ=255)
                                               aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACa                                                                                          >  7:182618/1‑145 (MQ=255)
                                                cTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGACCAAACAAAAACAAAAACCCGCCGAAGCAGCTTTTTAAGTAAATCAGGTGAAACTAACAGATAATCCGCGTTCTTTAGTGGACAGtca                                                                                      >  8:383089/1‑148 (MQ=14)
                                                   cgccaTACGCGGCTTTCATTTGTGACTTGCGTCGAGGCGTGTGATTCACCTATCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtcattc                                                                                   <  3:98187/144‑1 (MQ=17)
                                                     cGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCAATCATAATTCTGCCCCAACAAAAACAAAAAACCGCCGAAACGGGTTTTTAAGTAAATAAGGTGAAAATGACCGATAAGCCCGGTTATGTCGTGGACAGtcattaat                                                                                 >  8:18064/1‑148 (MQ=17)
                                                            ccGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                           >  1:262897/1‑147 (MQ=255)
                                                            ccGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCAGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                           >  2:113587/1‑147 (MQ=255)
                                                            ccGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCATTCATCTTTCTGCCCCTCCAAAATCAAAAACCCGCCGAAGAGGGTTTTTACGTAAATCACGTGAAACTGACCGTGAAGACGGGTTATGTGGTGAAAATTAATTCATATAGGc                                                                           >  6:120975/1‑147 (MQ=14)
                                                                   ctttTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATc                                                                   <  5:120975/148‑1 (MQ=255)
                                                                    ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                                                                                >  7:181391/1‑134 (MQ=255)
                                                                    ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATc                                                                                <  8:181391/134‑1 (MQ=255)
                                                                    ttttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTACAAAAGCAAAAACCCGCCGCAGCGGGTTTTTACGTAAATGAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGGCAGTCATTCATCTAGGCCAGCCATCg                                                                  >  1:374665/1‑148 (MQ=18)
                                                                      ttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAgg                                                                            <  6:199970/136‑1 (MQ=255)
                                                                      ttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAgg                                                                            >  5:199970/1‑136 (MQ=255)
                                                                          gTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAACAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                          >  5:65100/1‑149 (MQ=255)
                                                                           tataGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAATCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTgg                                                          <  7:18064/149‑1 (MQ=255)
                                                                               gcctcgcGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGc                                                         <  1:253159/143‑1 (MQ=255)
                                                                                    gtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagcag                                                 >  5:293084/1‑149 (MQ=255)
                                                                                      ggtatgTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAagcagc                                                <  7:276043/143‑1 (MQ=255)
                                                                                          aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGcc                                                                                                           <  2:272600/85‑1 (MQ=33)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                                                                                                                    >  1:115039/1‑71 (MQ=34)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAc                                                                                                                    <  2:115039/71‑1 (MQ=34)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTATGCCAGCACTCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                      >  6:57410/1‑149 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                      >  7:107140/1‑149 (MQ=255)
                                                                                               ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGt                                      <  2:262897/149‑1 (MQ=255)
                                                                                                  atCTTTCTGCCCCTCCAAAAGCAAAAGCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     <  5:271636/147‑1 (MQ=255)
                                                                                                  atCTTTCTGCCCCTCCAAAAGCAAAAGCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  6:271636/1‑147 (MQ=255)
                                                                                                  atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  2:204036/1‑147 (MQ=255)
                                                                                                  atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  7:223532/1‑147 (MQ=255)
                                                                                                     tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTc                                                       >  1:365127/1‑126 (MQ=255)
                                                                                                     tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTc                                                       <  2:365127/126‑1 (MQ=255)
                                                                                                     tttCTGCACCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTa                                >  1:363881/1‑149 (MQ=255)
                                                                                                       tCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGTGTTTTTACGTAAATCAGGTGAAACTGACCGATAATCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGt                                 <  2:374665/146‑1 (MQ=255)
                                                                                                       tCTGACCATCCAAAAGCATATACCCGCCGTTGGGGGTTTTTACGTAAATCAGGTGTAACTGACAGATAAGCGGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGt                                 <  7:383089/146‑1 (MQ=37)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                                                            >  8:372318/1‑68 (MQ=35)
                                                                                                          gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                                                            <  7:372318/68‑1 (MQ=35)
                                                                                                           ccccTCCAACAGCAATCTCTCGCGTAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACggg                            <  2:60077/147‑1 (MQ=38)
                                                                                                           ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGcc                          >  8:237187/1‑149 (MQ=255)
                                                                                                           ccccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGc                           <  8:107140/148‑1 (MQ=255)
                                                                                                             ccTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCg                         >  4:167272/1‑148 (MQ=255)
                                                                                                              cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTa                       >  4:125378/1‑149 (MQ=255)
                                                                                                               tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                         >  1:304009/1‑98 (MQ=255)
                                                                                                               tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCa                                                                         <  2:304009/98‑1 (MQ=255)
                                                                                                                  aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                            >  5:265126/1‑108 (MQ=255)
                                                                                                                  aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACt                                                            <  6:265126/108‑1 (MQ=255)
                                                                                                                  aaaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCt                    >  1:148466/1‑148 (MQ=255)
                                                                                                                    aaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTAt                 <  8:313147/149‑1 (MQ=255)
                                                                                                                     aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg                <  8:223532/149‑1 (MQ=255)
                                                                                                                     aGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATg                <  3:125378/149‑1 (MQ=255)
                                                                                                                      gtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     >  1:155710/3‑127 (MQ=38)
                                                                                                                      gtaaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGtt                                     <  2:155710/125‑1 (MQ=38)
                                                                                                                      gCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGa               <  7:237187/149‑1 (MQ=255)
                                                                                                                       cataaaCCCGCCGAAGCGGGTTTTTACGTAAATCATGTTAAACTGAGCGATAAGGCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa              <  2:363881/146‑1 (MQ=39)
                                                                                                                       cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTCCCCGGGTTCAGTACGGGCCGTACCTTATGac              >  6:180458/1‑148 (MQ=255)
                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGcc                          <  1:229206/136‑1 (MQ=255)
                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGcc                          >  2:229206/1‑136 (MQ=255)
                                                                                                                        aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc             <  6:293084/149‑1 (MQ=255)
                                                                                                                              ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt       <  1:204036/149‑1 (MQ=255)
                                                                                                                              ccgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAt       >  4:85801/1‑149 (MQ=255)
                                                                                                                                  cgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTgg   >  7:82578/1‑149 (MQ=255)
                                                                                                                                   gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGc  <  3:85801/149‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                                   
CAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270063‑3270342

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: