Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 251,476 A→G E193E (GAA→GAG dinP → damage‑inducible protein P

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913251,4760AG88.4% 671.4 / 68.9 241E193E (GAA→GAGdinPdamage‑inducible protein P
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base G (115/98);  total (131/110)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 245 total aligned reads displayed.

GATTACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGG  >  NC_000913/251386‑251566
                                                                                          |                                                                                          
  ttACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGa                                                                                        <  1:4806194/93‑1 (MQ=255)
  ttACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGa                                                                                        >  2:3283389/1‑93 (MQ=255)
     cGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATgg                                                                                     <  1:6643716/93‑1 (MQ=255)
      gCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATggg                                                                                    <  1:5752257/93‑1 (MQ=255)
       ccGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATgggg                                                                                   >  1:1004169/1‑93 (MQ=255)
       ccGGCAGAAGTGCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATgggg                                                                                   <  2:6441148/93‑1 (MQ=255)
        cGGCAGAGGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATTCCCTGCGTCGGCAACGTCTCAGCGGCAAAACTGGAGGCGATggggg                                                                                  >  1:5371235/1‑92 (MQ=255)
         ggCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGcg                                                                                 >  2:5585480/1‑92 (MQ=255)
         ggCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCt                                                                                 >  2:232447/1‑93 (MQ=255)
          gCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTg                                                                                <  1:2176492/93‑1 (MQ=255)
          gCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTg                                                                                >  1:4436061/1‑93 (MQ=255)
           cAGAAGTTCCGGCATTTTTACAAACCTTACTGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGc                                                                               <  1:535179/93‑1 (MQ=255)
           cAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGc                                                                               <  2:1143143/93‑1 (MQ=255)
            aGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCg                                                                              >  2:4051460/1‑93 (MQ=255)
             gAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCgg                                                                             >  2:4304112/1‑93 (MQ=255)
             gAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCgg                                                                             <  1:601999/93‑1 (MQ=255)
              aaGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGa                                                                            >  1:2055379/1‑93 (MQ=255)
               aGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGAc                                                                           >  2:126026/1‑93 (MQ=255)
               aGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGAc                                                                           >  1:3398031/1‑93 (MQ=255)
                gTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGAcc                                                                          <  2:2166458/93‑1 (MQ=255)
                gTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGAcc                                                                          <  2:3391963/93‑1 (MQ=255)
                 ttCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCt                                                                         <  1:2994274/93‑1 (MQ=255)
                  tCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTg                                                                        >  2:856072/1‑93 (MQ=255)
                     ggCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCgg                                                                     >  2:4280569/1‑93 (MQ=255)
                        aTTTTTACAAACGTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  <  1:2425979/93‑1 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGTCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  <  2:1422690/93‑1 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  <  2:5467780/93‑1 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  >  2:4000484/1‑93 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  >  2:2045587/1‑93 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGa                                                                  >  2:2657618/1‑93 (MQ=255)
                        aTTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGa                                                                  <  2:5847493/93‑1 (MQ=255)
                         tttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGAt                                                                 <  2:1671536/93‑1 (MQ=255)
                         tttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGAt                                                                 >  1:826463/1‑93 (MQ=255)
                          ttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATg                                                                <  1:1216595/93‑1 (MQ=255)
                          ttttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATg                                                                <  2:3122148/93‑1 (MQ=255)
                           tttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGt                                                               <  1:439186/93‑1 (MQ=255)
                           tttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGt                                                               <  2:2273550/93‑1 (MQ=255)
                              aCAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACa                                                            >  1:6322579/1‑93 (MQ=255)
                               cAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACaa                                                           >  1:2756573/1‑93 (MQ=255)
                               cAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACaa                                                           <  2:3693106/93‑1 (MQ=255)
                                  aCCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAg                                                        <  2:2798386/93‑1 (MQ=255)
                                   ccTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAgt                                                       <  2:3513411/93‑1 (MQ=255)
                                     ttACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAgtgt                                                     <  1:4642040/93‑1 (MQ=255)
                                     ttACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAgtgt                                                     <  2:2497532/93‑1 (MQ=255)
                                     ttACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAgtgt                                                     >  1:262608/1‑93 (MQ=255)
                                      tACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGtgtg                                                    >  1:6615275/1‑93 (MQ=255)
                                       aCCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGa                                                   <  2:1958075/93‑1 (MQ=255)
                                        ccGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGAt                                                  <  2:2288473/93‑1 (MQ=255)
                                        ccGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGAt                                                  >  2:393313/1‑93 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 <  1:1766192/93‑1 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 <  2:5114090/93‑1 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 <  1:1034747/93‑1 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 <  1:2494289/93‑1 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 >  1:3876785/1‑93 (MQ=255)
                                         cGCTGGCAAAAATCCCCGGCGTCGGCAAAGCCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATc                                                 <  1:1417015/93‑1 (MQ=255)
                                           cTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTg                                               >  1:349431/1‑93 (MQ=255)
                                              gCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTg                                            <  1:3280331/93‑1 (MQ=255)
                                               cAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGa                                           <  2:2137868/93‑1 (MQ=255)
                                                aaaaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAt                                          <  2:2218993/93‑1 (MQ=255)
                                                aaaaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAt                                          <  2:6004709/93‑1 (MQ=255)
                                                 aaaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtg                                         <  1:2043981/93‑1 (MQ=255)
                                                 aaaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtg                                         >  2:5160640/1‑93 (MQ=255)
                                                  acaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAGAGTGTGATCTGGTGAtgc                                        <  1:910015/91‑1 (MQ=255)
                                                  aaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATTGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtgc                                        >  2:4814833/1‑93 (MQ=255)
                                                  aaaTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtgc                                        >  1:3124809/1‑93 (MQ=255)
                                                    aTCCCCGGCGTCGGCAGAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtgctg                                      <  1:184553/93‑1 (MQ=255)
                                                    aTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtgctg                                      >  1:4034222/1‑93 (MQ=255)
                                                    aTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGAtgctg                                      >  2:3532821/1‑93 (MQ=255)
                                                    aTCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGCGATGTACAAAAGTGTGATCTGGTGAtgctg                                      >  2:3634604/1‑93 (MQ=255)
                                                      ccccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATgctgct                                    >  2:6337207/1‑93 (MQ=255)
                                                       cccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCtt                                   <  1:251430/93‑1 (MQ=255)
                                                       cccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCtt                                   >  2:5794390/1‑93 (MQ=255)
                                                        ccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACCAAAGTGTGATCTGGTGATGCTGCTTa                                  >  1:2572334/1‑93 (MQ=255)
                                                        ccGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTa                                  >  2:6457469/1‑93 (MQ=255)
                                                          ggCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTTCAAAAGTGTGATCTGGTGATGCTGCTTaaa                                >  2:6124313/1‑93 (MQ=255)
                                                          ggCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTaaa                                <  1:2035972/93‑1 (MQ=255)
                                                           gCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAAc                               >  1:4353251/1‑93 (MQ=255)
                                                           gCGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAAc                               >  2:1520176/1‑93 (MQ=255)
                                                            cGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACg                              <  2:1027408/93‑1 (MQ=255)
                                                            cGTCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACg                              >  1:1358511/1‑93 (MQ=255)
                                                              tCGGCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCt                            <  2:1779962/93‑1 (MQ=255)
                                                                ggCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCttt                          >  1:4175519/1‑93 (MQ=255)
                                                                 gCAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGGGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTg                         >  2:4904028/1‑93 (MQ=255)
                                                                  cAAAGTCTCAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTgg                        >  1:3670680/1‑93 (MQ=255)
                                                                       tctcAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAAt                   >  1:5214098/1‑93 (MQ=255)
                                                                       tctcAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAAt                   >  1:1688979/1‑93 (MQ=255)
                                                                         tcAGCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAAttt                 <  2:4935318/93‑1 (MQ=255)
                                                                            gCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGc              <  1:182363/93‑1 (MQ=255)
                                                                            gCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGc              <  1:5696541/93‑1 (MQ=255)
                                                                            gCGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGc              >  1:5850715/1‑93 (MQ=255)
                                                                             cGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGcc             <  1:2889109/93‑1 (MQ=255)
                                                                             cGGCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGGGATGCTGCTTAAACGCTTTGGCAAATTTGGcc             >  2:6292312/1‑93 (MQ=255)
                                                                              ggCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCg            >  1:4610784/1‑93 (MQ=255)
                                                                              ggCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCg            >  2:2893109/1‑93 (MQ=255)
                                                                              gccAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATTCTGCTTAATCGCTTTGGCAAATTTGGCCg            <  1:563282/91‑1 (MQ=255)
                                                                               gCAAAACTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGc           >  2:5405969/1‑93 (MQ=255)
                                                                                 aaaaCTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCAt         <  2:724515/93‑1 (MQ=255)
                                                                                  aaaCTGGAGGCGATAGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCAtt        >  1:6100157/1‑93 (MQ=255)
                                                                                   aaCTGGAGGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCAttt       >  2:4839997/1‑93 (MQ=255)
                                                                                        gAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTgg  >  2:4693616/1‑93 (MQ=255)
                                                                                          |                                                                                          
GATTACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGG  >  NC_000913/251386‑251566

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