Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,806,223 A→G I19M (ATA→ATG kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,2230AG98.7% 1703.5 / 3.3 536I19M (ATA→ATGkdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base G (266/263);  total (272/264)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 549 total aligned reads displayed.

AATCCCCAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGG  >  NC_000913/3806133‑3806312
                                                                                          |                                                                                         
aaTCCCCAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATTCg                                                                                         >  1:6243029/1‑93 (MQ=255)
aaTCCCCAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCg                                                                                         <  2:6604149/93‑1 (MQ=255)
    cccAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGctc                                                                                     <  1:487891/93‑1 (MQ=255)
    cccAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGctc                                                                                     <  1:1867188/93‑1 (MQ=255)
     ccAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGctct                                                                                    >  2:5399000/1‑93 (MQ=255)
     ccAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGctct                                                                                    >  1:7814396/1‑93 (MQ=255)
     ccAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGctct                                                                                    >  1:2969115/1‑93 (MQ=255)
       aGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTgg                                                                                  >  2:80300/1‑93 (MQ=255)
        gCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTggg                                                                                 >  2:2631071/1‑93 (MQ=255)
        gCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTggg                                                                                 >  1:5966077/1‑93 (MQ=255)
        gCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTggg                                                                                 >  1:8030638/1‑93 (MQ=255)
        gCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTggg                                                                                 >  2:3166554/1‑93 (MQ=255)
         cacaTCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGt                                                                                <  1:3827911/93‑1 (MQ=255)
           caTCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTgc                                                                              <  1:9258985/93‑1 (MQ=255)
             tCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTTGGTgcgc                                                                            >  1:3096402/1‑93 (MQ=255)
             tCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTgcgt                                                                            >  2:8196949/1‑92 (MQ=255)
             tCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTgcgc                                                                            <  2:3291566/93‑1 (MQ=255)
             tCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTgcgc                                                                            >  1:5425351/1‑93 (MQ=255)
               cATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCgg                                                                          <  1:6146596/93‑1 (MQ=255)
                 tAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGAc                                                                        >  1:7563554/1‑93 (MQ=255)
                 tAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCCCGCAc                                                                        >  1:6350726/1‑93 (MQ=255)
                   aGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGt                                                                      >  1:2677350/1‑93 (MQ=255)
                    gTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTa                                                                     <  2:4866928/93‑1 (MQ=255)
                     tCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTaa                                                                    <  1:1674605/93‑1 (MQ=255)
                     tCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTaa                                                                    <  2:8674132/93‑1 (MQ=255)
                     tCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTaa                                                                    <  1:1443192/93‑1 (MQ=255)
                     tCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTaa                                                                    >  1:8876184/1‑93 (MQ=255)
                     gcAGCTATTTACTAGCCTCGACTTGTTTTACCCCGCTCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTaa                                                                    <  2:3895094/92‑1 (MQ=255)
                       aGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAgg                                                                  >  2:6069152/1‑93 (MQ=255)
                       aGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAgg                                                                  >  1:9362647/1‑93 (MQ=255)
                       aGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTAAgg                                                                  >  1:9861944/1‑93 (MQ=255)
                        gCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGc                                                                 <  1:8083817/93‑1 (MQ=255)
                        gCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGc                                                                 <  1:6537427/93‑1 (MQ=255)
                         cTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCt                                                                >  2:2417356/1‑93 (MQ=255)
                          tATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTc                                                               >  2:8461671/1‑93 (MQ=255)
                          tATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTc                                                               >  1:2404527/1‑93 (MQ=255)
                           aTTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTcc                                                              >  1:6656141/1‑93 (MQ=255)
                           aTTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTcc                                                              <  2:3429973/93‑1 (MQ=255)
                           aTTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTcc                                                              <  2:703292/93‑1 (MQ=255)
                           aTTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTAAGGCTcc                                                              >  2:6999484/1‑93 (MQ=255)
                            tttACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCg                                                             >  1:1203503/1‑93 (MQ=255)
                              tACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGc                                                           >  2:1162106/1‑93 (MQ=255)
                               aCTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGcc                                                          >  1:4754969/1‑93 (MQ=255)
                               aCTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGcc                                                          <  1:7031375/93‑1 (MQ=255)
                                cTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCt                                                         <  1:2063124/93‑1 (MQ=255)
                                  aTGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTAt                                                       >  2:5697242/1‑93 (MQ=255)
                                  aTGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTAt                                                       >  1:9059886/1‑93 (MQ=255)
                                   tGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATc                                                      >  2:4253965/1‑93 (MQ=255)
                                   tGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATc                                                      >  1:9622719/1‑93 (MQ=255)
                                     cTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTACGGCTCCGGCCTATCGa                                                    >  2:6424000/1‑93 (MQ=255)
                                      tCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGaa                                                   >  2:7162698/1‑93 (MQ=255)
                                        gAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGaaaa                                                 <  2:492339/93‑1 (MQ=255)
                                        gAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGaaaa                                                 <  1:2550127/93‑1 (MQ=255)
                                        gAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGaaaa                                                 <  1:1024737/93‑1 (MQ=255)
                                         aaTTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAAc                                                <  1:6889586/93‑1 (MQ=255)
                                          aTTGCTTTACACCACCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACg                                               >  2:6863013/1‑93 (MQ=255)
                                           ttGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGc                                              <  2:9509569/93‑1 (MQ=255)
                                              cTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTgg                                           <  2:1048935/93‑1 (MQ=255)
                                                ttACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTgggg                                         <  2:305370/93‑1 (MQ=255)
                                                  acacCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTg                                       >  1:3986735/1‑93 (MQ=255)
                                                   cacCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGa                                      >  2:6818967/1‑93 (MQ=255)
                                                   cacCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGa                                      <  2:7398026/93‑1 (MQ=255)
                                                      cGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACg                                   <  2:3831709/93‑1 (MQ=255)
                                                       gCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGt                                  >  2:775363/1‑93 (MQ=255)
                                                        cccTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGtt                                 <  2:4719639/93‑1 (MQ=255)
                                                         ccTTCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTa                                <  1:5056416/93‑1 (MQ=255)
                                                           ttCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACg                              >  1:5251718/1‑93 (MQ=255)
                                                           ttCTCTACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACg                              >  1:1942889/1‑93 (MQ=255)
                                                             ctctACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGt                            <  1:7462949/93‑1 (MQ=255)
                                                                tACCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGtttt                         >  2:5321040/1‑93 (MQ=255)
                                                                 aCCTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGttttt                        >  2:142342/1‑93 (MQ=255)
                                                                   cTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTAc                      <  1:1128228/93‑1 (MQ=255)
                                                                   cTTATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTAc                      <  1:1308965/93‑1 (MQ=255)
                                                                    ttATTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTAcc                     <  1:3341685/93‑1 (MQ=255)
                                                                      aTTCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGc                   >  1:838893/1‑93 (MQ=255)
                                                                        tCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACTCTGGGGTGAAAGTTACGGTTTTTACCGCCa                 >  1:2227894/1‑93 (MQ=255)
                                                                        tCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCa                 <  1:3452745/93‑1 (MQ=255)
                                                                        tCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCa                 >  1:8561806/1‑93 (MQ=255)
                                                                        tCAGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCa                 <  2:1584491/93‑1 (MQ=255)
                                                                          aGCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATc               <  1:8675206/93‑1 (MQ=255)
                                                                          aGCCGCTGATCTGGATGCGGCTCTGGGTGCGAGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATc               <  1:56009/93‑1 (MQ=255)
                                                                           gCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAACCGCTGGGGGGAACGTTACGGTTTTTACCGCCATcc              >  1:8254096/1‑93 (MQ=255)
                                                                           gCCGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATcc              >  2:586697/1‑93 (MQ=255)
                                                                           gCCGCGGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATcc              <  2:2171623/93‑1 (MQ=255)
                                                                             cGCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGc            >  1:5710623/1‑93 (MQ=255)
                                                                              gCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCt           >  2:5105590/1‑93 (MQ=255)
                                                                              gCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCt           <  2:4047692/93‑1 (MQ=255)
                                                                              gCTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCt           <  2:3548550/93‑1 (MQ=255)
                                                                               cTGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTa          >  2:1766963/1‑93 (MQ=255)
                                                                               cTGATCTGGATGCGGCTCTGAGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTa          >  1:4595084/1‑93 (MQ=255)
                                                                                tGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTaa         <  1:447051/93‑1 (MQ=255)
                                                                                tGATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTaa         <  2:9213774/93‑1 (MQ=255)
                                                                                 gATCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTaaa        <  1:4274419/93‑1 (MQ=255)
                                                                                  aTCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTaaaa       <  1:7135983/93‑1 (MQ=255)
                                                                                   tCTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAAc      >  1:8820755/1‑93 (MQ=255)
                                                                                    cTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAAcc     <  2:4484172/93‑1 (MQ=255)
                                                                                    cTGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAAcc     <  2:228035/93‑1 (MQ=255)
                                                                                     tGGATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCa    <  2:2966245/93‑1 (MQ=255)
                                                                                       gATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAgg  <  2:7244169/93‑1 (MQ=255)
                                                                                       gATGCGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAgg  <  2:2741078/93‑1 (MQ=255)
                                                                                          |                                                                                         
AATCCCCAGCACATCCATAAGTCAGCTATTTACTATGCTCGAATTGCTTTACACCGCCCTTCTCTACCTTATTCAGCCGCTGATCTGGATACGGCTCTGGGTGCGCGGACGTAAGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGG  >  NC_000913/3806133‑3806312

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