Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 143,533 T→C 5.7% L252P (CTG→CCG)  yadG → putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913143,5330TC 5.7% 236.6 / 3.6 88L252P (CTG→CCG) yadGputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base C (3/2);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

GCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACTGGTTAA  >  NC_000913/143386‑143680
                                                                                                                                                   |                                                                                                                                                   
gCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTGATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTg                                                                                                                                                    <  1:36297/149‑1 (MQ=255)
  cGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGGTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGt                                                                                                                                                  >  1:1019062/1‑149 (MQ=255)
  cGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCCGCTGGCGAAGCTGAAATCGGAAACCTCTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCCGTAGGGACTGGt                                                                                                                                                  <  2:485127/149‑1 (MQ=255)
   gCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCCAGATGGCCATCCTTCTCGACCggcg                                                                                                                                                 >  2:433059/1‑147 (MQ=255)
    cAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTCCCGAAGCTCGATGGCTATCAGAACCGACGAGTCg                                                                                                                                                >  2:985850/1‑149 (MQ=255)
    cAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCg                                                                                                                                                >  1:408908/1‑149 (MQ=255)
         tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGCCTGGTCGATAcc                                                                                                                                           >  2:273121/1‑149 (MQ=255)
           ggCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGCCTGGTCGATACcgc                                                                                                                                         >  1:789146/1‑149 (MQ=255)
           ggCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACcgc                                                                                                                                         >  1:841/1‑149 (MQ=255)
                  ttCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTg                                                                                                                                  >  2:1017170/1‑149 (MQ=255)
                   tCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTgg                                                                                                                                 >  1:281687/1‑149 (MQ=255)
                     aaCACGGTGAGCTGGTGGCAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGaa                                                                                                                               >  2:743798/1‑149 (MQ=255)
                     aaCACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGaa                                                                                                                               >  1:751613/1‑149 (MQ=255)
                     aaCACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGaa                                                                                                                               >  2:743830/1‑149 (MQ=255)
                          ggTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCCCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGa                                                                                                                          <  2:306378/149‑1 (MQ=255)
                            tGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAg                                                                                                                        <  1:574621/149‑1 (MQ=255)
                             gAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTGTTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCCCGATGGCTATCAGTATCGACTGGTCGATACAGAGACGCGGGAAGTTGCAGt                                                                                                                       >  2:157831/1‑149 (MQ=255)
                                cTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACCAGTCGATActgata                                                                                                                                       <  1:324416/130‑6 (MQ=255)
                                cTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACCAGTCGATActgata                                                                                                                                       >  2:324416/1‑125 (MQ=255)
                                cTGGTGGAAAATACCTCGACGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACCAGTCGATActgata                                                                                                                                       <  1:324421/130‑6 (MQ=255)
                                cTGGTGGAAAATACCTCGACGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACCAGTCGATActgata                                                                                                                                       >  2:324421/1‑125 (MQ=255)
                                  ggtggAAAATACCTTGATGAAGGCGCTGCGGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGtgctgc                                                                                                                  <  1:14149/149‑1 (MQ=255)
                                  ggtggAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGtgctgc                                                                                                                  <  1:1076793/149‑1 (MQ=255)
                                    tggAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGt                                                                                                                <  2:1068520/149‑1 (MQ=255)
                                     ggAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTg                                                                                                               >  1:652510/1‑149 (MQ=255)
                                     ggAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATAGGCAACCTTTATTCTCGATCTTGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTg                                                                                                               >  1:274688/1‑149 (MQ=255)
                                            cccTCGAGGACGCCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAggg                                                                                                        <  2:315940/148‑1 (MQ=255)
                                            aCCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAggg                                                                                                        <  2:653187/149‑1 (MQ=255)
                                            aCCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCGCGCACCGAAAAGCCCGTTACCGAAGCTCGCTGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAggg                                                                                                        <  1:731757/149‑1 (MQ=255)
                                              cTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGACGTGCTACGTGAGCCGGGGa                                                                                                      >  2:340192/1‑149 (MQ=255)
                                              cTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGa                                                                                                      <  2:1045186/149‑1 (MQ=255)
                                              cTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGa                                                                                                      >  1:1064119/1‑149 (MQ=255)
                                               tCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAATCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGAt                                                                                                     >  1:451775/1‑149 (MQ=255)
                                                 gATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCa                                                                                                   <  2:1060687/149‑1 (MQ=255)
                                                  agtaaGGCGGTGTTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATGTCGCACCGAAAAGCCCGTTACCGACGCTTCACGGCGATCAGTATCGACTGGTCGAGACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCaa                                                                                                  <  2:924262/146‑1 (MQ=255)
                                                    gAAGTCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACa                                                                                                <  1:1049742/149‑1 (MQ=255)
                                                    gAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACGTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACa                                                                                                <  1:478444/149‑1 (MQ=255)
                                                        gcgcTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGt                                                                                            >  1:729316/1‑149 (MQ=255)
                                                        gcgcTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGt                                                                                            <  2:408908/149‑1 (MQ=255)
                                                           ctgctgGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTAtt                                                                                         <  1:1017170/149‑1 (MQ=255)
                                                            tgctgGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTAttt                                                                                        >  1:50853/1‑149 (MQ=255)
                                                            tcttgGCCAAGCTGAAATCGGAAACCTTTAGTCCCGATCCCGCACCGAAAAGCCCGTTACCGAAGCCCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTAttt                                                                                        <  1:204578/146‑1 (MQ=255)
                                                               tgGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACg                                                                                     <  1:804572/149‑1 (MQ=255)
                                                                      gcggAAATCGGAAACCTTGATCCGCGGCCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGCCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTaa                                                                              <  1:991515/146‑1 (MQ=255)
                                                                      gCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTaa                                                                              <  1:676845/149‑1 (MQ=255)
                                                                      ccTGAAATCCGAAACCCTTCTTCTCGACCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGGTaa                                                                              <  2:1056277/148‑1 (MQ=255)
                                                                        tGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGACGCTCGATGGCTATCAGTATCGCCTGGTCGCTACCGCGAGGCTGGACGTTGCAGTGCTGCGTGAGCAGGGGATCAACAGCGGCTTTCCGCACATAAGt                                                                            >  2:9902/1‑149 (MQ=255)
                                                                        tGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGt                                                                            <  1:1093109/149‑1 (MQ=255)
                                                                        tGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGt                                                                            <  1:709976/149‑1 (MQ=255)
                                                                         gaatgcGGAAACCTTCATGCTGGATCTCGCACGGAAAAGCCCGTTACCGAAGCGCGATGGCTATAAGTATGGACTGGACGCTCACGTGGCGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTg                                                                           <  1:273121/144‑1 (MQ=255)
                                                                         gAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCACCAGCGTATTTACGAAGTTAGGTg                                                                           >  2:911912/1‑149 (MQ=255)
                                                                            aTCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGATCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGc                                                                        <  1:743798/149‑1 (MQ=255)
                                                                            aTCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGCGCTGCGTGACCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGc                                                                        <  1:743830/149‑1 (MQ=255)
                                                                                gAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGGGCCggg                                                                    >  1:1119145/1‑149 (MQ=255)
                                                                                gAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAggg                                                                    >  1:474783/1‑149 (MQ=255)
                                                                                     cTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTc                                                               >  2:134460/1‑149 (MQ=255)
                                                                                      tttATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGGCGGTACGTGGGCAGGGCATTCa                                                              >  2:956935/1‑149 (MQ=255)
                                                                                       ttATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAg                                                             <  1:444924/149‑1 (MQ=255)
                                                                                       ttATGCCCGGTCTCGCACCGAAAAGCCCGGTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAg                                                             <  1:348796/149‑1 (MQ=255)
                                                                                                aTCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGACGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGt                                                    >  1:624246/1‑149 (MQ=255)
                                                                                                   tcGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGGTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCAGTAAGCTATTAAGTATg                                                 >  1:912293/1‑149 (MQ=255)
                                                                                                     gCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTAGAAGTGCTGCGAGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAAGGATTAAGttagcg                                               >  2:954736/1‑144 (MQ=255)
                                                                                                     cccacGAAAGGCCCGTTCCCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCg                                               <  2:502497/145‑1 (MQ=255)
                                                                                                      cACCGAAAAGCCCTTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGt                                              <  2:281612/149‑1 (MQ=255)
                                                                                                      cACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGt                                              <  1:60015/149‑1 (MQ=255)
                                                                                                      cACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGt                                              >  1:495968/1‑149 (MQ=255)
                                                                                                      cACCGAAAAGCCCGGTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAACTGAGCAGGGCATTCAGGTATTAAGTATGCgc                                              <  2:843889/149‑2 (MQ=255)
                                                                                                       aCCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTa                                             >  2:733506/1‑149 (MQ=255)
                                                                                                         cGAAAAGCCCGTTACCGAAGCTCGATGGCTATGCGTAGCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAAc                                           <  2:919884/149‑1 (MQ=255)
                                                                                                            aaaGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACaaa                                        <  1:244594/149‑1 (MQ=255)
                                                                                                             aaGCCCGTTACCGAAGCTCGATGGCGATCATTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAg                                       >  1:596862/1‑149 (MQ=255)
                                                                                                                cccGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCCGGTATTAAGTATGCGTACAAAAGCTa                                    >  2:92553/1‑149 (MQ=255)
                                                                                                                 ccGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTaa                                   <  1:161004/149‑1 (MQ=255)
                                                                                                                       ccGAAGCTCGATGGCTAGCAGTATCGACTGGGCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCt                             <  1:342486/149‑1 (MQ=255)
                                                                                                                         gAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGGCGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGAAGTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTgg                           >  2:862276/1‑149 (MQ=255)
                                                                                                                           aGCTCGAGGGCTATCAGTATCCGCTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTTTTTACGCAGTTAAGTGCGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCAAACGGTCTGGaa                         <  2:1056779/149‑1 (MQ=255)
                                                                                                                             cTCGATGGCTATCAGGAGCGACTGGTCGATACCGGGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAaga                       <  2:281687/149‑1 (MQ=255)
                                                                                                                                gATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCt                    <  1:253907/149‑1 (MQ=255)
                                                                                                                                      tATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCtgtttgt              >  2:881529/1‑149 (MQ=255)
                                                                                                                                         cAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAACAGCTGTTTGTTTc           >  2:264596/1‑149 (MQ=255)
                                                                                                                                            tATCGACTGGTCGATACCGCGACGCTGGAGGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGAAGTTAAGTGAGCAGGGCCTTCTGGGATTAGGTATGCGGTACAATGCTAACCGTCTGGGAGGGCTGTTTGTTTCACt        >  2:634466/1‑149 (MQ=255)
                                                                                                                                            tATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGGGCTGTTTGTTTCACt        >  1:295671/1‑149 (MQ=255)
                                                                                                                                            tATCGACTGGTCGAGACCGCGACGCTGGAAGGTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACt        <  2:715970/149‑1 (MQ=255)
                                                                                                                                              tCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACtag      >  1:1011890/1‑147 (MQ=255)
                                                                                                                                                 aCTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACTGGTTa   >  2:1005758/1‑149 (MQ=255)
                                                                                                                                                 aCTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGAAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGCCTCGCAGAGCTGTTTGTTTCACTGGTTa   >  1:711305/1‑149 (MQ=255)
                                                                                                                                                  cTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACTGGTTaa  >  2:92058/1‑149 (MQ=255)
                                                                                                                                                  cTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACTGGTTaa  >  2:808012/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAGCCCGTTACCGAAGCTCGATGGCTATCAGTATCGACTGGTCGATACCGCGACGCTGGAAGTTGAAGTGCTGCGTGAGCAGGGGATCAACAGCGTATTTACGCAGTTAAGTGAGCAGGGCATTCAGGTATTAAGTATGCGTAACAAAGCTAACCGTCTGGAAGAGCTGTTTGTTTCACTGGTTAA  >  NC_000913/143386‑143680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑