Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,482,227 C→A 8.7% intergenic (+26/+28) 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,2270CA 8.7% 301.6 / 13.5 92intergenic (+26/+28)ynbD/azoRputative phosphatase/FMN dependent NADH:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base A (3/5);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

TTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAAG  >  NC_000913/1482079‑1482375
                                                                                                                                                    |                                                                                                                                                    
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: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                                                                                               >  2:1452182/1‑150 (MQ=255)
                                                      gCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGgct                                                                                               >  1:384880/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                                          >  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: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                                     <  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                      >  1:1038990/1‑135 (MQ=255)
                                                                                                                                              gAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTca                      <  2:1038990/135‑1 (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)
                                                                                                                                                  atCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAACGCCGATaaa   >  2:32631/1‑150 (MQ=255)
                                                                                                                                                   tCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAGTTTCACATCGGTTATGCCGATAAAg  >  2:864742/1‑150 (MQ=255)
                                                                                                                                                   tCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGCATTACCGAtgtga  <  2:94173/150‑5 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGTTATGTTACGGACACTGTAAAACCGTTAATGAAGCGATTAGCTATATTCGAGCCAGACGCCCGCAGATTGTGCTGACAGACGAGCACAAAGCGATGCTGAGATTATGGGAAAACAGGTAAGTGGATTGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTTTCGCTGCCATTTCCGGACCGTATGCGATCCCTTCGGCGAAGACAAATTTCACATCGGTAATGCCGATAAAG  >  NC_000913/1482079‑1482375

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