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,3390AG99.8% 1773.5 / NA 544E122G (GAA→GGA) yjgLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (272/271);  minor base C (1/0);  total (273/271)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 561 total aligned reads displayed.

ATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAT  >  NC_000913/4473249‑4473428
                                                                                          |                                                                                         
 tatGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACg                                                                                        >  2:9093180/1‑93 (MQ=255)
 tatGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACg                                                                                        <  2:6847879/93‑1 (MQ=255)
  atgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGa                                                                                       >  1:6593014/1‑93 (MQ=255)
   tgatgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAg                                                                                      <  1:9657970/93‑1 (MQ=255)
     atgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCg                                                                                    <  1:9845304/93‑1 (MQ=255)
      tgGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCgg                                                                                   <  1:7003211/93‑1 (MQ=255)
       ggatCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGa                                                                                  <  1:5729597/90‑1 (MQ=255)
       ggCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGa                                                                                  <  1:7511194/93‑1 (MQ=255)
         cTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGAtt                                                                                <  2:8482713/93‑1 (MQ=255)
          tCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTc                                                                               >  1:527484/1‑93 (MQ=255)
          tCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTc                                                                               >  2:2869299/1‑93 (MQ=255)
            atatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAc                                                                             <  2:5599904/93‑1 (MQ=255)
             tatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAcc                                                                            >  1:1107520/1‑93 (MQ=255)
             tatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAcc                                                                            <  1:7180732/93‑1 (MQ=255)
             tatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAcc                                                                            >  1:8412406/1‑93 (MQ=255)
             tatTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCAcc                                                                            >  2:3373774/1‑93 (MQ=255)
               ttAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCa                                                                          >  1:413113/1‑93 (MQ=255)
               ttAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCa                                                                          <  2:3864143/93‑1 (MQ=255)
                tAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCaa                                                                         <  2:1871785/93‑1 (MQ=255)
                tAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCaa                                                                         >  1:1405503/1‑93 (MQ=255)
                 aaGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAg                                                                        <  2:2665709/93‑1 (MQ=255)
                 aaGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAg                                                                        <  2:3602772/93‑1 (MQ=255)
                  aGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAgg                                                                       >  2:301596/1‑93 (MQ=255)
                  aGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAgg                                                                       <  1:1333885/93‑1 (MQ=255)
                     ccTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGa                                                                    >  2:2483046/1‑93 (MQ=255)
                     ccTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGa                                                                    >  2:10023864/1‑93 (MQ=255)
                     ccTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGa                                                                    >  2:9370313/1‑93 (MQ=255)
                     ccTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGa                                                                    <  2:1235200/93‑1 (MQ=255)
                        taTCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATta                                                                 <  2:4253874/92‑1 (MQ=255)
                         aTCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtat                                                                <  1:6747726/93‑1 (MQ=255)
                         aTCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtat                                                                <  1:1240131/93‑1 (MQ=255)
                          tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               >  1:3960727/1‑93 (MQ=255)
                          tCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATtata                                                               <  1:8212443/93‑1 (MQ=255)
                              aTTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGaat                                                           >  1:563690/1‑92 (MQ=255)
                               ttCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAAt                                                          >  2:6985272/1‑93 (MQ=255)
                                 cTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAAttt                                                        <  2:4855511/93‑1 (MQ=255)
                                     aGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACt                                                    <  2:1979100/93‑1 (MQ=255)
                                         aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTcc                                                >  1:1715387/1‑93 (MQ=255)
                                         aCTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTcc                                                <  1:4930554/93‑1 (MQ=255)
                                          cTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCa                                               >  1:7324241/1‑93 (MQ=255)
                                           tAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCAc                                              >  2:7740011/1‑93 (MQ=255)
                                             aTGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTa                                            >  2:7232824/1‑93 (MQ=255)
                                              tGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTAt                                           >  2:756383/1‑93 (MQ=255)
                                               tgCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACCATc                                          <  2:1332452/92‑1 (MQ=255)
                                               ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATc                                          >  1:5935406/1‑93 (MQ=255)
                                               ggCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATc                                          <  1:9210928/93‑1 (MQ=255)
                                                 cAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCac                                        <  2:4752851/93‑1 (MQ=255)
                                                    tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactac                                     >  1:932858/1‑93 (MQ=255)
                                                    tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactac                                     >  2:6607512/1‑93 (MQ=255)
                                                    tGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactac                                     <  1:3441221/93‑1 (MQ=255)
                                                     gAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGGCCGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCactact                                    >  2:5121908/1‑93 (MQ=255)
                                                      aaaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCActacta                                   <  1:100084/93‑1 (MQ=255)
                                                      aaaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCActacta                                   <  1:6180626/93‑1 (MQ=255)
                                                       aaaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTAt                                  >  2:3998515/1‑93 (MQ=255)
                                                        aaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATg                                 >  1:3721384/1‑93 (MQ=255)
                                                        aaaaTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCAAAGGTGATTATAGAAATTTCACTTCCCACTATCACTACTATg                                 <  2:3437482/93‑1 (MQ=255)
                                                           aTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATgggg                              >  1:196912/1‑93 (MQ=255)
                                                           aTGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATgggg                              >  1:8835526/1‑93 (MQ=255)
                                                            tGAAGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGa                             >  2:2166421/1‑93 (MQ=255)
                                                               aGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTc                          <  2:1398514/93‑1 (MQ=255)
                                                               aGTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTc                          >  2:6010796/1‑93 (MQ=255)
                                                                gTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCg                         >  1:956654/1‑93 (MQ=255)
                                                                gTTTTACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCg                         >  2:6344948/1‑93 (MQ=255)
                                                                  tttACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTc                       <  1:244826/93‑1 (MQ=255)
                                                                  tgtACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTc                       <  1:5532750/91‑1 (MQ=255)
                                                                   ttACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCa                      <  2:1776708/93‑1 (MQ=255)
                                                                   ttACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCa                      <  1:7784596/93‑1 (MQ=255)
                                                                    tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  2:8692561/93‑1 (MQ=255)
                                                                    tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  2:5034235/93‑1 (MQ=255)
                                                                    tACATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCaa                     <  2:5790208/93‑1 (MQ=255)
                                                                     aCATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                    <  2:4608425/93‑1 (MQ=255)
                                                                     aCATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAg                    <  1:4439059/93‑1 (MQ=255)
                                                                      cATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAgg                   <  1:2985199/93‑1 (MQ=255)
                                                                      cATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAgg                   <  2:6159148/93‑1 (MQ=255)
                                                                      cATGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAgg                   >  1:6046961/1‑93 (MQ=255)
                                                                       aggAATAATCCTGTAGGGGGACGATCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCTCTACTATGGGGACTCGTCAAGGa                  <  2:2064632/91‑1 (MQ=255)
                                                                       aTGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGa                  >  1:2236186/1‑93 (MQ=255)
                                                                        tGAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGAc                 >  1:4062212/1‑93 (MQ=255)
                                                                         gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                <  1:6405009/93‑1 (MQ=255)
                                                                         gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                >  1:4231148/1‑93 (MQ=255)
                                                                         gAATAATCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACa                <  2:1067561/93‑1 (MQ=255)
                                                                          aataatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACAt               <  1:2784679/93‑1 (MQ=255)
                                                                           ttaaaCCTGTAGGGGGGCGAGGGGGTTTACCCCAGGGGATTATAGAAAATTTACTTTCCCCTATCACTACTATGGGGACTCGTCAAGGAAATa              >  2:8203636/2‑93 (MQ=255)
                                                                            taatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATAc             <  2:2762512/93‑1 (MQ=255)
                                                                             aatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa            <  2:5328878/93‑1 (MQ=255)
                                                                             aatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa            <  2:2302438/93‑1 (MQ=255)
                                                                             aatCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACa            >  2:9723567/1‑93 (MQ=255)
                                                                              atCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAg           >  2:4318405/1‑93 (MQ=255)
                                                                               tCCTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGc          >  1:8651109/1‑93 (MQ=255)
                                                                                 cTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTGCTATGGGGACTCGTCAAGGACATACAGCCa        >  2:9822790/1‑93 (MQ=255)
                                                                                 cTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCa        >  1:3564867/1‑93 (MQ=255)
                                                                                 cTGTAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCa        >  1:4631474/1‑93 (MQ=255)
                                                                                    tAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta     >  2:7915012/1‑93 (MQ=255)
                                                                                    tAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta     >  2:6710322/1‑93 (MQ=255)
                                                                                    tAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta     <  2:8824351/93‑1 (MQ=255)
                                                                                    tAGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATta     <  1:8031429/93‑1 (MQ=255)
                                                                                     aGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat    <  2:8379697/93‑1 (MQ=255)
                                                                                     aGGGGGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtat    >  1:6467347/1‑93 (MQ=255)
                                                                                     aGGGGGACGAGCGGAATCACCCAAGGTGATTATAGAAATTTCACTATCCACTATCACTACTATAGGGACTCGTCAAGGACATACAGCCATtaa    <  1:8930240/93‑2 (MQ=255)
                                                                                      ggggGACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATtata   >  2:8033734/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAAAAATGAAGTTTTACATGAATAATCCTGTAGGGGAACGAGCGGATTCACCCAAGGTGATTATAGAAATTTCACTTTCCACTATCACTACTATGGGGACTCGTCAAGGACATACAGCCATTATAT  >  NC_000913/4473249‑4473428

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