Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,017,557 A→G E116G (GAG→GGG)  yigP → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,5570AG99.4% 1126.4 / ‑5.0 356E116G (GAG→GGG) yigPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (148/206);  total (148/208)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 376 total aligned reads displayed.

AGCTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGAGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTCC  >  NC_000913/4017467‑4017646
                                                                                          |                                                                                         
  cTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTc                                                                                       <  2:3092346/93‑1 (MQ=255)
  cTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTc                                                                                       >  1:5926175/1‑93 (MQ=255)
  cTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTc                                                                                       >  2:8455116/1‑93 (MQ=255)
   ttACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCg                                                                                      <  2:9305462/93‑1 (MQ=255)
   ttACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCg                                                                                      <  1:3448787/93‑1 (MQ=255)
    tACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGa                                                                                     <  2:207412/93‑1 (MQ=255)
     tccGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAc                                                                                    <  1:8837900/92‑1 (MQ=255)
     aCCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAc                                                                                    <  2:8352004/93‑1 (MQ=255)
     aCCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAc                                                                                    <  2:67903/93‑1 (MQ=255)
     aCCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAc                                                                                    >  1:1222694/1‑93 (MQ=255)
       cGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAccc                                                                                  <  2:559803/93‑1 (MQ=255)
       cGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGAccc                                                                                  <  1:3470858/93‑1 (MQ=255)
         cACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTg                                                                                >  1:10000329/1‑93 (MQ=255)
         cACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTg                                                                                >  1:3868800/1‑93 (MQ=255)
          aCTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGTGATATTCAGGTGGTCCAACTCTTCGTAGCGCGGGCAGATCTGGCAGGGTTCGGCCCTGc                                                                               <  2:4411232/93‑1 (MQ=255)
           cTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCg                                                                              >  1:5804369/1‑93 (MQ=255)
            tGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCgg                                                                             >  1:624768/1‑93 (MQ=255)
               ttCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAc                                                                          >  1:1669905/1‑93 (MQ=255)
               ttCGCAGTGATGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAc                                                                          <  2:2439180/93‑1 (MQ=255)
                tCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAct                                                                         >  1:4021336/1‑93 (MQ=255)
                tCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAAct                                                                         >  2:5918430/1‑93 (MQ=255)
                   cAGTGGTGAGATGGAGGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCGGATCTGGCAGGGTTCGACCCTGCGGAActgct                                                                      <  2:8642067/93‑1 (MQ=255)
                     gtggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTgg                                                                    >  1:2607279/1‑93 (MQ=255)
                     gtggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTgg                                                                    >  1:9801421/1‑93 (MQ=255)
                     gtggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTgg                                                                    <  1:8205419/93‑1 (MQ=255)
                      tggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGc                                                                   >  2:3376544/1‑93 (MQ=255)
                      tggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGc                                                                   >  2:2715954/1‑93 (MQ=255)
                      tggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGc                                                                   >  1:8919040/1‑93 (MQ=255)
                      tggtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGCGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGc                                                                   >  2:3677999/1‑93 (MQ=255)
                        gtgAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGccc                                                                 <  1:9638572/93‑1 (MQ=255)
                        gtgAGCTGGAAGTGCAGGGCGATATTCAGGGGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGccc                                                                 >  1:1564985/1‑93 (MQ=255)
                           aGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCtt                                                              <  2:2488962/93‑1 (MQ=255)
                           aGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCtt                                                              <  2:536474/93‑1 (MQ=255)
                            ggtGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTta                                                             <  1:485020/91‑1 (MQ=255)
                            gCTGGAAGTGCAGGGCGATATTCAGGTTGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGCTTCGACCCTGCGGAACTGCTGCCCCCTta                                                             <  2:1749172/93‑1 (MQ=255)
                            gCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTta                                                             <  2:9446953/93‑1 (MQ=255)
                              tGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGTCCCCTtata                                                           <  2:9895099/93‑1 (MQ=255)
                              tGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTtata                                                           <  1:834849/93‑1 (MQ=255)
                              gggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGAGTTCGACCCTGCGGAACTGCTGGCCCCTtata                                                           <  1:6034316/92‑1 (MQ=255)
                               ggAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATAc                                                          >  1:9130204/1‑93 (MQ=255)
                                  aGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCgg                                                       <  2:5013753/93‑1 (MQ=255)
                                   gTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGt                                                      >  1:2339091/1‑93 (MQ=255)
                                   gTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGt                                                      >  2:724469/1‑93 (MQ=255)
                                     gCAGGGCGATATTCAGGTGGTGCAAAGCTTCGTTGCGTTGGCAGATCTGGCATGTTTCGACCCTGCGGAACTGCTGGGACCTTATACCGGTGa                                                    <  2:6310723/93‑1 (MQ=255)
                                     gCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGa                                                    >  1:3734682/1‑93 (MQ=255)
                                     gCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGa                                                    <  1:4996155/93‑1 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                 >  2:8868542/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                 >  1:2939930/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                 <  1:9443558/93‑1 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                 >  2:6655135/1‑93 (MQ=255)
                                        gggCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGatat                                                 <  1:3550274/93‑1 (MQ=255)
                                            gATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCATATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgcc                                             <  2:4589882/93‑1 (MQ=255)
                                            gATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgcc                                             >  2:4910804/1‑93 (MQ=255)
                                             atatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccg                                            <  1:8754285/93‑1 (MQ=255)
                                             atatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccg                                            >  1:7553430/1‑93 (MQ=255)
                                              tatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccgc                                           <  1:5659645/93‑1 (MQ=255)
                                              tatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccgc                                           >  2:6868572/1‑93 (MQ=255)
                                              tatTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATcgccgc                                           <  1:4351495/93‑1 (MQ=255)
                                                ttCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTg                                         <  2:4526226/93‑1 (MQ=255)
                                                 tCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGa                                        <  2:6015849/93‑1 (MQ=255)
                                                 tCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGa                                        <  2:3842502/93‑1 (MQ=255)
                                                  cAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGACCTGCTGGCCCCTTATACCGGTGATATCGCCGCTGaa                                       <  1:1428998/93‑1 (MQ=255)
                                                  cAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGaa                                       <  2:8181593/93‑1 (MQ=255)
                                                  cAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGaa                                       <  2:6963394/93‑1 (MQ=255)
                                                  cAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGaa                                       <  2:1622656/93‑1 (MQ=255)
                                                   aGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAg                                      <  2:5164407/93‑1 (MQ=255)
                                                        gtgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATc                                 <  1:659385/93‑1 (MQ=255)
                                                        gtgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATc                                 <  2:6213873/93‑1 (MQ=255)
                                                         tgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCa                                <  2:9657258/93‑1 (MQ=255)
                                                         tgCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCa                                <  2:1299323/93‑1 (MQ=255)
                                                          gCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAg                               <  1:5308003/93‑1 (MQ=255)
                                                          gCAAAACTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAg                               >  1:3000941/1‑93 (MQ=255)
                                                             aaaCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCaa                            <  1:7134287/93‑1 (MQ=255)
                                                             aaaCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCaa                            >  1:6813974/1‑93 (MQ=255)
                                                               aCTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAg                          <  2:2575159/93‑1 (MQ=255)
                                                                cTTCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGc                         >  1:6611101/1‑93 (MQ=255)
                                                                 ttCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGcc                        <  1:3281505/93‑1 (MQ=255)
                                                                  tCGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCa                       <  1:8544607/93‑1 (MQ=255)
                                                                   cGTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCAt                      <  2:5727408/93‑1 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                     <  2:1988742/93‑1 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                     <  2:2974462/93‑1 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                     >  2:6589579/1‑93 (MQ=255)
                                                                    gTTGCGCTGGCAGATCTGGCAGGGTGCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                     <  2:7411534/93‑1 (MQ=255)
                                                                    gTTGCGCTGGCAGAGCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATCCCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATg                     <  2:4477086/93‑1 (MQ=255)
                                                                       gcgcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATgcgc                  <  1:8701205/93‑1 (MQ=255)
                                                                        cgcTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGcgcg                 <  2:4640740/93‑1 (MQ=255)
                                                                          cTGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGTCCCCTTATACCGGTGATCTCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGa               <  1:1827519/93‑1 (MQ=255)
                                                                           tGGCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAg              >  2:1465385/1‑93 (MQ=255)
                                                                            ggCAGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAgg             <  2:4337302/93‑1 (MQ=255)
                                                                               aGATCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGgcgc          >  2:7989166/1‑93 (MQ=255)
                                                                                  tCTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCaaa       <  1:8314654/93‑1 (MQ=255)
                                                                                   cTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAg      >  2:76832/1‑93 (MQ=255)
                                                                                   cTGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAg      >  2:1851418/1‑93 (MQ=255)
                                                                                    tGGCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGt     <  2:9161924/93‑1 (MQ=255)
                                                                                     ggCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGtt    >  1:537108/1‑93 (MQ=255)
                                                                                      gCAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTc   >  1:4556691/1‑93 (MQ=255)
                                                                                       cAGGGTTCGACCCTGCGGAACTGCTGGCGCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTcc  <  1:7966239/93‑1 (MQ=255)
                                                                                       cAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTcc  <  1:209185/93‑1 (MQ=255)
                                                                                       cAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTcc  <  2:7083385/93‑1 (MQ=255)
                                                                                       cAGGGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTcc  <  2:6559088/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGCTTACCGCACTGATTCGCAGTGGTGAGCTGGAAGTGCAGGGCGATATTCAGGTGGTGCAAAACTTCGTTGCGCTGGCAGATCTGGCAGAGTTCGACCCTGCGGAACTGCTGGCCCCTTATACCGGTGATATCGCCGCTGAAGGAATCAGCAAAGCCATGCGCGGAGGCGCAAAGTTCC  >  NC_000913/4017467‑4017646

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