Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,577,595 C→A 6.9% V9F (GTC→TTC)  yddB ← putative TonB‑dependent outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,577,5950CA 6.9% 260.6 / 2.2 93V9F (GTC→TTC) yddBputative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (39/46);  new base A (4/4);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGT  >  NC_000913/1577448‑1577744
                                                                                                                                                   |                                                                                                                                                     
gtAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGCCCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGc                                                                                                                                                     >  2:1045967/1‑150 (MQ=255)
gtAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGc                                                                                                                                                     >  2:1316538/1‑150 (MQ=255)
    ccGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAgg                                                                                                                                                 <  2:606023/149‑1 (MQ=255)
    aCGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAgg                                                                                                                                                 <  1:525462/150‑1 (MQ=255)
       cTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAAt                                                                                                                                              >  2:587331/1‑150 (MQ=255)
       cTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAAt                                                                                                                                              >  1:1092791/1‑150 (MQ=255)
         cGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACATGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGTCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATaa                                                                                                                                            <  2:973502/150‑1 (MQ=255)
         cGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATaa                                                                                                                                            >  2:675274/1‑150 (MQ=255)
          gATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAg                                                                                                                                           <  1:638333/150‑1 (MQ=255)
          gATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAg                                                                                                                                           <  2:435074/150‑1 (MQ=255)
              cTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCAGGAATAAGAACTc                                                                                                                                      >  2:101241/1‑150 (MQ=255)
                   aCTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCtt                                                                                                                                  >  1:345459/1‑150 (MQ=255)
                      ggTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCat                                                                                                                               >  2:279230/1‑150 (MQ=255)
                       gtgtAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatt                                                                                                                              >  2:763697/1‑150 (MQ=255)
                       gtgtAATGCGTGTGGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatt                                                                                                                              <  1:359570/150‑1 (MQ=255)
                        tgtAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGCACTCGCTTCattt                                                                                                                             >  2:1258107/1‑150 (MQ=255)
                          tAATGCGTGTTGCCGTTTACTGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATTCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatttat                                                                                                                           <  2:805193/150‑1 (MQ=255)
                          tAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatttat                                                                                                                           >  2:84275/1‑150 (MQ=255)
                            aTGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatttattt                                                                                                                         <  2:414166/150‑1 (MQ=255)
                            aTGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCatttattt                                                                                                                         <  1:1103358/150‑1 (MQ=255)
                                gtgtTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTatat                                                                                                                     <  2:777835/150‑1 (MQ=255)
                                 tgtTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATAtc                                                                                                                    <  1:1249690/150‑1 (MQ=255)
                                     gCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                                                                                                                >  1:569521/1‑150 (MQ=255)
                                     gCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                                                                                                                >  2:440803/1‑150 (MQ=255)
                                     gCCGATTACGGGTACAGGGACAGTAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCg                                                                                                                <  1:495090/150‑1 (MQ=255)
                                      ccGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGt                                                                                                               >  1:516825/1‑150 (MQ=255)
                                           tACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTAccgc                                                                                                          >  2:95612/1‑150 (MQ=255)
                                           tACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTAcagc                                                                                                          <  2:473288/150‑1 (MQ=255)
                                           tACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGAAGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTAcagc                                                                                                          >  2:721293/1‑150 (MQ=255)
                                           tACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGAAGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTAcagc                                                                                                          >  2:726052/1‑150 (MQ=255)
                                             cGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTAcagcag                                                                                                        <  2:1231013/150‑1 (MQ=255)
                                             cGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGAAGACAGGAATAAGAACTCGCTTCATTTATTTATATATCCGTAcagcag                                                                                                        >  2:853654/1‑150 (MQ=255)
                                                gTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGc                                                                                                     >  2:306013/1‑150 (MQ=255)
                                                   cAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCttt                                                                                                  <  2:52407/150‑1 (MQ=255)
                                                   cAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCttt                                                                                                  >  1:1252529/1‑150 (MQ=255)
                                                    aGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCtttt                                                                                                 >  1:540357/1‑150 (MQ=255)
                                                    aGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCtttt                                                                                                 <  2:22927/150‑1 (MQ=255)
                                                     gggACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGAAGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTa                                                                                                >  2:985265/1‑150 (MQ=255)
                                                      ggACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTAt                                                                                               <  1:1234765/150‑1 (MQ=255)
                                                      ggACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTAt                                                                                               <  2:650536/150‑1 (MQ=255)
                                                        aCAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTg                                                                                             <  2:550658/150‑1 (MQ=255)
                                                          aGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGtt                                                                                           <  2:978246/150‑1 (MQ=255)
                                                           ggAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGttt                                                                                          <  1:500025/150‑1 (MQ=255)
                                                            gAGCATAGACCGTCATCTCTTCAAAAAAATTCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGtttt                                                                                         <  1:985265/150‑1 (MQ=255)
                                                            gAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGtttt                                                                                         <  2:564679/150‑1 (MQ=255)
                                                               cATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTtat                                                                                      <  1:864518/150‑1 (MQ=255)
                                                                aTAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTtata                                                                                     <  2:1043384/150‑1 (MQ=255)
                                                                 tAGACCTTCATCTCTTCAAAAAAATGATTGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGACACATAAAATGACGCCAGTAATAAGTACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATaa                                                                                    <  1:266438/150‑1 (MQ=255)
                                                                   gACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATaaaa                                                                                  >  1:991109/1‑150 (MQ=255)
                                                                    aCCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAAc                                                                                 <  2:1267915/150‑1 (MQ=255)
                                                                     ccGTCATCTCTTCAAAAAAATGCAAGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGCAATAAGAACTCGCTTCATTTATTTATATATCCGAACAGCAGGCCTTTTATTGTTTTTATAAAAca                                                                                >  1:1042553/1‑149 (MQ=255)
                                                                         cATCACTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCgcgc                                                                            <  1:379896/150‑1 (MQ=255)
                                                                          aTCTCTTCATAATAATGGATGTACATGTTTTTGTCATCGACGTCCTTCGCCTCATCATAGCCACATAAATTGAATCGAGGAATAAGAAGTAGCTTCATTTATTTATATCTCCGTACAGCATGGCTTTTATTGTTTTTATAAAACCGCGCt                                                                           <  2:266072/150‑1 (MQ=255)
                                                                          aTCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCt                                                                           >  1:1237718/1‑150 (MQ=255)
                                                                          aTCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCt                                                                           <  1:1338596/150‑1 (MQ=255)
                                                                               ttCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGt                                                                      >  2:444679/1‑150 (MQ=255)
                                                                                  aaaaaaaTGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTcac                                                                   >  2:1056616/1‑150 (MQ=255)
                                                                                        aTGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAatat                                                             <  1:1077551/150‑1 (MQ=255)
                                                                                           cATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATatcat                                                          >  1:867225/1‑150 (MQ=255)
                                                                                            aTGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATatcatc                                                         <  1:864980/150‑1 (MQ=255)
                                                                                              gTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCgg                                                       <  1:1353801/150‑1 (MQ=255)
                                                                                               tACATGTTTTTGTCATCGACGGCCTGCGCCCCATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGc                                                      >  1:391459/1‑150 (MQ=255)
                                                                                                aCATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACATCAGGGCTTTTATTGTTTTTATAAACCAGCGCTAATGTCACAAATATCATCGGcc                                                     >  1:673807/1‑150 (MQ=255)
                                                                                                aCATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGcc                                                     <  2:749977/150‑1 (MQ=255)
                                                                                                 aaTGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCa                                                    >  1:1344964/2‑150 (MQ=255)
                                                                                                  aTGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGTACTCGCTTCATTTATTTATATCTCCGTACATCAGGGCTTTTATTGTTTTTATAAAACCGCTCTAATGTCACAAATATCATCGGCCAg                                                   <  2:403513/150‑1 (MQ=255)
                                                                                                   tGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAgg                                                  <  1:494189/150‑1 (MQ=255)
                                                                                                        ttGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCa                                             <  1:1046926/150‑1 (MQ=255)
                                                                                                         tGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCAATTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAg                                            <  1:141860/150‑1 (MQ=255)
                                                                                                         tGTCATAGACGGCCTGAGCCACATCAGCGCCACATAAAATGAAGACATGAATAAGAACTCGATTCATTTATTTATATCTCCGTACAGCAGGGTTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAg                                            <  1:23896/150‑1 (MQ=255)
                                                                                                          gTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAgc                                           >  2:108789/1‑149 (MQ=255)
                                                                                                          gTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGa                                           >  2:630591/1‑150 (MQ=255)
                                                                                                            cATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTAATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGaca                                         >  2:1297801/1‑150 (MQ=255)
                                                                                                                gACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCgg                                     >  1:124344/1‑150 (MQ=255)
                                                                                                                caCGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCgg                                     <  2:464151/149‑1 (MQ=255)
                                                                                                                        gcgcCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGAtgt                             <  1:1316538/150‑1 (MQ=255)
                                                                                                                         cgcCACATCAGCGCCACATTAAATGACGCCTGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGATAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGAtgtg                            <  1:856246/150‑1 (MQ=255)
                                                                                                                          gcCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATgtgt                           >  1:1004238/1‑150 (MQ=255)
                                                                                                                           ccACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTa                          >  1:978447/1‑150 (MQ=255)
                                                                                                                           ccACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCGAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTa                          <  2:1069931/150‑1 (MQ=255)
                                                                                                                            cacaTCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTaa                         >  2:1217891/1‑150 (MQ=255)
                                                                                                                            cacaTCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTaa                         >  2:2171/1‑150 (MQ=255)
                                                                                                                             acaTCAGCGCCACATAAAATGCCGCCAGGAATACGAACTCGCTTCATTTATTTATATCTCCGTCAAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGCCACCGGGTTGATGAGTAAc                        >  2:460068/1‑150 (MQ=255)
                                                                                                                                 cAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATa                    >  1:907240/1‑150 (MQ=255)
                                                                                                                                 aaGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATa                    <  1:84275/149‑1 (MQ=255)
                                                                                                                                     gcCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGa                <  2:1252529/150‑1 (MQ=255)
                                                                                                                                       cacaTCAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGAcg              <  1:95612/150‑1 (MQ=255)
                                                                                                                                       cacaTAAAATGAAGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGAcg              <  1:721293/150‑1 (MQ=255)
                                                                                                                                       cacaTAAAATGAAGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGAcg              <  1:726052/150‑1 (MQ=255)
                                                                                                                                              aaTGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTtg       >  2:561967/1‑150 (MQ=255)
                                                                                                                                              aaTGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTtg       <  1:1045967/150‑1 (MQ=255)
                                                                                                                                              aaTGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTtg       >  2:987387/1‑150 (MQ=255)
                                                                                                                                                  aCGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCCTCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTggg   <  1:1058228/150‑1 (MQ=255)
                                                                                                                                                   cGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGt  >  2:1245044/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAAACGCTCGATGCTTTCACTGGTGTAATGCGTGTTGCCGTTTACGGGTACAGGGACAGGAGCATAGACCGTCATCTCTTCAAAAAAATGCATGTACATGTTTTTGTCATCGACGGCCTGCGCCACATCAGCGCCACATAAAATGACGCCAGGAATAAGAACTCGCTTCATTTATTTATATCTCCGTACAGCAGGGCTTTTATTGTTTTTATAAAACCGCGCTAATGTCACAAATATCATCGGCCAGGTTCCAGACACCGGGTTGATGTGTAACCATAATGACGCCGCTTGTGGGT  >  NC_000913/1577448‑1577744

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