Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,256 T→A 6.2% I272N (ATT→AAT)  cmoB → tRNA U34 carboxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,2560TA 6.2% 356.4 / 10.4 112I272N (ATT→AAT) cmoBtRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (50/55);  new base A (3/4);  total (53/59)
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.71e-01

Only 100 of 114 total aligned reads displayed.

GTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCA  >  NC_000913/1954108‑1954405
                                                                                                                                                    |                                                                                                                                                     
gTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAtt                                                                                                                                                      >  1:1006840/1‑150 (MQ=255)
 tGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATg                                                                                                                                                     <  1:1497747/150‑1 (MQ=255)
 tGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATg                                                                                                                                                     <  1:1491745/150‑1 (MQ=255)
  gCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                    >  2:749972/1‑150 (MQ=255)
  gCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGc                                                                                                                                                    <  2:619171/150‑1 (MQ=255)
   cTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCg                                                                                                                                                   >  2:151163/1‑150 (MQ=255)
    tGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCgg                                                                                                                                                  <  2:42788/150‑1 (MQ=255)
     ggAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                                 <  1:722188/150‑1 (MQ=255)
     ggAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                                 >  2:1365290/1‑150 (MQ=255)
       aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtg                                                                                                                                               <  2:607752/150‑1 (MQ=255)
       aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtg                                                                                                                                               <  1:1182399/150‑1 (MQ=255)
         aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgtg                                                                                                                                             <  1:134720/150‑1 (MQ=255)
          cGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                            >  1:1548153/1‑150 (MQ=255)
          cGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGa                                                                                                                                            <  2:293394/150‑1 (MQ=255)
           gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGTCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAg                                                                                                                                           >  2:1493976/1‑150 (MQ=255)
           gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAtatcta                                                                                                                                           <  1:1556506/150‑6 (MQ=255)
            cTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGc                                                                                                                                          >  1:321089/1‑150 (MQ=255)
             tGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCg                                                                                                                                         >  2:1509609/1‑150 (MQ=255)
             tGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCTGAtatctaca                                                                                                                                         <  2:1606045/150‑8 (MQ=255)
              ggTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCATATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGt                                                                                                                                        >  1:336562/1‑150 (MQ=255)
                ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTa                                                                                                                                      <  2:631610/150‑1 (MQ=255)
                ttATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTCGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTa                                                                                                                                      >  2:283964/1‑150 (MQ=255)
                 tATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTac                                                                                                                                     <  1:1101206/150‑1 (MQ=255)
                  aTTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                                    >  1:800884/1‑150 (MQ=255)
                   ttGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacca                                                                                                                                   <  1:626366/150‑1 (MQ=255)
                   ttGATGGTGATGAGAACACGCTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTACCGCGCTGGCGCTGAAAAACTGGATGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacca                                                                                                                                   >  1:1387051/1‑150 (MQ=255)
                    tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccac                                                                                                                                  >  1:1321559/1‑150 (MQ=255)
                    tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAACGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccac                                                                                                                                  <  1:899226/150‑1 (MQ=255)
                       tGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACcaccac                                                                                                                               >  1:311296/1‑150 (MQ=255)
                       tGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACcaccac                                                                                                                               <  1:1306230/150‑1 (MQ=255)
                        ggTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACcaca                                                                                                                              <  2:1211438/150‑1 (MQ=255)
                             tgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaaga                                                                                                                         >  2:8320/1‑150 (MQ=255)
                             tgaGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaaga                                                                                                                         <  1:1437716/150‑1 (MQ=255)
                               agaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagc                                                                                                                       <  2:54114/150‑1 (MQ=255)
                                gaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                                      >  1:742979/1‑150 (MQ=255)
                                  acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcagc                                                                                                                    >  2:361269/1‑150 (MQ=255)
                                  acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcagc                                                                                                                    >  2:254446/1‑150 (MQ=255)
                                   cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCg                                                                                                                   <  2:581688/150‑1 (MQ=255)
                                   cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCg                                                                                                                   <  2:391457/150‑1 (MQ=255)
                                     cGGTGCTGGTACCACGCGAGCGGTACGCGCAAATGCGTAATGTCTCTTTCATCCCATCCGCGCTGGCGCTGAAAACCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCCTTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGAc                                                                                                                 <  2:174746/150‑1 (MQ=255)
                                       gTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGc                                                                                                               >  1:176886/1‑150 (MQ=255)
                                         gCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCAc                                                                                                             >  1:1461289/1‑150 (MQ=255)
                                             gTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                         >  2:324140/1‑150 (MQ=255)
                                              tACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAat                                                                                                        <  1:283964/150‑1 (MQ=255)
                                                 cAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTAGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgga                                                                                                     >  1:403790/1‑150 (MQ=255)
                                                        tCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTAc                                                                                              <  2:1043301/150‑1 (MQ=255)
                                                          gTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATAGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCg                                                                                            >  1:1504623/1‑150 (MQ=255)
                                                             aCGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAAt                                                                                         >  2:1061722/1‑150 (MQ=255)
                                                              cgcgCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATc                                                                                        >  1:82466/1‑150 (MQ=255)
                                                                   aaaTGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTgg                                                                                   <  2:318845/150‑1 (MQ=255)
                                                                              gTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCtt                                                                        <  1:652201/150‑2 (MQ=255)
                                                                                    ttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGAcccc                                                                  >  2:379479/1‑149 (MQ=255)
                                                                                       aTCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACAt                                                               >  1:742105/1‑150 (MQ=255)
                                                                                       aTCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACAt                                                               >  2:1050311/1‑150 (MQ=255)
                                                                                          ccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGAt                                                            >  1:1531663/1‑150 (MQ=255)
                                                                                           cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATc                                                           >  1:1332085/1‑150 (MQ=255)
                                                                                               cgcgCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCggg                                                       <  1:324140/150‑1 (MQ=255)
                                                                                                 cgctggcgctgAAAAACTGTCTTAAGAATTGTGTTTTTTTATATATACTAATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCTATTTCCTCGACCCACATGATCCGGGTa                                                     <  1:493468/150‑1 (MQ=255)
                                                                                                 cgctggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTa                                                     <  2:82466/150‑1 (MQ=255)
                                                                                                 cgctggcgctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTa                                                     <  2:1461289/150‑1 (MQ=255)
                                                                                                         ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTg                                             <  2:1429834/150‑1 (MQ=255)
                                                                                                         ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTg                                             <  1:913726/150‑1 (MQ=255)
                                                                                                           gAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGa                                           >  2:1153347/1‑150 (MQ=255)
                                                                                                             aaaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAg                                         <  2:800124/150‑1 (MQ=255)
                                                                                                              aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAgg                                        <  2:336562/150‑1 (MQ=255)
                                                                                                              aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAgg                                        <  2:742105/150‑1 (MQ=255)
                                                                                                              aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAgg                                        <  2:5715/150‑1 (MQ=255)
                                                                                                               aaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGt                                       <  1:1453914/150‑1 (MQ=255)
                                                                                                                aCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                                      >  1:601199/1‑150 (MQ=255)
                                                                                                                aCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                                      >  2:766407/1‑150 (MQ=255)
                                                                                                                      tGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCgg                                >  1:199233/1‑150 (MQ=255)
                                                                                                                      tGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCgg                                >  2:463868/1‑150 (MQ=255)
                                                                                                                       gaagaaGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                               >  1:1188564/1‑150 (MQ=255)
                                                                                                                          gaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCAcc                            >  2:568666/1‑150 (MQ=255)
                                                                                                                          gaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCAcc                            <  1:749972/150‑1 (MQ=255)
                                                                                                                          gaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTATACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCAcc                            <  2:311296/150‑1 (MQ=255)
                                                                                                                             gtgtGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTaa                         >  1:901342/1‑150 (MQ=255)
                                                                                                                              tgtgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTaaa                        <  1:524509/150‑1 (MQ=255)
                                                                                                                                tgGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcg                      >  1:1405750/1‑150 (MQ=255)
                                                                                                                                 ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgc                     >  2:315441/1‑150 (MQ=255)
                                                                                                                                    tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATTGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAc                       >  1:271971/1‑145 (MQ=255)
                                                                                                                                    tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATTGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAc                       <  2:271971/145‑1 (MQ=255)
                                                                                                                                     ttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCTATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcgcgt                 <  2:1387051/150‑2 (MQ=255)
                                                                                                                                      tGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGt                <  1:1061722/150‑1 (MQ=255)
                                                                                                                                        ctcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc              >  1:24457/6‑150 (MQ=255)
                                                                                                                                        ctcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc              >  1:228502/6‑150 (MQ=255)
                                                                                                                                        ctcacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGc              >  1:1105432/6‑150 (MQ=255)
                                                                                                                                         aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCt             <  2:1006840/150‑1 (MQ=255)
                                                                                                                                         aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCt             <  1:1153380/150‑1 (MQ=255)
                                                                                                                                          gATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTg            <  2:18442/150‑1 (MQ=255)
                                                                                                                                            tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGAt          <  1:568666/150‑1 (MQ=255)
                                                                                                                                             atCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGAtt         <  1:528807/150‑1 (MQ=255)
                                                                                                                                             atCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGAtt         >  2:652236/1‑150 (MQ=255)
                                                                                                                                              tCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTg        >  1:1049097/1‑150 (MQ=255)
                                                                                                                                              tCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTg        <  1:361269/150‑1 (MQ=255)
                                                                                                                                                cGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcg      >  1:956606/1‑150 (MQ=255)
                                                                                                                                                 gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcgc     <  1:1365290/150‑1 (MQ=255)
                                                                                                                                                  cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcgcg    <  1:1493976/150‑1 (MQ=255)
                                                                                                                                                  cATAGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTgcgcg    <  2:1504623/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCa  >  2:1369659/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGCTGGAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATTGCGCGCA  >  NC_000913/1954108‑1954405

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