Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,229 A→T 8.9% intergenic (+28/+26) 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,2290AT 8.9% 261.2 / 13.5 90intergenic (+28/+26)ynbD/azoRputative phosphatase/FMN dependent NADH:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (45/37);  new base T (3/5);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

TGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAAGC  >  NC_000913/1482080‑1482376
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                          <  1:645095/136‑1 (MQ=255)
                         cGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGTAGATATTCAGTCCACATc                                                                                                                                          >  2:645095/1‑136 (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:1039226/1‑150 (MQ=255)
                                          gCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGc                                                                                                           >  2:1284722/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                                                                                                >  1:384880/1‑150 (MQ=255)
                                                     gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                                >  2:1452182/1‑150 (MQ=255)
                                                       cAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgc                                                                                              >  2:5571/1‑150 (MQ=255)
                                                       cAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCC‑GTTATGCAGAAACAATGCTGTCGATGgctgct                                                                                             >  1:1042901/1‑150 (MQ=255)
                                                        aGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgctgct                                                                                             <  2:1047657/150‑1 (MQ=255)
                                                                gCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGta                                                                                     >  2:1439603/1‑149 (MQ=255)
                                                                gCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCC‑GTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCa                                                                                    <  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                                          >  2:386108/1‑150 (MQ=255)
                                                                                                           tGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCt                                          >  1:582556/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                                 <  2:1192172/150‑1 (MQ=255)
                                                                                                                    aGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAg                                 <  1:1460567/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:100008/1‑150 (MQ=255)
                                                                                                                                             gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGa        <  2:384880/150‑1 (MQ=255)
                                                                                                                                              aaTATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGAt       >  2:1309680/1‑150 (MQ=255)
                                                                                                                                               atatCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTTTGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATa      <  2:365640/150‑1 (MQ=255)
                                                                                                                                                 atCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAACGCCGATaaa    >  2:32631/1‑150 (MQ=255)
                                                                                                                                                  tCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAGTTTCACATCGGTTATGCCGATAAAg   >  2:864742/1‑150 (MQ=255)
                                                                                                                                                  tCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGCATTACCGAtgtga   <  2:94173/150‑5 (MQ=255)
                                                                                                                                                   cgaCAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAAAGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCCCTGCCATTTCCGGCCCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATGGGGAATGCCGATAAAGc  >  2:466227/3‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAAGC  >  NC_000913/1482080‑1482376

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: ‑