Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,224 T→G 8.9% intergenic (+23/+31) ynbD → / ← azoR putative phosphatase/FMN dependent NADH:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,482,2240TG 8.9% 253.3 / 14.5 88intergenic (+23/+31)ynbD/azoRputative phosphatase/FMN dependent NADH:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (43/37);  new base G (3/5);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATA  >  NC_000913/1482076‑1482372
                                                                                                                                                    |                                                                                                                                                    
tGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAata                                                                                                                                                     <  1:280952/150‑1 (MQ=255)
gggTTGTTATTTTACGGACACTGTAAAACCGTTTCTGTAGCGATTAGCTATATTCGAGGCAGAGGGCCGCAGATTGTGCTCACCGACGATCACATAGCGATGCTGTGATTATGTGATTACAGTTAAGTGGATTGAGATGTGGTCTGAata                                                                                                                                                     <  2:944422/149‑1 (MQ=255)
   ttgttATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCt                                                                                                                                                  <  1:374802/150‑1 (MQ=255)
    tgttatgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTa                                                                                                                                                 <  2:286777/150‑1 (MQ=255)
     gttatgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTAc                                                                                                                                                >  1:1192172/1‑150 (MQ=255)
     gttatgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTAc                                                                                                                                                >  1:1259337/1‑150 (MQ=255)
       tatgttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAg                                                                                                                                              <  2:1101301/150‑1 (MQ=255)
          gttaCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTcc                                                                                                                                           >  2:1424178/1‑150 (MQ=255)
             aCGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCcaca                                                                                                                                        >  2:941859/1‑150 (MQ=255)
              cGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACAt                                                                                                                                       >  1:1441855/1‑150 (MQ=255)
                gACACTGTAAAACCGTTAATGAAGCTATTAGCTATATTCGAGCCAGACGCCCGCATTTTGTGCTGACAGACGAGCACAAAGCGATTCTGACATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATct                                                                                                                                     <  2:1281287/150‑2 (MQ=255)
                gACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCa                                                                                                                                     <  2:1441855/150‑1 (MQ=255)
                gACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCa                                                                                                                                     <  2:1165267/150‑1 (MQ=255)
                gACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCa                                                                                                                                     <  1:1430552/150‑1 (MQ=255)
                 acacTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATctc                                                                                                                                    <  1:1102539/150‑3 (MQ=255)
                  cacTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATctca                                                                                                                                   <  2:1271081/150‑4 (MQ=255)
                    cTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGAATGAATATCTACAGTCAACATCAAgcc                                                                                                                                 >  2:793024/1‑148 (MQ=255)
                     tGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGAcc                                                                                                                                >  2:263848/1‑150 (MQ=255)
                      gTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCg                                                                                                                               >  1:146970/1‑150 (MQ=255)
                        aaaaCCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCgtg                                                                                                                             <  2:1191318/150‑1 (MQ=255)
                             cGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATc                                                                                                                                      >  2:645095/1‑136 (MQ=255)
                             cGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATc                                                                                                                                      <  1:645095/136‑1 (MQ=255)
                             cGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCgg                                                                                                                        <  2:91738/150‑1 (MQ=255)
                              gTTAATGTAGCTATTAGCTATATTCTAGCCAGACGCCCGCAGTTTGTGCTGAGAGACGAGCTCAAAGCGATGCTGAGTTTATGGGAAATCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGt                                                                                                                       <  2:439565/150‑1 (MQ=255)
                              gTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGt                                                                                                                       >  2:1427946/1‑150 (MQ=255)
                                tAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACCTCAAGACCGTGTCCGGTTa                                                                                                                     >  2:727885/1‑150 (MQ=255)
                                    gAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCa                                                                                                                 >  2:1423353/1‑150 (MQ=255)
                                      aGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGTAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGa                                                                                                               <  2:181366/150‑1 (MQ=255)
                                       gCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGGGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGaa                                                                                                              <  2:1185562/150‑1 (MQ=255)
                                         gATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAAc                                                                                                            <  1:800006/150‑1 (MQ=255)
                                             aGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTAc                                                                                                                                                >  1:949530/1‑110 (MQ=255)
                                             aGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTAc                                                                                                                                                <  2:949530/110‑1 (MQ=255)
                                             aGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATg                                                                                                        >  1:747643/1‑150 (MQ=255)
                                              gCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGc                                                                                                       >  2:1284722/1‑150 (MQ=255)
                                              gCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGc                                                                                                       >  2:1039226/1‑150 (MQ=255)
                                              gCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACGGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGc                                                                                                       >  1:35368/1‑150 (MQ=255)
                                               cTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAATACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATggt                                                                                                      >  2:347541/1‑148 (MQ=255)
                                                tataTTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCCCAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTg                                                                                                     >  2:934979/1‑150 (MQ=255)
                                                      cGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATg                                                                                               >  2:475488/1‑150 (MQ=255)
                                                       gAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATgg                                                                                              >  2:1130154/1‑150 (MQ=255)
                                                         gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                            >  1:706583/1‑150 (MQ=255)
                                                         gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                            >  1:501561/1‑150 (MQ=255)
                                                         gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                            >  2:1452182/1‑150 (MQ=255)
                                                         gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                            >  1:384880/1‑150 (MQ=255)
                                                           cAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgc                                                                                          >  2:5571/1‑150 (MQ=255)
                                                           cAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGTTATGCAGAAACAATGCTGTCGATGgctgct                                                                                         >  1:1042901/1‑150 (MQ=255)
                                                            aGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgct                                                                                         <  2:1047657/150‑1 (MQ=255)
                                                                    gCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGta                                                                                 >  2:1439603/1‑149 (MQ=255)
                                                                    gCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCa                                                                                <  2:1042901/150‑1 (MQ=255)
                                                                      aGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAg                                                                               <  1:654801/150‑1 (MQ=255)
                                                                              cTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGAGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCt                                                                       >  1:1245022/1‑150 (MQ=255)
                                                                               tGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCtt                                                                      >  1:900322/1‑150 (MQ=255)
                                                                                 aCAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTc                                                                    <  1:784127/150‑1 (MQ=255)
                                                                                 aCAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATGGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGGTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGGCAGACTGTGCTTTc                                                                    <  1:784129/150‑1 (MQ=255)
                                                                                       gAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGTGTCAGACTGTGCTTTCGCTGcc                                                              <  1:475488/150‑1 (MQ=255)
                                                                                            caAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTc                                                         <  1:793024/150‑1 (MQ=255)
                                                                                              aaGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCg                                                       <  1:1284722/150‑1 (MQ=255)
                                                                                               aGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCgg                                                      <  1:1509628/150‑1 (MQ=255)
                                                                                                        gagaTTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGc                                             <  1:1423353/150‑1 (MQ=255)
                                                                                                            ttATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATc                                         <  1:263848/150‑1 (MQ=255)
                                                                                                              aTGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATccc                                       <  1:1130154/150‑1 (MQ=255)
                                                                                                              aTGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATccc                                       >  2:388446/1‑150 (MQ=255)
                                                                                                               tGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCt                                      >  1:582556/1‑150 (MQ=255)
                                                                                                               tGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCt                                      >  2:386108/1‑150 (MQ=255)
                                                                                                                gggAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCtt                                     >  1:1381728/1‑150 (MQ=255)
                                                                                                                     aaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCg                                <  2:146970/150‑1 (MQ=255)
                                                                                                                     aaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCg                                <  2:1259337/150‑1 (MQ=255)
                                                                                                                        aGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAg                             <  2:747643/150‑1 (MQ=255)
                                                                                                                        aGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAg                             <  1:1460567/150‑1 (MQ=255)
                                                                                                                        aGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAg                             <  2:1192172/150‑1 (MQ=255)
                                                                                                                         ggTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGa                            >  1:378662/1‑150 (MQ=255)
                                                                                                                                ggATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAAttt                     >  1:1289999/1‑150 (MQ=255)
                                                                                                                                  aTTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTca                   <  1:712357/150‑1 (MQ=255)
                                                                                                                                   ttGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGTTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTcac                  <  1:316919/150‑1 (MQ=255)
                                                                                                                                   ttGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTcac                  <  1:941859/150‑1 (MQ=255)
                                                                                                                                     gagaTGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACAt                >  1:1406905/1‑150 (MQ=255)
                                                                                                                                      agaTGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATc               >  2:969161/1‑150 (MQ=255)
                                                                                                                                         tgtgGACTGTAGATATTCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGt            >  2:288306/1‑150 (MQ=255)
                                                                                                                                         tgtgGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGt            <  2:1289999/150‑1 (MQ=255)
                                                                                                                                          gtgGACTGTAGATATTCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCgg                                  >  1:583372/1‑127 (MQ=255)
                                                                                                                                          gtgGACTGTAGATATTCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCgg                                  <  2:583372/127‑1 (MQ=255)
                                                                                                                                           tgGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTaa          >  1:1198039/1‑150 (MQ=255)
                                                                                                                                              aCTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGc       <  1:1452182/150‑1 (MQ=255)
                                                                                                                                               cTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGcc      >  1:797147/1‑150 (MQ=255)
                                                                                                                                                 gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTca                   <  2:1038990/135‑1 (MQ=255)
                                                                                                                                                 gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTca                   >  1:1038990/1‑135 (MQ=255)
                                                                                                                                                 gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa    >  2:558161/1‑150 (MQ=255)
                                                                                                                                                 gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa    <  2:384880/150‑1 (MQ=255)
                                                                                                                                                 gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa    >  2:100008/1‑150 (MQ=255)
                                                                                                                                                  aaTATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGAt   >  2:1309680/1‑150 (MQ=255)
                                                                                                                                                   atatCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTTTGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATa  <  2:365640/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATA  >  NC_000913/1482076‑1482372

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: