Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,155,322 G→C 5.3% D67H (GAT→CAT)  tmk → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,155,3220GC 5.3% 247.8 / 4.0 75D67H (GAT→CAT) tmkhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/39);  new base C (2/2);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

CGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATT  >  NC_000913/1155177‑1155468
                                                                                                                                                 |                                                                                                                                                  
cgcgcgTAATGTGGTGGTTGAGACGCTCGAGCACCTGGGTATCCGCGACATGGTTTTCCCTCGGGAACCTGGCGGGACGCAACTTGCCGAAAAGTTAAGCAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATaa                                                                                                                                                >  2:438283/1‑150 (MQ=255)
      tAATGTGGTGGTTGATACGCTCGAGCAACTGGGTATCCGCGACATTGTTTTCACTCGGGAACGTGGCGGTACGCAACTTGCCGAAAAGTTAAGAATACTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCTATAAAGCCGa                                                                                                                                          <  2:321466/150‑1 (MQ=255)
      tAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGa                                                                                                                                          <  2:407204/150‑1 (MQ=255)
           gggtggtTGAGACGCTCGAGCAACTGGGTTTCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGAAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTc                                                                                                                                     <  2:1381576/149‑1 (MQ=255)
            ggtggtTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCt                                                                                                                                    >  2:362987/1‑150 (MQ=255)
             gtggtTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTg                                                                                                                                   <  2:361662/150‑1 (MQ=255)
             gtggtTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTg                                                                                                                                   <  2:446723/150‑1 (MQ=255)
             gtggtTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTg                                                                                                                                   >  1:195900/1‑150 (MQ=255)
               ggtTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGAt                                                                                                                                 <  2:643479/150‑1 (MQ=255)
                  tGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCGTGGTGCTGGATATCAAATCGGTAGGGGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGtt                                                                                                                              >  2:77896/1‑150 (MQ=255)
                   gagaCGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCCACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGttt                                                                                                                             >  1:158278/1‑150 (MQ=255)
                     gaCGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTa                                                                                                                           >  1:693783/1‑150 (MQ=255)
                       cGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATg                                                                                                                         <  2:987793/150‑1 (MQ=255)
                             agcatgtGGGTATCCGGGACATGGTTTTCACTGGGGAAGCTGGCGGTTCGGTACTTGCGGAATTGTTTAGATTCCTGGTGCTGGATATCAATTCGGTAGGCGATGAAGTCATTACCGATAAAGGCGAAGTTCTGATGTTTTATGCcgcgc                                                                                                                   <  2:461277/144‑1 (MQ=255)
                              gCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCcgcgcg                                                                                                                  <  1:130295/150‑1 (MQ=255)
                               cAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCCATAAAGCCGAAGTTCTGATGTTTTATGCcgcgcgc                                                                                                                 >  2:298562/1‑150 (MQ=255)
                               cAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCCATAAAGCCGAAGTTCTGATGTTTTATGCcgcgcgc                                                                                                                 >  2:304420/1‑150 (MQ=255)
                                  cTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGtt                                                                                                              <  2:89666/150‑1 (MQ=255)
                                   tGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTc                                                                                                             <  1:1263817/150‑1 (MQ=255)
                                     ggTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCaa                                                                                                           <  2:691573/150‑1 (MQ=255)
                                        aTCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTACCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTg                                                                                                        <  1:355273/150‑1 (MQ=255)
                                         tCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCTGCGCGCGTTCAACTgg                                                                                                       <  2:158278/150‑1 (MQ=255)
                                         tCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTgg                                                                                                       <  2:731231/150‑1 (MQ=255)
                                         tCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTgg                                                                                                       <  1:1296918/150‑1 (MQ=255)
                                          ccGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGt                                                                                                      <  2:1075858/150‑1 (MQ=255)
                                           cgcgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTa                                                                                                     >  1:1108448/1‑150 (MQ=255)
                                           cgcgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTa                                                                                                     <  2:40784/150‑1 (MQ=255)
                                            gcgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCCATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAg                                                                                                    <  1:298562/150‑1 (MQ=255)
                                            gcgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCCATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAg                                                                                                    <  1:304420/150‑1 (MQ=255)
                                             cgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGGCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGa                                                                                                   <  2:448177/150‑1 (MQ=255)
                                             cgACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGGTAAGAAGCCGGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTGCTGATGCTTTATGCCGCGCGCGTGCAACTGGTAGa                                                                                                   <  2:580475/150‑1 (MQ=255)
                                              gACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGaa                                                                                                  >  2:831492/1‑150 (MQ=255)
                                                cATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAAc                                                                                                >  1:89866/1‑150 (MQ=255)
                                                   ggTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGt                                                                                             >  2:670583/1‑150 (MQ=255)
                                                    gTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGtc                                                                                            <  2:1392905/150‑1 (MQ=255)
                                                     ttttCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGtca                                                                                           >  2:975005/1‑150 (MQ=255)
                                                      tttCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGtcat                                                                                          <  2:1416322/150‑1 (MQ=255)
                                                       ttCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGtcatc                                                                                         <  1:496755/150‑1 (MQ=255)
                                                        tCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGtcatca                                                                                        >  2:246355/1‑150 (MQ=255)
                                                         cACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATATAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCaa                                                                                       >  2:938911/1‑150 (MQ=255)
                                                          aCTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCaaa                                                                                      >  1:835633/1‑150 (MQ=255)
                                                           cTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGAGCGCGTTCAACTGGTAGAAACGGTCATCAAAc                                                                                     >  2:1428026/1‑150 (MQ=255)
                                                                   ccTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTg                                                                             <  1:537727/150‑1 (MQ=255)
                                                                    cTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTgg                                                                            <  1:670583/150‑1 (MQ=255)
                                                                     tGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGc                                                                           <  2:54683/150‑1 (MQ=255)
                                                                      ggCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCt                                                                          >  1:1342729/1‑150 (MQ=255)
                                                                           tACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACgg                                                                     <  1:1401389/150‑1 (MQ=255)
                                                                             cGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCa                                                                   >  1:1275722/1‑150 (MQ=255)
                                                                                   ttGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTGCCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTggg                                                             >  2:1354093/1‑150 (MQ=255)
                                                                                   ttGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTggg                                                             <  2:428717/150‑1 (MQ=255)
                                                                                    tGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGt                                                            <  2:1029079/150‑1 (MQ=255)
                                                                                     gCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTg                                                           >  1:406333/1‑150 (MQ=255)
                                                                                         aaaaGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTg                                                       <  2:1108448/150‑1 (MQ=255)
                                                                                           aaGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGc                                                     <  2:131202/150‑1 (MQ=255)
                                                                                             gTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGa                                                   <  1:1088491/150‑1 (MQ=255)
                                                                                                  gaagCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGcc                                              >  2:1432646/1‑150 (MQ=255)
                                                                                                   aagCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCa                                             <  2:1275722/150‑1 (MQ=255)
                                                                                                   aagCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTAGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCa                                             >  1:297280/1‑150 (MQ=255)
                                                                                                   aagCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTAGTAGAAACGGTCAACAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCa                                             >  1:303066/1‑150 (MQ=255)
                                                                                                        tGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAACCCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGAtc                                        >  2:1220389/1‑150 (MQ=255)
                                                                                                               ggATATCATATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGATGTTCTGATGTTTTATGCCGCGCGCGTTCTACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTTACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCAc                                 <  1:1068769/150‑1 (MQ=255)
                                                                                                               ggATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTAGCGATCGCCACGATCTCTCCAc                                 >  2:425178/1‑150 (MQ=255)
                                                                                                                gATATCAAATCGGTAGGCGCTGAAGGCATTACCGATACAGCCGAAGTTCTGATGTTTTATCCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACt                                <  2:1566800/150‑1 (MQ=255)
                                                                                                                gATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACt                                >  1:1174806/1‑150 (MQ=255)
                                                                                                                  tatCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCGCCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCa                              >  1:714908/1‑150 (MQ=255)
                                                                                                                     cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGc                           >  2:1560208/1‑150 (MQ=255)
                                                                                                                     cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGc                           >  1:1594790/1‑150 (MQ=255)
                                                                                                                       aaTCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGt                         >  1:913808/1‑150 (MQ=255)
                                                                                                                          cGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATc                      <  2:89866/150‑1 (MQ=255)
                                                                                                                               ggCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggc                 <  2:1342729/150‑1 (MQ=255)
                                                                                                                                 cGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggcgg               <  2:779074/150‑1 (MQ=255)
                                                                                                                                   aTGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCATCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGgcggcg             <  2:693783/150‑1 (MQ=255)
                                                                                                                                        gTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGt        >  2:99326/1‑150 (MQ=255)
                                                                                                                                          cATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg      <  2:195900/150‑1 (MQ=255)
                                                                                                                                              aCCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTAtt  >  1:281003/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGCGCGTAATGTGGTGGTTGAGACGCTCGAGCAACTGGGTATCCGCGACATGGTTTTCACTCGGGAACCTGGCGGTACGCAACTTGCCGAAAAGTTAAGAAGCCTGGTGCTGGATATCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATT  >  NC_000913/1155177‑1155468

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