Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,391,480 A→G E84G (GAG→GGG)  b1329 → putative transport periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,391,4800AG87.3% 587.8 / 74.2 237E84G (GAG→GGG) b1329putative transport periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (9/21);  new base G (89/118);  total (98/139)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

Only 100 of 254 total aligned reads displayed.

TCGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGATG  >  NC_000913/1391389‑1391569
                                                                                           |                                                                                         
 cGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGggga                                                                                         <  2:6711200/93‑1 (MQ=255)
 cGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGggga                                                                                         <  1:7195365/93‑1 (MQ=255)
 cGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGgaga                                                                                         <  1:1649130/93‑1 (MQ=255)
 cGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCCATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGggga                                                                                         >  1:7376877/1‑93 (MQ=255)
     ggATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGt                                                                                     >  2:3102584/1‑93 (MQ=255)
        tCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTccc                                                                                  <  1:2650079/93‑1 (MQ=255)
        tCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTccc                                                                                  <  2:6601476/93‑1 (MQ=255)
        tCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTccc                                                                                  <  1:255672/93‑1 (MQ=255)
         cccGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTcccc                                                                                 <  1:4884644/93‑1 (MQ=255)
           cGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCgg                                                                               <  2:4989684/93‑1 (MQ=255)
            gCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGc                                                                              >  1:5654764/1‑93 (MQ=255)
             cTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCg                                                                             >  1:1775535/1‑93 (MQ=255)
              tAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGt                                                                            <  2:578768/93‑1 (MQ=255)
              tAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGt                                                                            >  1:4310399/1‑93 (MQ=255)
                 aGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGc                                                                         >  2:4716915/1‑93 (MQ=255)
                 aGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGc                                                                         <  1:607016/93‑1 (MQ=255)
                  gCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCg                                                                        >  1:3914587/1‑93 (MQ=255)
                  gCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCg                                                                        >  1:1953953/1‑93 (MQ=255)
                    cGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGAc                                                                      <  2:5639998/93‑1 (MQ=255)
                    cGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGAc                                                                      >  1:775326/1‑93 (MQ=255)
                       gggCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCa                                                                   <  2:1217059/93‑1 (MQ=255)
                         gCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGt                                                                 >  2:1692066/1‑93 (MQ=255)
                          ccTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTg                                                                <  1:3031464/93‑1 (MQ=255)
                             gCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGaa                                                             >  2:4085659/1‑93 (MQ=255)
                               cAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGaaaa                                                           <  2:13550/93‑1 (MQ=255)
                                 gagaTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGt                                                         >  1:6721069/1‑93 (MQ=255)
                                  agaTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTa                                                        <  2:1569917/93‑1 (MQ=255)
                                   gaTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTaa                                                       <  1:4562630/93‑1 (MQ=255)
                                   gaTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTaa                                                       >  1:5738540/1‑93 (MQ=255)
                                   gaTTCACGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTaa                                                       >  2:7297176/1‑93 (MQ=255)
                                    aTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAAt                                                      <  2:5403104/93‑1 (MQ=255)
                                    aTTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAAt                                                      >  2:5239440/1‑93 (MQ=255)
                                     ttCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATg                                                     >  1:7001937/1‑93 (MQ=255)
                                      tCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGa                                                    >  1:1384369/1‑93 (MQ=255)
                                      tCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGa                                                    >  2:3590509/1‑93 (MQ=255)
                                       cAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGAc                                                   <  1:1194070/93‑1 (MQ=255)
                                       cAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGAc                                                   <  2:4962708/93‑1 (MQ=255)
                                        aGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACa                                                  >  1:2399821/1‑93 (MQ=255)
                                        aGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACa                                                  >  2:6207610/1‑93 (MQ=255)
                                        aGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACa                                                  <  1:5866530/93‑1 (MQ=255)
                                         ggTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACaa                                                 >  2:4636636/1‑93 (MQ=255)
                                         ggTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACaa                                                 <  1:1594739/93‑1 (MQ=255)
                                         ggTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACaa                                                 <  1:620027/93‑1 (MQ=255)
                                          gTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAAc                                                <  2:4298699/93‑1 (MQ=255)
                                          gTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAAc                                                >  1:1425931/1‑93 (MQ=255)
                                            cATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCg                                              <  1:4077794/93‑1 (MQ=255)
                                              ttCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTa                                            >  2:893711/1‑93 (MQ=255)
                                              ttCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTa                                            >  1:4554978/1‑93 (MQ=255)
                                              ttCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTa                                            >  1:3613044/1‑93 (MQ=255)
                                               tCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTAt                                           <  1:5071218/93‑1 (MQ=255)
                                               tCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTAt                                           <  2:7523146/93‑1 (MQ=255)
                                                cgcgATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATc                                          <  1:6349557/93‑1 (MQ=255)
                                                cgcgATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATc                                          <  2:4471408/93‑1 (MQ=255)
                                                 gcgATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCt                                         >  2:2060642/1‑93 (MQ=255)
                                                  cgATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTg                                        <  2:2495982/93‑1 (MQ=255)
                                                  cgATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTg                                        <  1:5415501/93‑1 (MQ=255)
                                                   gATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTgg                                       >  1:6752836/1‑93 (MQ=255)
                                                   gATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTgg                                       >  2:6785753/1‑93 (MQ=255)
                                                   gATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTgg                                       >  2:4224938/1‑93 (MQ=255)
                                                   gATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTgg                                       >  2:2716722/1‑93 (MQ=255)
                                                   gATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTgg                                       >  1:2373262/1‑93 (MQ=255)
                                                    aTCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGa                                      <  2:2544948/93‑1 (MQ=255)
                                                    aTCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGa                                      <  2:2969233/93‑1 (MQ=255)
                                                     tCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGAc                                     >  2:7047818/1‑93 (MQ=255)
                                                     tCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGAc                                     <  1:5622621/93‑1 (MQ=255)
                                                       tGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACtt                                   >  2:224014/1‑93 (MQ=255)
                                                       tGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACtt                                   >  1:2850078/1‑93 (MQ=255)
                                                           tGAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAc                               <  1:5127829/93‑1 (MQ=255)
                                                            gAAGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAcc                              >  1:5358857/1‑93 (MQ=255)
                                                             aaGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAccc                             >  2:1288682/1‑93 (MQ=255)
                                                             aaGGTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTAccc                             >  2:4825886/1‑93 (MQ=255)
                                                               ggTCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAACGACAACCGTATCTGGACTTTTACCCTg                           <  2:6404346/93‑1 (MQ=255)
                                                                 tCTGGTGAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTgcg                         <  1:2054161/93‑1 (MQ=255)
                                                                       gAATCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATaa                   <  2:3896803/93‑1 (MQ=255)
                                                                        aaTCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAAc                  <  2:274900/93‑1 (MQ=255)
                                                                         aTCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACg                 >  1:3531653/1‑93 (MQ=255)
                                                                         aTCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACg                 >  2:5086955/1‑93 (MQ=255)
                                                                          tCAGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGc                >  1:5498448/1‑93 (MQ=255)
                                                                            aGAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaa              <  2:7432591/93‑1 (MQ=255)
                                                                             gAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaaa             >  1:2722010/1‑93 (MQ=255)
                                                                             gAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaaa             >  1:4389758/1‑93 (MQ=255)
                                                                             gAACGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaaa             <  2:1035481/93‑1 (MQ=255)
                                                                              aaCGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaaaa            <  2:3440843/93‑1 (MQ=255)
                                                                              aaCGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCaaaa            <  1:4114599/93‑1 (MQ=255)
                                                                               aCGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAAt           <  1:6534641/93‑1 (MQ=255)
                                                                                cGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATg          <  1:3333491/93‑1 (MQ=255)
                                                                                cGAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATg          <  2:2633447/93‑1 (MQ=255)
                                                                                 gAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATgg         >  2:461971/1‑93 (MQ=255)
                                                                                 gAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATgg         <  1:4673005/93‑1 (MQ=255)
                                                                                 gAAAAAGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATgg         <  2:7169123/93‑1 (MQ=255)
                                                                                  aaaaaGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATggg        <  1:4150421/93‑1 (MQ=255)
                                                                                  aaaaaGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATggg        <  2:6928207/93‑1 (MQ=255)
                                                                                   aaaaGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGc       <  1:3866078/93‑1 (MQ=255)
                                                                                     aaGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCgg     >  1:2088776/1‑93 (MQ=255)
                                                                                      aGGGGGGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGa    <  2:1448533/93‑1 (MQ=255)
                                                                                       ggggggATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGAt   <  1:1082459/93‑1 (MQ=255)
                                                                                       ggggggATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGAt   <  2:3360176/93‑1 (MQ=255)
                                                                                       ggggggATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGAt   >  2:4159817/1‑93 (MQ=255)
                                                                                       ggggAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGAt   <  2:5996719/93‑1 (MQ=255)
                                                                                        ggggGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGATg  <  1:3218592/93‑1 (MQ=255)
                                                                                           |                                                                                         
TCGCTGGATCCCGCTAAAGCCGTGGGCCTGCCAGAGATTCAGGTCATTCGCGATCTGTTTGAAGGTCTGGTGAATCAGAACGAAAAAGGGGAGATTGTCCCCGGCGTTGCGACTCAGTGGAAAAGTAATGACAACCGTATCTGGACTTTTACCCTGCGCGATAACGCAAAATGGGCGGATG  >  NC_000913/1391389‑1391569

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