Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,339 A→G E122G (GAA→GGA)  yjgL → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,3390AG98.1% 492.6 / ‑3.6 155E122G (GAA→GGA) yjgLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (68/84);  total (70/85)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAT  >  NC_000913/4473250‑4473428
                                                                                         |                                                                                         
tatGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACg                                                                                        >  2:1582944/1‑93 (MQ=255)
 atgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGa                                                                                       <  2:120359/93‑1 (MQ=255)
 atgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGa                                                                                       >  1:1789602/1‑93 (MQ=255)
  tgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAg                                                                                      >  2:4593803/1‑93 (MQ=255)
  tgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAg                                                                                      >  2:1697748/1‑93 (MQ=255)
    atgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCg                                                                                    <  2:2973724/93‑1 (MQ=255)
     tgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCgg                                                                                   >  2:4126462/1‑93 (MQ=255)
      ggCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGa                                                                                  <  1:1605477/93‑1 (MQ=255)
       gCTCTTATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCTATGACAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGAAt                                                                                 <  2:4226168/93‑1 (MQ=255)
        cTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGAtt                                                                                >  2:2355386/1‑93 (MQ=255)
          cATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAGAATGAAGTTTTACATGAATAATCGTGTAGGGGGACGAGCGAATTCa                                                                              >  1:4679558/1‑93 (MQ=255)
            tatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAcc                                                                            <  2:2964241/93‑1 (MQ=255)
             atTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAccc                                                                           >  1:4336544/1‑93 (MQ=255)
              ttAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCa                                                                          >  2:2507687/1‑93 (MQ=255)
               tAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCaa                                                                         <  2:2702305/93‑1 (MQ=255)
                aaGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAg                                                                        <  1:2231129/93‑1 (MQ=255)
                aaGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAg                                                                        <  2:4102593/93‑1 (MQ=255)
                 aGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAgg                                                                       >  1:3886931/1‑93 (MQ=255)
                 aGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAgg                                                                       <  1:1148583/93‑1 (MQ=255)
                     cTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGAt                                                                   <  2:547823/93‑1 (MQ=255)
                     cTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGAt                                                                   >  1:1861743/1‑93 (MQ=255)
                     cTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAAGGAAGTTTTACATGAATAAACCTGTAGGGGGACGAGCGGATTCACCCAAGGTGAt                                                                   <  2:4054929/93‑1 (MQ=255)
                      tCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGAtt                                                                  >  2:1813791/1‑93 (MQ=255)
                      tCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGAtt                                                                  >  1:3900101/1‑93 (MQ=255)
                       cATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATta                                                                 <  1:630056/93‑1 (MQ=255)
                       cATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATta                                                                 >  1:3848102/1‑93 (MQ=255)
                         tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               <  1:1045853/93‑1 (MQ=255)
                         tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               >  2:3378376/1‑93 (MQ=255)
                         tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               <  1:2610366/93‑1 (MQ=255)
                         tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               >  1:2929396/1‑93 (MQ=255)
                         tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               <  1:1948823/93‑1 (MQ=255)
                          cTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAg                                                              <  2:2219046/93‑1 (MQ=255)
                             aTTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGaaa                                                           >  2:4731017/1‑93 (MQ=255)
                              ttCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAAt                                                          >  2:2289795/1‑93 (MQ=255)
                               tCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAAtt                                                         >  1:2232300/1‑93 (MQ=255)
                                cTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAAttt                                                        <  1:3069391/93‑1 (MQ=255)
                                   aaGATACTAATGGCACTGACAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCAc                                                     <  2:4203974/93‑1 (MQ=255)
                                   aaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCAc                                                     <  2:3536578/93‑1 (MQ=255)
                                   aaGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCAc                                                     <  2:1667343/93‑1 (MQ=255)
                                    aGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACt                                                    >  1:2502416/1‑93 (MQ=255)
                                     gATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACtt                                                   <  1:1243944/93‑1 (MQ=255)
                                     gATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACtt                                                   <  1:3652833/93‑1 (MQ=255)
                                     gATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACtt                                                   <  1:1156478/93‑1 (MQ=255)
                                     gATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACtt                                                   >  2:4350572/1‑93 (MQ=255)
                                      aTACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACttt                                                  <  1:2097876/93‑1 (MQ=255)
                                       tACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTc                                                 >  2:2765251/1‑93 (MQ=255)
                                        aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTcc                                                >  1:4284227/1‑93 (MQ=255)
                                        aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTcc                                                <  1:3664930/93‑1 (MQ=255)
                                         cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCa                                               >  2:1633393/1‑93 (MQ=255)
                                          tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCAc                                              <  1:1120997/93‑1 (MQ=255)
                                            aTGGCAATGAAAAAATTAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTa                                            >  1:2362935/1‑93 (MQ=255)
                                             tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTAt                                           <  2:1908659/93‑1 (MQ=255)
                                             tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTAt                                           >  1:94254/1‑93 (MQ=255)
                                             tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTAt                                           <  2:1432641/93‑1 (MQ=255)
                                              ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATc                                          <  1:4554162/93‑1 (MQ=255)
                                                cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCac                                        <  2:3617973/93‑1 (MQ=255)
                                                 aaTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCact                                       >  2:1694538/1‑93 (MQ=255)
                                                 aaTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCact                                       <  1:709320/93‑1 (MQ=255)
                                                  aTGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCacta                                      >  2:4678206/1‑93 (MQ=255)
                                                   tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactac                                     >  2:1278001/1‑93 (MQ=255)
                                                   agAAAAAATGAAGTCTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactac                                     <  1:1898440/92‑1 (MQ=255)
                                                     aaaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCActacta                                   <  1:3850358/93‑1 (MQ=255)
                                                      aaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTAt                                  <  1:2678821/93‑1 (MQ=255)
                                                        aaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATgg                                >  1:3918795/1‑93 (MQ=255)
                                                         aaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATggg                               <  2:522301/93‑1 (MQ=255)
                                                         aaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATggg                               >  2:2861632/1‑93 (MQ=255)
                                                          aTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCGCCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATgggg                              >  1:1191727/1‑93 (MQ=255)
                                                          aTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATgggg                              <  1:498076/93‑1 (MQ=255)
                                                           tGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGa                             >  1:1386015/1‑93 (MQ=255)
                                                           tGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGa                             >  1:1386016/1‑93 (MQ=255)
                                                           agAAGTTTTAAATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCCCTTTCCACTATCACTACTATGGGGa                             <  1:31106/92‑1 (MQ=255)
                                                            gAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGAc                            <  1:2386535/93‑1 (MQ=255)
                                                              aGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTc                          <  1:3279422/93‑1 (MQ=255)
                                                                  ttACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCa                      >  2:4311677/1‑93 (MQ=255)
                                                                   tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  1:4112629/93‑1 (MQ=255)
                                                                   tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  2:2599530/93‑1 (MQ=255)
                                                                   tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  1:3247422/93‑1 (MQ=255)
                                                                   tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  1:3779342/93‑1 (MQ=255)
                                                                    aCATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                    >  2:3235331/1‑93 (MQ=255)
                                                                      aTGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGa                  <  2:1323378/93‑1 (MQ=255)
                                                                       tGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTAGTCAAGGAc                 >  2:2694352/1‑93 (MQ=255)
                                                                       tGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACCTTCCACTATCACTACTATGGGGACTCGTCAAGGAc                 <  2:3977797/93‑1 (MQ=255)
                                                                       agAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGAc                 <  1:4189133/92‑1 (MQ=255)
                                                                        gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                >  2:4400326/1‑93 (MQ=255)
                                                                        gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                <  1:1559702/93‑1 (MQ=255)
                                                                        gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                <  2:3694838/93‑1 (MQ=255)
                                                                        gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                >  2:207790/1‑93 (MQ=255)
                                                                         aataatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACAt               <  2:778223/93‑1 (MQ=255)
                                                                          ataatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATa              >  1:1070540/1‑93 (MQ=255)
                                                                           taatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATAc             >  1:3462664/1‑93 (MQ=255)
                                                                            aatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa            <  2:3779754/93‑1 (MQ=255)
                                                                              tCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc          >  2:4593534/1‑93 (MQ=255)
                                                                              tCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc          >  1:4099074/1‑93 (MQ=255)
                                                                              tCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc          >  1:2997330/1‑93 (MQ=255)
                                                                               ccTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGcc         <  2:3251532/93‑1 (MQ=255)
                                                                                 tGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCAt       >  2:4727148/1‑93 (MQ=255)
                                                                                    aGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat    <  1:1944182/93‑1 (MQ=255)
                                                                                     ggggGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtata   >  1:4273536/1‑93 (MQ=255)
                                                                                     ggggAACGAGCGGGTTCACCCAAGGTGATTATAGAAATTTCACTTTTCACTATCACTACTATCGGGACTCGTCAAAGACATACCGCCATtata   >  2:3861676/1‑93 (MQ=255)
                                                                                      ggggACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTatat  >  2:2780055/1‑93 (MQ=255)
                                                                                         |                                                                                         
TATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAT  >  NC_000913/4473250‑4473428

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