Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 891,504 A→G Q105Q (CAA→CAG rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913891,5040AG98.7% 480.5 / ‑3.0 152Q105Q (CAA→CAGrimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (83/67);  total (84/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 157 total aligned reads displayed.

GGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGG  >  NC_000913/891414‑891593
                                                                                          |                                                                                         
ggcACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGct                                                                                         >  1:3188139/1‑93 (MQ=255)
ggcACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGct                                                                                         <  1:4706456/93‑1 (MQ=255)
 gcACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGctg                                                                                        >  2:331530/1‑93 (MQ=255)
  cACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGctgc                                                                                       <  2:4509549/93‑1 (MQ=255)
    cTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGctgctg                                                                                     >  2:1414568/1‑93 (MQ=255)
    cTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGctgctg                                                                                     >  1:4453733/1‑93 (MQ=255)
    cTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGctgctg                                                                                     >  1:1648386/1‑93 (MQ=255)
     tGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTgg                                                                                    <  1:1407863/93‑1 (MQ=255)
     tGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTgg                                                                                    >  2:1368344/1‑93 (MQ=255)
      gcgcCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGgc                                                                                   <  2:317248/93‑1 (MQ=255)
       cgcCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGgcg                                                                                  <  2:3715803/93‑1 (MQ=255)
            gTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCa                                                                             >  1:1776578/1‑93 (MQ=255)
            gTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCa                                                                             >  1:3877668/1‑93 (MQ=255)
            gTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCa                                                                             >  1:2147869/1‑93 (MQ=255)
            gTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCa                                                                             >  2:2993863/1‑93 (MQ=255)
               cGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAggg                                                                          >  2:1395932/1‑93 (MQ=255)
                gagaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGc                                                                         >  1:2062884/1‑93 (MQ=255)
                gagaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGc                                                                         >  2:4144111/1‑93 (MQ=255)
                 agaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCa                                                                        <  1:594134/93‑1 (MQ=255)
                 agaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATGGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCa                                                                        <  1:594133/93‑1 (MQ=255)
                  gctGCTGGGGAGCTGTCCCCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCAt                                                                       <  2:3977796/91‑1 (MQ=255)
                  gaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCAt                                                                       >  1:3305973/1‑93 (MQ=255)
                  gaTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCAt                                                                       >  2:565071/1‑93 (MQ=255)
                   aTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCTCGTCAGGGCATc                                                                      <  1:3361766/93‑1 (MQ=255)
                   aTGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATc                                                                      <  1:821845/93‑1 (MQ=255)
                    tGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCg                                                                     >  1:1294769/1‑93 (MQ=255)
                       tGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcc                                                                  <  2:2560194/93‑1 (MQ=255)
                        ggggAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcct                                                                 <  2:3630973/93‑1 (MQ=255)
                        ggggAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcct                                                                 <  2:1421833/93‑1 (MQ=255)
                        ggggAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcct                                                                 >  1:3578400/1‑93 (MQ=255)
                            aGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcctgcct                                                             <  2:4563289/93‑1 (MQ=255)
                            aGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcctgcct                                                             <  1:68494/93‑1 (MQ=255)
                             gCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcctgcctg                                                            <  1:3480468/93‑1 (MQ=255)
                             gCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGAcctgcctg                                                            <  1:87078/93‑1 (MQ=255)
                                aTCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCa                                                         >  2:1854547/1‑93 (MQ=255)
                                  ccGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACg                                                       <  2:2604630/93‑1 (MQ=255)
                                  ccGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACg                                                       <  1:4748210/93‑1 (MQ=255)
                                  ccGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATACAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACg                                                       <  2:1642551/93‑1 (MQ=255)
                                    gCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACggg                                                     >  2:2027196/1‑93 (MQ=255)
                                    gCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACggg                                                     >  1:1856032/1‑93 (MQ=255)
                                      tCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCa                                                   >  1:4570873/1‑93 (MQ=255)
                                       cAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCAt                                                  >  1:4114488/1‑93 (MQ=255)
                                        aaTGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCAtt                                                 >  2:2657386/1‑93 (MQ=255)
                                        aaTGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCAtt                                                 >  1:2735866/1‑93 (MQ=255)
                                         aTGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTg                                                <  1:1393899/93‑1 (MQ=255)
                                            aGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTgcgc                                             >  2:4745983/1‑93 (MQ=255)
                                             gtcggtcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCa                                            >  2:163099/1‑93 (MQ=255)
                                             gtcggtcgCCATTGCCCGGGCGCGTGACAAATTGCGGTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCa                                            <  1:3259331/93‑1 (MQ=255)
                                                ggtcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTc                                         >  1:3377609/1‑93 (MQ=255)
                                                  tcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGc                                       <  2:127575/93‑1 (MQ=255)
                                                  tcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGc                                       >  2:868924/1‑93 (MQ=255)
                                                  tcgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGc                                       >  1:4075877/1‑93 (MQ=255)
                                                   cgCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGcc                                      <  1:3278051/93‑1 (MQ=255)
                                                      cATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGga                                   >  2:4426726/1‑93 (MQ=255)
                                                        ttGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGgatg                                 >  2:1381382/1‑93 (MQ=255)
                                                        ttGCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGgatg                                 <  2:21079/93‑1 (MQ=255)
                                                          gCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGgatgat                               <  2:2459007/93‑1 (MQ=255)
                                                          gCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGgatgat                               <  1:1530893/93‑1 (MQ=255)
                                                          gCCCGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGgatgat                               <  2:2459008/93‑1 (MQ=255)
                                                           cccGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATa                              >  1:2962873/1‑93 (MQ=255)
                                                             cGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATagc                            <  2:1174388/93‑3 (MQ=255)
                                                             cGGGCGCGTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATAcc                            <  1:2968014/93‑1 (MQ=255)
                                                              gggCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCACTCGCCGGATGATACCa                           >  1:4276444/1‑93 (MQ=255)
                                                                  gcgTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGa                       <  1:4071734/93‑1 (MQ=255)
                                                                  gcgTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGa                       >  1:3653080/1‑93 (MQ=255)
                                                                  gcgTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGa                       >  2:2795194/1‑93 (MQ=255)
                                                                   cgTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAt                      >  2:3593688/1‑93 (MQ=255)
                                                                    gacaCAAATTGCCTTCCATGCAGCTGCTGGCGCGGCAGGGCATCGCCCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     <  1:1932456/90‑1 (MQ=255)
                                                                    gTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     >  2:3070921/1‑93 (MQ=255)
                                                                    gTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     >  2:4092903/1‑93 (MQ=255)
                                                                    gTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     >  1:3790645/1‑93 (MQ=255)
                                                                    gTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     >  1:257398/1‑93 (MQ=255)
                                                                    gTGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAtt                     <  1:2082083/93‑1 (MQ=255)
                                                                     tGACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGAttt                    <  2:5295/93‑1 (MQ=255)
                                                                      gACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTa                   >  2:761746/1‑93 (MQ=255)
                                                                      gACAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTa                   <  1:635127/93‑1 (MQ=255)
                                                                       aCAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTaa                  >  1:929310/1‑93 (MQ=255)
                                                                       aCAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTaa                  >  1:3077211/1‑93 (MQ=255)
                                                                       aCAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTaa                  <  2:1094832/93‑1 (MQ=255)
                                                                       aCAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTaa                  >  1:4440053/1‑93 (MQ=255)
                                                                        cAAATTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAAt                 >  1:1517852/1‑93 (MQ=255)
                                                                         aaaaTGCGTTCCATGTAGCTAAGGCCGCGTTAGGGTATCGCCCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATc                <  2:3969947/93‑1 (MQ=255)
                                                                          aaTTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCg               >  1:1363624/1‑93 (MQ=255)
                                                                          aaTTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCg               >  2:4641805/1‑93 (MQ=255)
                                                                          aaTTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCg               >  2:3438819/1‑93 (MQ=255)
                                                                          aaTTGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCg               >  1:2896558/1‑93 (MQ=255)
                                                                            ttGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGAc             <  1:966508/93‑1 (MQ=255)
                                                                             tGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACa            <  2:1776467/93‑1 (MQ=255)
                                                                             tGCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACa            <  2:773740/93‑1 (MQ=255)
                                                                              gCGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACAt           <  2:2493494/93‑1 (MQ=255)
                                                                               cGTTCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATg          <  1:2839744/93‑1 (MQ=255)
                                                                                 ttCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGt        <  1:4025308/93‑1 (MQ=255)
                                                                                  tCCATGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTc       >  1:3929176/1‑93 (MQ=255)
                                                                                     aTGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCTGATGATACCAGCGATTTAATCGACATGGTCggt    <  2:2537628/93‑1 (MQ=255)
                                                                                     aTGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggt    <  1:2792219/93‑1 (MQ=255)
                                                                                      tGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggtg   >  1:2109993/1‑93 (MQ=255)
                                                                                      tGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggtg   >  1:38405/1‑93 (MQ=255)
                                                                                      tGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggtg   >  1:1915945/1‑93 (MQ=255)
                                                                                      tGCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggtg   <  2:2779504/93‑1 (MQ=255)
                                                                                       gCAGCTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCggtgg  <  2:1689050/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGCACTGCGCCAGTTCGAGATGCTGGGGAGCTATCCGCTCAATGAGTCGGTCGCCATTGCCCGGGCGCGTGACAAATTGCGTTCCATGCAACTGCTGGCGCGTCAGGGCATCGACCTGCCTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGG  >  NC_000913/891414‑891593

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