Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,862,057 T→C Y210H (TAT→CAT)  ygbJ → putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,862,0570TC99.0% 401.8 / ‑3.9 104Y210H (TAT→CAT) ygbJputative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (51/52);  total (52/52)
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.97e-01

Only 100 of 106 total aligned reads displayed.

GCCGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTG  >  NC_000913/2861910‑2862203
                                                                                                                                                   |                                                                                                                                                  
gCCGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATg                                                                                                                                                 >  4:343489/1‑151 (MQ=255)
 ccGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTTTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGa                                                                                                                                                <  3:241574/151‑1 (MQ=255)
 ccGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGa                                                                                                                                                <  2:196530/151‑1 (MQ=255)
    gAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgt                                                                                                                                             <  2:1028233/151‑1 (MQ=255)
    gAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgt                                                                                                                                             >  2:1082971/1‑151 (MQ=255)
     aaaaGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgtc                                                                                                                                            >  4:199079/1‑151 (MQ=255)
      aaaGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgtcg                                                                                                                                           >  2:638936/1‑151 (MQ=255)
      aaaGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgtcg                                                                                                                                           >  3:250429/1‑151 (MQ=255)
        aGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTg                                                                                                                                         >  2:888470/1‑151 (MQ=255)
        aGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTg                                                                                                                                         <  1:459999/151‑1 (MQ=255)
         gTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTg                                                                                                                                         >  1:351359/1‑150 (MQ=255)
         gTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGa                                                                                                                                        <  3:114126/151‑1 (MQ=255)
          tttATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGAc                                                                                                                                       >  1:1201522/1‑151 (MQ=255)
            tATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGAcc                                                                                                                                      <  4:435553/150‑1 (MQ=255)
            tATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGAcc                                                                                                                                      >  2:341231/1‑150 (MQ=255)
             aTCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCa                                                                                                                                     >  2:1043654/1‑150 (MQ=255)
                    aGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATgccgcc                                                                                                                             <  4:413432/151‑1 (MQ=255)
                     ggTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGccgccg                                                                                                                            <  4:415441/151‑1 (MQ=255)
                      gTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCgg                                                                                                                           <  1:1158345/151‑1 (MQ=255)
                        gCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGaa                                                                                                                         <  2:114449/151‑1 (MQ=255)
                        gCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGa                                                                                                                          <  4:263613/150‑1 (MQ=255)
                         cAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGaaa                                                                                                                        <  2:1204709/151‑1 (MQ=255)
                          aGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCGTGACGTCGTGACCAATGCCGCCGGAAAt                                                                                                                       >  1:1219992/1‑151 (MQ=255)
                          aGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAAt                                                                                                                       >  1:413152/1‑151 (MQ=255)
                             aCCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTcc                                                                                                                    >  1:34421/1‑151 (MQ=255)
                               tgggACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTg                                                                                                                  >  2:1263569/2‑151 (MQ=255)
                                 ggACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTg                                                                                                                  >  2:59748/1‑149 (MQ=255)
                                  gagtAGGTTCGACCGTAAAAATTATTCACCAGGTATTAGCGGGAGTACATATTGCTGCCTGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACTAACGCCGCCGGAAATTCCTGGAt                                                                                                               <  2:590351/148‑1 (MQ=255)
                                    cTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGt                                                                                                             >  3:21829/1‑151 (MQ=255)
                                        gTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCg                                                                                                          >  2:241651/1‑150 (MQ=255)
                                        gTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGa                                                                                                         <  2:944123/151‑1 (MQ=255)
                                          tCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGaa                                                                                                        >  4:49414/1‑150 (MQ=255)
                                          tCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGaa                                                                                                        >  1:419495/1‑150 (MQ=255)
                                           cGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGaaaa                                                                                                      <  1:548296/151‑1 (MQ=255)
                                            gACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAAc                                                                                                     >  3:315726/1‑151 (MQ=255)
                                            gACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAAc                                                                                                     >  2:37047/1‑151 (MQ=255)
                                             aCCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAAcc                                                                                                    <  4:108613/151‑1 (MQ=255)
                                               cGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCgg                                                                                                  >  1:277303/1‑151 (MQ=255)
                                                gTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGa                                                                                                 <  2:1171581/151‑1 (MQ=255)
                                                  aaaaaTTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGAt                                                                                                <  2:137861/150‑1 (MQ=255)
                                                    aaaTTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCg                                                                                             >  1:654304/1‑151 (MQ=255)
                                                    aaaTTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCg                                                                                             <  4:122928/151‑1 (MQ=255)
                                                     aaTTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGt                                                                                            >  2:1179155/1‑151 (MQ=255)
                                                        tattCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCAt                                                                                         >  2:458788/1‑151 (MQ=255)
                                                         attCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCAtg                                                                                        >  1:6526/1‑151 (MQ=255)
                                                         attCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCAtg                                                                                        <  4:128902/151‑1 (MQ=255)
                                                           tCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCAtgt                                                                                       <  3:371183/150‑1 (MQ=255)
                                                             aCCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATgtggt                                                                                    >  1:1026936/1‑151 (MQ=255)
                                                              ccAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCAtgtg                                                                                      >  1:459031/1‑148 (MQ=255)
                                                              ccAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATgtggtg                                                                                   >  1:714160/1‑151 (MQ=255)
                                                               cAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGtggtgg                                                                                  <  4:349668/151‑1 (MQ=255)
                                                                  ttgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATg                                                                               >  1:784064/1‑151 (MQ=255)
                                                                   tgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATg                                                                               <  1:939695/150‑1 (MQ=255)
                                                                    gttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGc                                                                             >  4:348622/1‑151 (MQ=255)
                                                                      tAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCg                                                                            <  1:250880/150‑1 (MQ=255)
                                                                          gggCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAATCCGGATGCGTCATGTGGTGGATGGCGATTac                                                                       <  1:514630/151‑1 (MQ=255)
                                                                          gggCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTac                                                                       >  2:837414/1‑151 (MQ=255)
                                                                           ggCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTac                                                                       >  2:1080117/1‑150 (MQ=255)
                                                                            gCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTacac                                                                     <  1:163910/151‑1 (MQ=255)
                                                                             cGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcc                                                                    >  2:512940/1‑151 (MQ=255)
                                                                               tACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcccc                                                                  >  2:908841/1‑151 (MQ=255)
                                                                                aCATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcccc                                                                  <  1:297388/150‑1 (MQ=255)
                                                                                aCATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcccc                                                                  <  1:1126322/150‑1 (MQ=255)
                                                                                aCATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcccc                                                                  >  2:371274/1‑150 (MQ=255)
                                                                                 cATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGc                                                                >  1:365583/1‑151 (MQ=255)
                                                                                      tgctgcCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTc                                                            <  2:1014565/150‑1 (MQ=255)
                                                                                      tgctgcCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTc                                                            >  1:447489/1‑150 (MQ=255)
                                                                                          gcCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCg                                                       >  3:86397/1‑151 (MQ=255)
                                                                                             ggAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCg                                                    <  2:436779/151‑1 (MQ=255)
                                                                                              gAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCGATGCCGCCGGAAATTCCTGGATGTTCGAAGACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGa                                                   <  3:402192/151‑1 (MQ=255)
                                                                                                gCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGata                                                 <  1:1197749/151‑1 (MQ=255)
                                                                                                gCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGACGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGata                                                 <  1:202095/151‑1 (MQ=255)
                                                                                                 ccGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGatat                                                >  3:448095/1‑151 (MQ=255)
                                                                                                   gAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATAttt                                              <  2:886267/151‑1 (MQ=255)
                                                                                                      gCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATAttttt                                            >  2:684734/1‑150 (MQ=255)
                                                                                                        gATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATCTTTTTGtt                                         >  3:123506/1‑151 (MQ=255)
                                                                                                        gATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGtt                                         >  2:368531/1‑151 (MQ=255)
                                                                                                         aTGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTa                                        <  1:1065144/151‑1 (MQ=255)
                                                                                                           ggCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTa                                        >  1:1119133/1‑149 (MQ=255)
                                                                                                             cACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGa                                    <  1:1043612/151‑1 (MQ=255)
                                                                                                             cACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGa                                    <  1:1169722/151‑1 (MQ=255)
                                                                                                              aCTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGAt                                   <  2:413081/151‑1 (MQ=255)
                                                                                                              aCTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGAt                                   <  3:129257/151‑1 (MQ=255)
                                                                                                                    aGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGt                             <  4:225721/151‑1 (MQ=255)
                                                                                                                       ccGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTg                          <  4:365042/151‑1 (MQ=255)
                                                                                                                        cGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTgg                         >  2:888692/1‑151 (MQ=255)
                                                                                                                        cGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTgg                         <  1:1083175/151‑1 (MQ=255)
                                                                                                                        cGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTgg                         <  2:288925/151‑1 (MQ=255)
                                                                                                                         gTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGt                        >  2:176418/1‑151 (MQ=255)
                                                                                                                         gTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGt                        >  1:593642/1‑151 (MQ=255)
                                                                                                                          tGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGtt                       >  3:357104/1‑151 (MQ=255)
                                                                                                                            cGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGc                     >  1:1069185/1‑151 (MQ=255)
                                                                                                                             ggggATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGcc                    <  2:463732/151‑1 (MQ=255)
                                                                                                                                   cccGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATAc               <  2:391217/150‑1 (MQ=255)
                                                                                                                                    ccGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAg             <  2:367957/151‑1 (MQ=255)
                                                                                                                                     cGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGc            <  1:498144/151‑1 (MQ=255)
                                                                                                                                      gCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGACATTCCTGGATGTTCGAAAACCGGATGCGCCATGTGGCGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGcc           >  2:1155584/1‑151 (MQ=255)
                                                                                                                                         ggATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCaaa        <  1:214349/151‑1 (MQ=255)
                                                                                                                                         ggATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCaa         <  4:225880/150‑1 (MQ=255)
                                                                                                                                               gatgcATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGTATTCAGCCGTCGATATTTTTGTTAAGGAGCTTGGTCTGGTTGCCGATACAGCCAAAGCCCt   <  1:730928/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCTG  >  NC_000913/2861910‑2862203

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