Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,782,160 T→C V369A (GTG→GCG)  ydiD → putative ligase/synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,782,1600TC99.0% 292.3 / ‑5.0 97V369A (GTG→GCG) ydiDputative ligase/synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (37/59);  total (38/59)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTTCGCCGCATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGAT  >  NC_000913/1782070‑1782248
                                                                                          |                                                                                        
gTTCGCCGCATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCgg                                                                                        <  1:837891/93‑1 (MQ=255)
 ttCGCCGCATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGt                                                                                       >  1:1610681/1‑93 (MQ=255)
  tCGCCGCATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTc                                                                                      <  1:1376698/93‑1 (MQ=255)
    gccgcATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGa                                                                                    <  1:1100739/93‑1 (MQ=255)
    gccgcATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGa                                                                                    <  2:4217899/93‑1 (MQ=255)
     ccgcATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGAt                                                                                   >  2:2170603/1‑93 (MQ=255)
       gcATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGa                                                                                 <  1:1518710/93‑1 (MQ=255)
       gcATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGa                                                                                 <  1:4128331/93‑1 (MQ=255)
        cATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGac                                                                                >  1:3918072/1‑93 (MQ=255)
        cATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGac                                                                                <  1:660485/93‑1 (MQ=255)
         aTGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGccg                                                                               >  1:856916/1‑93 (MQ=255)
          tGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGacgc                                                                              <  2:209222/93‑1 (MQ=255)
          tGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGacgc                                                                              >  1:4343451/1‑93 (MQ=255)
           gCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGacgca                                                                             <  1:664306/93‑1 (MQ=255)
            cGGTGGTGATTCTCGATGATCCTTTGTCGCGCTTTCTGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTTGATGacgcac                                                                            <  2:2301819/93‑1 (MQ=255)
            cGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGacgcac                                                                            <  1:4372336/93‑1 (MQ=255)
              gtggtgAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTGAAGCGGTCGATGacgccggc                                                                          >  2:856394/1‑89 (MQ=255)
                  tgAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGa                                                                      <  2:385369/93‑1 (MQ=255)
                   gAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGAc                                                                     <  1:491982/93‑1 (MQ=255)
                   gAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGAc                                                                     >  1:2972075/1‑93 (MQ=255)
                   gAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGAc                                                                     >  1:2726343/1‑93 (MQ=255)
                     aTCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCt                                                                   <  2:4009304/93‑1 (MQ=255)
                      tctcGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCtt                                                                  >  1:2699560/1‑93 (MQ=255)
                      tctcGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCtt                                                                  >  2:2316517/1‑93 (MQ=255)
                      tctcGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCtt                                                                  <  1:268506/93‑1 (MQ=255)
                       ctcGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTa                                                                 <  2:3915030/93‑1 (MQ=255)
                          gatgatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTACCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCg                                                              >  1:2514564/1‑93 (MQ=255)
                           attatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGc                                                             <  2:1536139/90‑1 (MQ=255)
                           atgatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGc                                                             <  2:4032065/93‑1 (MQ=255)
                             gatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCa                                                           >  1:93927/1‑93 (MQ=255)
                              atCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAg                                                          <  2:2814222/93‑1 (MQ=255)
                                ccTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGt                                                        >  2:1736146/1‑93 (MQ=255)
                                 cTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGtt                                                       <  1:220005/93‑1 (MQ=255)
                                  tttGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTg                                                      >  1:1882157/1‑93 (MQ=255)
                                      tCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGaa                                                  >  2:1590647/1‑93 (MQ=255)
                                      tCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGaa                                                  <  2:2843175/93‑1 (MQ=255)
                                       cgcgCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAg                                                 <  2:325232/93‑1 (MQ=255)
                                        gcgcTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAgg                                                >  2:2888652/1‑93 (MQ=255)
                                         cgcTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGCGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGt                                               <  2:3069758/93‑1 (MQ=255)
                                           cTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTga                                             <  1:2619285/93‑1 (MQ=255)
                                             ttATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTgaag                                           <  1:1393333/93‑1 (MQ=255)
                                              tATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTgaaga                                          <  2:430305/93‑1 (MQ=255)
                                               aTGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTgaagaa                                         <  1:2010677/93‑1 (MQ=255)
                                               aTGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTgaagaa                                         <  1:1648131/93‑1 (MQ=255)
                                                tGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGaagaag                                        <  1:2808857/93‑1 (MQ=255)
                                                   acacCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCt                                     <  1:2116314/93‑1 (MQ=255)
                                                   acacCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCt                                     <  2:836997/93‑1 (MQ=255)
                                                    cacCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTc                                    <  2:2428092/93‑1 (MQ=255)
                                                     acCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcg                                   >  1:3719517/1‑93 (MQ=255)
                                                      ccGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgc                                  <  2:972391/93‑1 (MQ=255)
                                                       cGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgcg                                 >  1:2887083/1‑93 (MQ=255)
                                                       cGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgcg                                 <  1:4502523/93‑1 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgcgg                                >  1:3077964/1‑92 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgcgg                                >  1:2163404/1‑92 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGt                                >  1:4293151/1‑93 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGt                                >  1:4145408/1‑93 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGt                                >  2:1814245/1‑93 (MQ=255)
                                                        gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGt                                >  1:783764/1‑93 (MQ=255)
                                                         aTGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTg                               >  1:473111/1‑93 (MQ=255)
                                                         aTGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTg                               >  1:3046148/1‑93 (MQ=255)
                                                         aTGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTg                               >  2:1994255/1‑93 (MQ=255)
                                                         aTGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTg                               <  2:4420458/93‑1 (MQ=255)
                                                          tGGTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTgg                              >  1:2965128/1‑93 (MQ=255)
                                                           ggTTACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAATACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGc                             <  1:3300082/93‑1 (MQ=255)
                                                             ttACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGccc                           <  1:2885784/93‑1 (MQ=255)
                                                             ttACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGccc                           <  1:1774785/93‑1 (MQ=255)
                                                              tACGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGcccc                          <  2:3956073/93‑1 (MQ=255)
                                                               aCGCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCa                         <  2:3287687/93‑1 (MQ=255)
                                                                 gCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                       >  2:4271375/1‑93 (MQ=255)
                                                                 gCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                       <  2:512189/93‑1 (MQ=255)
                                                                 gCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                       <  2:73155/93‑1 (MQ=255)
                                                                 gCTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                       <  1:2381370/93‑1 (MQ=255)
                                                                  cTGCGGCAGGTTTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtg                      <  2:2047149/93‑1 (MQ=255)
                                                                  cTGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtg                      >  1:2744101/1‑93 (MQ=255)
                                                                   tGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgt                     <  1:2963182/93‑1 (MQ=255)
                                                                   tGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgt                     >  1:4258151/1‑93 (MQ=255)
                                                                   tGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgt                     <  2:3635379/93‑1 (MQ=255)
                                                                   tGCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgt                     <  2:557088/93‑1 (MQ=255)
                                                                    gCCGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgtg                    <  2:3992379/93‑1 (MQ=255)
                                                                    gCCGCAGGTGCAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgtg                    >  2:3395829/1‑93 (MQ=255)
                                                                     ccGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATgtgt                   <  2:4083261/93‑1 (MQ=255)
                                                                      cGCAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGtt                  <  1:3341967/93‑1 (MQ=255)
                                                                        catgtgtAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                <  2:525730/90‑1 (MQ=255)
                                                                        cAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                <  1:2664029/93‑1 (MQ=255)
                                                                        cAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                <  2:1458512/93‑1 (MQ=255)
                                                                        cAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                >  2:2944791/1‑93 (MQ=255)
                                                                        cAGGTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                <  2:4530625/93‑1 (MQ=255)
                                                                          ggTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATg              <  2:2676483/93‑1 (MQ=255)
                                                                          ggTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATg              <  1:1997783/93‑1 (MQ=255)
                                                                          ggTGTAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGACTCGCGTGGCCCCAATGTGTTTATg              >  1:1637687/1‑93 (MQ=255)
                                                                            tgtAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATggg            >  2:2846882/1‑93 (MQ=255)
                                                                              tAGAGATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGt          >  2:343082/1‑93 (MQ=255)
                                                                               agagATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTa         <  1:2015021/93‑1 (MQ=255)
                                                                               agagATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTa         <  1:1700833/93‑1 (MQ=255)
                                                                               agagATTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGACTCGCGTGGCCCCAATGTGTTTATggggga         >  1:3050015/1‑91 (MQ=255)
                                                                                 agaTTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtt       >  2:2249688/1‑93 (MQ=255)
                                                                                   aTTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtttt     <  2:4654353/93‑1 (MQ=255)
                                                                                   aTTAAAGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtttt     >  1:214838/1‑93 (MQ=255)
                                                                                      aaaGCGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgat  <  2:3536757/93‑1 (MQ=255)
                                                                                          |                                                                                        
GTTCGCCGCATGCGGTGGTGAATCTCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGAT  >  NC_000913/1782070‑1782248

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