Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,421,374 C→T S17L (TCA→TTA)  yhdY → putative ABC transporter membrane subunit YhdY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,3740CT98.4% 217.3 / NA 62S17L (TCA→TTA) yhdYputative ABC transporter membrane subunit YhdY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/32);  minor base G (0/1);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGC  >  NC_000913/3421134‑3421619
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
atgacctctatgagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCcgcgcg                                                                                                                                                                                                                                             <  1:533162‑M2/18‑1 (MQ=255)
       ctatgagctggtactggctgaagtagaacngcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGgt                                                                                                                                                                                                                                        <  2:288919‑M2/23‑1 (MQ=255)
           gagctggtactggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGt                                                                                                                                                                                                                                  <  2:7602‑M2/29‑1 (MQ=255)
                   actggctgaagtagaacagcccctgttgganatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaaa                                                                                                                                                                                                                           <  2:253188‑M2/36‑1 (MQ=255)
                    ctggctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGaaaaa                                                                                                                                                                                                                          >  1:223975‑M2/214‑250 (MQ=255)
                       gctgaagtagaacagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCt                                                                                                                                                                                                                      >  1:457637‑M2/211‑251 (MQ=255)
                                  acagcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCTTAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCt                                                                                                                                                                                                           <  2:452098‑M2/52‑1 (MQ=255)
                                     gcccctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGa                                                                                                                                                                                                        <  1:494585‑M2/55‑1 (MQ=255)
                                        cctgttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgtgcatcaaccgtggtacgctgcgtaaaaaagtgaaaaaatacggcatgaactaattcaggttggctaaatgcttgattaaaaaggcgctactcggcatggggaagcgcctgttttataggtttcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCa                                                                                                                                                                                                     <  2:194303‑M2/58‑1 (MQ=255)
                                           gttggacatggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAAt                                                                                                                                                                                                   <  1:249536‑M2/60‑1 (MQ=255)
                                                 catggtgatgcaatacacccgtggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgct                                                                                                                                                                                           >  1:13896‑M2/185‑251 (MQ=255)
                                                                 ccccgtggtaaccagacccgtgctgcgctgangatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCat                                                                                                                                                                            <  2:293789‑M2/83‑1 (MQ=255)
                                                                 ccccgtggtaaccagacccgtgcggcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCat                                                                                                                                                                            <  2:457637‑M2/83‑1 (MQ=255)
                                                                       ggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGt                                                                                                                                                                       <  1:216299‑M2/88‑1 (MQ=255)
                                                                       ggtaaccagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattaaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGtt                                                                                                                                                                      <  2:136886‑M2/89‑1 (MQ=255)
                                                                             cagacccgtgctgcgcttatgatgggcatcanccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcgtgattaaaaaggcgctactcggcatggggaagcgcctttgttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCATCTGGAGCAATAGCCTGCTGTCTATTGGCTGCATATGGTTGAtgtg                                                                                                                                                                <  2:273380‑M2/95‑1 (MQ=255)
                                                                             cagacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGAtgtg                                                                                                                                                                >  1:510157‑M2/157‑251 (MQ=255)
                                                                               gacccgtgctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaagacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                                                                                                                                              <  1:104395‑M2/97‑1 (MQ=255)
                                                                                    gtgatgcgatgatgatgggcatcaaccgtggtacgttgcgtaaaaaattgaaaaaatacggcatgaactaattcgggttagctaaatgcttgattaaaaaggcgctactcggcattgagaagcgccttttgtataggtgtcacaaaggaCATACCTGAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTg                                                                                                                                                         <  2:336064‑M2/102‑1 (MQ=255)
                                                                                      gctgcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATc                                                                                                                                                      <  2:13896‑M2/104‑1 (MQ=255)
                                                                                         gcgctgatgatgggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATTTGGGAACTGATc                                                                                                                                                      <  2:218573‑M2/105‑1 (MQ=255)
                                                                                                 gatgggcatcaaccgtggtacgctgcgtaaaaaagtgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttgtttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTg                                                                                                                                             <  2:25417‑M2/114‑1 (MQ=255)
                                                                                                    gggcatcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTg                                                                                                                                          <  1:569216‑M2/117‑1 (MQ=255)
                                                                                                       catcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggctagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGGGTGGGTGCGAAACAATCTGTTCTCCAGCTGGGGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGGTGTGGGTACTGTACCCACAGTTGCTTAACt                                                                                                                                      >  2:3722‑M2/131‑251 (MQ=255)
                                                                                                        atcaaccgtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTg                                                                                                                                     <  1:151237‑M2/122‑1 (MQ=255)
                                                                                                               gtggtacgctgcgtaaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCAttt                                                                                                                              <  2:223975‑M2/129‑1 (MQ=255)
                                                                                                                             aaaaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGCGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTg                                                                                                                >  1:245753‑M2/109‑251 (MQ=255)
                                                                                                                               aaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTg                                                                                                                                                                     >  1:368642‑M2/107‑196 (MQ=255)
                                                                                                                               aaaattgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTg                                                                                                                                                                     <  2:368642‑M2/90‑1 (MQ=255)
                                                                                                                                  attgaaaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTg                                                                                                           <  1:425810‑M2/148‑1 (MQ=255)
                                                                                                                                       aaaaatacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTc                                                                                                       >  1:483373‑M2/99‑250 (MQ=255)
                                                                                                                                           atacggcatgaactaattcaggttagctaaatgcttgattaaacaggcgctactcggcatggggaagcgccttttttataggtgtcacaacggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATACCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAAcgc                                                                                                  >  1:208810‑M2/95‑251 (MQ=255)
                                                                                                                                           atacggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAAcgc                                                                                                  >  1:130621‑M2/95‑251 (MQ=255)
                                                                                                                                              cggcatgaactaattcaggttagctaaatgcttgattaaaaaggcgctactaggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAAATGGGTTGGCTCAACGCGTg                                                                                               >  1:69960‑M2/92‑251 (MQ=255)
                                                                                                                                                  atgaactaattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGCCTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTACCTGGGTTGGCTCAACGCGTGCCg                                                                                            >  2:62597‑M2/88‑250 (MQ=255)
                                                                                                                                                        taattcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGCCTGcac                                                                                     >  1:390949‑M2/82‑251 (MQ=255)
                                                                                                                                                            gcaggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgtcttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACaaa                                                                                  <  1:494288‑M2/173‑1 (MQ=255)
                                                                                                                                                              aggttagctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGc                                                                               >  1:517356‑M2/76‑251 (MQ=255)
                                                                                                                                                                gttagctacatgcttgattaaacaggcgctactcggcatggggaagagccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGATGATGTGGGAACTGATCCCACCTTTGCTGAACTAGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCg                                                                             >  2:449809‑M2/74‑251 (MQ=255)
                                                                                                                                                                   agctaaatgcttgattaataaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTg                                                                          >  2:481877‑M2/71‑251 (MQ=255)
                                                                                                                                                                    gctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCCGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTAAACGCGTGCCGACTGCACAAAAGCCGGTGc                                                                         >  1:316731‑M2/70‑251 (MQ=255)
                                                                                                                                                                     ctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGAc                                                                                          >  2:205431‑M2/69‑233 (MQ=255)
                                                                                                                                                                     ctaaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGAc                                                                                          <  1:205431‑M2/165‑1 (MQ=255)
                                                                                                                                                                       aaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGcc                                                                        <  1:563321‑M2/183‑1 (MQ=255)
                                                                                                                                                                       aaatgcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCt                                                                       <  2:447168‑M2/184‑1 (MQ=255)
                                                                                                                                                                           gcttgattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTgg                                                                  >  1:226916‑M2/63‑251 (MQ=255)
                                                                                                                                                                             ttgattaaaaaggcgctactcggcctggggaagcgccttttttataggtgtcacacaggaCATAACTTAAGCAGCGCGATGGAGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGAATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTAGGCTCAACGCGTGCCGACTGCACAAAAGACGGTGCCTGTtgggg                                                                >  1:459908‑M2/61‑250 (MQ=255)
                                                                                                                                                                               gattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGc                                                                                                                                                                                 <  2:291496‑M2/78‑1 (MQ=255)
                                                                                                                                                                               gattaaaaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGc                                                                                                                                                                                 >  1:291496‑M2/59‑136 (MQ=255)
                                                                                                                                                                                     aaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATcc                                                        >  2:353529‑M2/53‑251 (MQ=255)
                                                                                                                                                                                     aaaggcgctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATcc                                                        >  1:372283‑M2/53‑251 (MQ=255)
                                                                                                                                                                                           gctactcggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCacgaac                                                  >  1:417123‑M2/47‑251 (MQ=255)
                                                                                                                                                                                                  ggcatggggaagcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTgg                                           >  1:417905‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                             gcgccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgt                                >  1:312890‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                ccttttttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTtatgtatg                             <  1:212518‑M2/226‑1 (MQ=255)
                                                                                                                                                                                                                     tttataggtgtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCt                         <  1:94385‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                                              gtcacaacggaCATAAATTAAGCCGCGCGCTGGTGTGGGTGCGAAAAAATCTGTTATCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGTTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGGGCCGACTGCACAAAAGACGGTGCCTGTTGGGGCTTCATCCACGTACGATNTGGTCAGATTAATAATGGGCTTTACCAacac               >  2:271814‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                              gtcacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCTAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTACCTTGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTttaccctca                <  1:51892‑M2/239‑4 (MQ=255)
                                                                                                                                                                                                                                cacaaaggaCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTATCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGGGGGAACTGAACCCACCGTTGCTGAACTGGGCATTTTTGCATGCTAAATGGGTTGGCTCAACGCGTGCCGACGGCACAAAAGCCGATGCCTGGTGGGTCTTCATCCACGAACTATTTGGTCAGGTTATGTAGGGgctttaccccaacgc             >  2:548116‑M2/10‑245 (MQ=255)
                                                                                                                                                                                                                                        aCATAACTTACGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAGCTGATCCCCCCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCt     <  1:4547‑M2/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aCATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCt     >  1:12047‑M2/2‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGc  <  1:341277/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
CTCTTGCCGCCGCTATTGGCTATCCCGATATGGTTTCGCTGTTTGCCGGCACCGTGCTGAATCAGACGGGGCAAGCCATCGAGACGATAGCCATGACCATGTCGGTCTATCTGATTATCAGCCTGACTATCTCGCTGCTGATGAATATCTATAACCGCCGCATCGCGATCGTTGAACGCTAAGGAGTCATGATGACAAAAGTATTGCTGTCTCATCCCCCGCGCCCGGCGAGCCATAACTCAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGCAGGCTAACTGGGTTGGCTCAACGCGTGCCGACTGCACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGC  >  NC_000913/3421134‑3421619

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