Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,117,945 G→A P451S (CCG→TCG)  yghK ← putative permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,117,9450GA86.2% 562.1 / 59.5 195P451S (CCG→TCG) yghKputative permease
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base A (115/53);  total (127/68)
Fisher's exact test for biased strand distribution p-value = 2.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 222 total aligned reads displayed.

GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC  >  NC_000913/3117853‑3118034
                                                                                            |                                                                                         
gtgcGATTGCAGTGAACCAAACGGGGCATAAGGGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACgg                                                                                           <  1:2384772/90‑1 (MQ=25)
ggTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACgg                                                                                           >  2:2473156/1‑93 (MQ=255)
   ggtttGCAGTAACCAAAACGGGGCGTAAGAGGAGGTGCCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGaa                                                                                        <  2:5777063/90‑1 (MQ=21)
   cGATTGCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGTTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGaa                                                                                        <  1:5343256/93‑1 (MQ=25)
    gATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAAc                                                                                       >  1:968672/1‑93 (MQ=255)
     aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACa                                                                                      >  2:2969266/1‑93 (MQ=255)
     aTTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACa                                                                                      >  1:1902201/1‑93 (MQ=255)
       tGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATc                                                                                    >  1:1900505/1‑93 (MQ=255)
        gCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                                                   >  1:603504/1‑93 (MQ=255)
        gCAGTGAACAAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCa                                                                                   <  1:203284/93‑1 (MQ=255)
         cAGTGACCAAAACGGGCCGTTAGAGGAGGTGCCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCAc                                                                                  <  1:6403471/93‑1 (MQ=25)
            tGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGcc                                                                               >  1:4506266/1‑93 (MQ=255)
            gaacccAAACGGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGcc                                                                               <  2:3814939/89‑1 (MQ=255)
              aaCCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTg                                                                             >  2:3443462/1‑93 (MQ=255)
                ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTa                                                                           >  2:124541/1‑93 (MQ=255)
                ccAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTa                                                                           >  2:4132774/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  2:6301940/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  1:4662038/1‑93 (MQ=255)
                 cAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAc                                                                          >  2:4318449/1‑93 (MQ=255)
                  aaaaAGGGTGTTAGGGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTAcc                                                                         <  1:1958748/93‑1 (MQ=25)
                   aaCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCt                                                                        >  2:6010800/1‑93 (MQ=255)
                    acgggggGTAAGAGGAGGGGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTg                                                                       <  2:2557488/90‑1 (MQ=255)
                    aCAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTg                                                                       >  1:4996463/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCa                                                                    >  2:889654/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCa                                                                    >  1:5443636/1‑93 (MQ=255)
                       gggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACAGCACGCCTGTACCTGCCa                                                                    >  1:3720094/1‑93 (MQ=255)
                        ggCGTTAGCGGAGGTGTCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAg                                                                   <  2:854557/93‑1 (MQ=37)
                        ggCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAg                                                                   >  1:3105561/1‑93 (MQ=255)
                           gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                                                                >  1:109029/1‑93 (MQ=255)
                           gTTAGAGGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAc                                                                >  2:3614032/1‑93 (MQ=255)
                             tgggggaggTTCCCGAGCGGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAcca                                                              <  2:1459488/89‑1 (MQ=255)
                              agagGAGGTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccag                                                             >  2:4980090/1‑93 (MQ=255)
                               gaggagGTGCCGGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagc                                                            <  2:5394554/93‑1 (MQ=255)
                                aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATGACGCCTGTACCTGCCAGTAccagcg                                                           >  1:5513416/1‑93 (MQ=255)
                                aggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcg                                                           >  1:4464902/1‑93 (MQ=255)
                                 ggaggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgc                                                          >  1:2752082/1‑93 (MQ=255)
                                  gaggTACCCGACCCGTTAAGAAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAAAACATCACGCCTGTACCTGCCAGTAccagcgcc                                                         <  1:3803865/93‑1 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  2:5483829/1‑93 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  1:3610822/1‑93 (MQ=255)
                                    ggTGTCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccag                                                       >  1:4335208/1‑93 (MQ=255)
                                     gtgtCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTAccagcgccagc                                                      >  1:2960456/1‑93 (MQ=255)
                                        tCCGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtg                                                   >  2:1511734/1‑93 (MQ=255)
                                        gcggAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtg                                                   <  2:4475702/90‑1 (MQ=255)
                                          cGAGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt                                                 <  1:6149973/93‑1 (MQ=255)
                                          cGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCgtggt                                                 <  1:4099929/93‑1 (MQ=255)
                                            aGCCGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGtggtgg                                               >  1:659453/1‑93 (MQ=255)
                                              cgggTAAGGAATGCGCACGGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAc                                             <  1:743423/91‑1 (MQ=255)
                                              ccGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAc                                             >  2:4654746/1‑93 (MQ=255)
                                               cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGACCGAGAACATCACGCTTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa                                            <  2:5914295/93‑1 (MQ=255)
                                               cGGTAAGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCGGCGCCAGCGTGGTGGACa                                            >  1:5908432/1‑93 (MQ=255)
                                               cGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACa                                            <  2:6466579/93‑1 (MQ=255)
                                                 gTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAcgtc                                          >  2:512525/1‑90 (MQ=255)
                                                    aGGAATGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCa                                       <  2:2680041/93‑1 (MQ=255)
                                                      gAATCCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGa                                     <  1:3205736/93‑1 (MQ=255)
                                                      gAATCCCCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGa                                     <  1:1684850/93‑1 (MQ=38)
                                                       aaggcGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGaa                                    <  2:1130456/89‑1 (MQ=255)
                                                       aaTGCGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGaa                                    >  2:2883552/1‑93 (MQ=255)
                                                         gccGCCGAGCAAGCCGAGAATCGCTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGTACATGCCAGAATa                                  <  2:2171877/91‑1 (MQ=255)
                                                          gcGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAg                                 <  2:3204526/92‑1 (MQ=255)
                                                             cccAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTg                              >  1:4215726/1‑93 (MQ=255)
                                                                agccagccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTg                           >  2:4178306/1‑93 (MQ=255)
                                                                 gccagccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGa                          >  1:5392979/1‑93 (MQ=255)
                                                                   cagccGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACg                        >  2:369176/1‑93 (MQ=255)
                                                                    agccGAGAAACGGTGAGAAGAACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGa                       <  1:1606470/93‑1 (MQ=255)
                                                                    agccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGa                       >  2:6169871/1‑93 (MQ=255)
                                                                     gcggagaAAGGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:1195833/90‑1 (MQ=255)
                                                                     gccGAGAAAGGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACgca                      <  1:3465831/93‑3 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      >  1:2085463/1‑93 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:1860469/93‑1 (MQ=255)
                                                                     gccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGaa                      <  2:696140/93‑1 (MQ=255)
                                                                      ccGAGAAAGTGTGAGAAGAACGAGAACCTCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  2:5592849/93‑1 (MQ=255)
                                                                      ccGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  1:6367213/93‑1 (MQ=255)
                                                                      ccGAGAAACGCTGAGAAGAGCGAGACCTTCCCGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAg                     <  1:5455319/93‑1 (MQ=255)
                                                                       cGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTCCCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAgg                    <  2:5477515/93‑1 (MQ=255)
                                                                       cGAGAAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGGGACGAAgg                    >  2:656329/1‑93 (MQ=255)
                                                                        gagaAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGc                   >  1:413429/1‑93 (MQ=255)
                                                                        gagaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGc                   <  2:4915628/93‑1 (MQ=255)
                                                                         agaAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCg                  <  1:1141775/93‑1 (MQ=255)
                                                                         agaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCg                  >  1:4001152/1‑93 (MQ=255)
                                                                          gaAACGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGa                 >  2:4550675/1‑93 (MQ=255)
                                                                           aaaCGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGaa                <  1:4292591/93‑1 (MQ=255)
                                                                           aaaCGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGaa                >  1:5841928/1‑93 (MQ=255)
                                                                            aaCGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAAc               >  2:441329/1‑93 (MQ=255)
                                                                              cGGTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGc             <  2:5193632/93‑1 (MQ=255)
                                                                               ggTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGcc            >  2:169296/1‑93 (MQ=255)
                                                                               ggTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGcc            <  1:6357012/93‑1 (MQ=255)
                                                                                gTGAGAAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGCACGCCa           >  2:1400391/1‑93 (MQ=255)
                                                                                 ggagaAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAg          <  2:2739730/92‑1 (MQ=255)
                                                                                  gagaAGAACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGc         >  2:4378743/1‑93 (MQ=255)
                                                                                   agaagaACGAGACCATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCa        <  1:1158615/93‑1 (MQ=255)
                                                                                   agaagaACGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCa        >  1:4969150/1‑93 (MQ=255)
                                                                                    gaagaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAc       <  2:5373781/93‑1 (MQ=255)
                                                                                    aaagaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAc       <  2:2629817/92‑1 (MQ=255)
                                                                                     aagaaCGGGAACATCCCGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAcc      <  2:4274395/93‑1 (MQ=255)
                                                                                     aagaaCGGAAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCAcc      <  2:6397433/93‑1 (MQ=255)
                                                                                      agaaCGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCa     >  2:5663711/1‑93 (MQ=255)
                                                                                      agaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCa     >  2:4844919/1‑93 (MQ=255)
                                                                                      agaaCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGAAATGCCAGAATAGTTGGTGACGAAGGCGCACGCCCGCAccc     >  1:5158560/1‑92 (MQ=255)
                                                                                         aCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGc  <  2:330031/93‑1 (MQ=255)
                                                                                         aCGAGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGc  >  1:250222/1‑93 (MQ=255)
                                                                                            |                                                                                         
GGTCGATTGCAGTGAACCAAACAGGGCGTTAGAGGAGGTGTCCGAGCCGGTAAGGAATACGCCCAGCCAGCCGAGAAACGGTGAGAAGAACGGGAACATCACGCCTGTACCTGCCAGTACCAGCGCCAGCGTGGTGGACATGCCAGAATAGTTGGTGACGAAGGCGAACGCCAGCACCATGC  >  NC_000913/3117853‑3118034

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