Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,222 A→T 9.0% intergenic (+21/+33) 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,2220AT 9.0% 265.3 / 13.5 89intergenic (+21/+33)ynbD/azoRputative phosphatase/FMN dependent NADH:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (42/39);  new base T (3/5);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCATGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGAT  >  NC_000913/1482073‑1482371
                                                                                                                                                     |                                                                                                                                                     
gcATGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGa                                                                                                                                                       >  1:1330484/1‑150 (MQ=255)
   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                                                                                                                                               <  2:949530/110‑1 (MQ=255)
                                                aGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTAc                                                                                                                                               >  1:949530/1‑110 (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:384880/1‑150 (MQ=255)
                                                            gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                           >  1:706583/1‑150 (MQ=255)
                                                            gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                           >  2:1452182/1‑150 (MQ=255)
                                                            gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                           >  1:501561/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                                      >  2:388446/1‑150 (MQ=255)
                                                                                                                 aTGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATccc                                      <  1:1130154/150‑1 (MQ=255)
                                                                                                                  tGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCt                                     >  2:386108/1‑150 (MQ=255)
                                                                                                                  tGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCt                                     >  1:582556/1‑150 (MQ=255)
                                                                                                                   gggAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCtt                                    >  1:1381728/1‑150 (MQ=255)
                                                                                                                        aaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCg                               <  2:1259337/150‑1 (MQ=255)
                                                                                                                        aaCAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCg                               <  2:146970/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                                 <  2:583372/127‑1 (MQ=255)
                                                                                                                                             gtgGACTGTAGATATTCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCgg                                 >  1:583372/1‑127 (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:384880/150‑1 (MQ=255)
                                                                                                                                                    gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa   >  2:558161/1‑150 (MQ=255)
                                                                                                                                                    gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa   >  2:100008/1‑150 (MQ=255)
                                                                                                                                                     aaTATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGAt  >  2:1309680/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCATGGTTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGAT  >  NC_000913/1482073‑1482371

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: