Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,248,121 T→C intergenic (+68/‑161) malM → / → b4559 MalM/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,248,1210TC83.8% 540.7 / 101.0 241intergenic (+68/‑161)malM/b4559MalM/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base C (98/104);  total (117/124)
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.99e-01

Only 100 of 262 total aligned reads displayed.

AGCAGTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCTTTTCCT  >  NC_000913/4248030‑4248212
                                                                                           |                                                                                           
 gcagTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGcttccctc                                                                                           >  1:2033318/1‑93 (MQ=255)
  cagTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCttccttcc                                                                                          <  2:7497969/93‑1 (MQ=255)
  cagTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCttccctcc                                                                                          <  2:5975117/93‑1 (MQ=255)
  cagTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCttccctcc                                                                                          >  2:5436813/1‑93 (MQ=255)
  cagTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCttccctcc                                                                                          <  2:5067635/93‑1 (MQ=255)
   agTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTtccctcct                                                                                         >  2:2733131/1‑93 (MQ=255)
    gtgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACTTTTTTTTATGCTTCCCTCCTg                                                                                        <  2:5295357/93‑1 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTgg                                                                                       >  2:6080431/1‑93 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTgg                                                                                       >  2:35707/1‑93 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTgg                                                                                       >  1:7448131/1‑93 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTgg                                                                                       >  1:958706/1‑93 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTgg                                                                                       >  1:32808/1‑93 (MQ=255)
     tgtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTgg                                                                                       >  2:4826478/1‑93 (MQ=255)
      gtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTggg                                                                                      >  2:1952226/1‑93 (MQ=255)
      gtAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTggg                                                                                      >  2:6886741/1‑93 (MQ=255)
           aGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATg                                                                                 <  2:1417983/93‑1 (MQ=255)
            ggCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGa                                                                                <  1:4346442/93‑1 (MQ=255)
              cAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGc                                                                              <  2:4502438/93‑1 (MQ=255)
              cAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGc                                                                              <  1:204140/93‑1 (MQ=255)
               aaGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCg                                                                             >  2:5586380/1‑93 (MQ=255)
                 ggggTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGAt                                                                           >  1:49080/1‑93 (MQ=255)
                  gggTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGAtt                                                                          <  1:287022/93‑1 (MQ=255)
                    gTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGAtttt                                                                        >  2:1329888/1‑93 (MQ=255)
                     tAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGAttttt                                                                       >  1:4121829/1‑93 (MQ=255)
                     tAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGAttttt                                                                       >  1:4633502/1‑93 (MQ=255)
                       aTTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTta                                                                     <  1:5964706/93‑1 (MQ=255)
                           cGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGt                                                                 <  2:5251567/93‑1 (MQ=255)
                              cccACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAAc                                                              >  2:871199/1‑93 (MQ=255)
                               ccACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACt                                                             >  1:3527810/1‑93 (MQ=255)
                               ccACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACt                                                             >  1:6748659/1‑93 (MQ=255)
                                cacaGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTc                                                            >  1:3357070/1‑93 (MQ=255)
                                   aGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCAct                                                         <  2:7484568/93‑1 (MQ=255)
                                      gCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtct                                                      >  1:1469074/1‑93 (MQ=255)
                                      gCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtct                                                      <  1:3775422/93‑1 (MQ=255)
                                       cTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtctt                                                     >  2:5856609/1‑93 (MQ=255)
                                       cTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtctt                                                     <  2:1531337/93‑1 (MQ=255)
                                        tGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttc                                                    >  1:4226677/1‑93 (MQ=255)
                                         gATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttct                                                   >  1:7063186/1‑93 (MQ=255)
                                         gATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttct                                                   >  2:6576252/1‑93 (MQ=255)
                                          aTTTTGCAACAACTGGTGCGTCTCCTGGCGCTCCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttctt                                                  <  1:2365262/93‑1 (MQ=255)
                                          aTTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttctt                                                  >  2:1629721/1‑93 (MQ=255)
                                          aTTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCActtcttctt                                                  <  1:5755585/93‑1 (MQ=255)
                                            tttGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCa                                                <  1:2283060/93‑1 (MQ=255)
                                              tGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACt                                              <  1:4887945/93‑1 (MQ=255)
                                               gCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTa                                             <  1:7128346/93‑1 (MQ=255)
                                               gCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTa                                             >  1:6262835/1‑93 (MQ=255)
                                                caacaaCTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTaa                                            <  2:4428620/93‑1 (MQ=255)
                                                 aacaacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAg                                           >  2:6838451/1‑93 (MQ=255)
                                                 aacaacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAg                                           <  2:7010024/93‑1 (MQ=255)
                                                  acaacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGa                                          <  2:1252837/93‑1 (MQ=255)
                                                   caacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGaa                                         <  2:5837026/93‑1 (MQ=255)
                                                    aacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAat                                        <  1:1193851/93‑1 (MQ=255)
                                                    aacTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAat                                        >  2:4984668/1‑93 (MQ=255)
                                                     acTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAata                                       <  2:5704078/93‑1 (MQ=255)
                                                     acTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAata                                       <  1:523685/93‑1 (MQ=255)
                                                     acTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAata                                       >  1:1024966/1‑93 (MQ=255)
                                                     acTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAata                                       >  2:560979/1‑93 (MQ=255)
                                                      cTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAatat                                      <  1:1835438/93‑1 (MQ=255)
                                                       tGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATc                                     >  2:3248175/1‑93 (MQ=255)
                                                       tGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATc                                     <  2:7339969/93‑1 (MQ=255)
                                                        ggTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATcc                                    <  1:550495/93‑1 (MQ=255)
                                                        ggTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATcc                                    <  2:4695072/93‑1 (MQ=255)
                                                         gTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCa                                   <  2:1760830/93‑1 (MQ=255)
                                                         gTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCa                                   <  1:6140435/93‑1 (MQ=255)
                                                          tGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCAt                                  <  1:1925733/93‑1 (MQ=255)
                                                           gCGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCAtt                                 <  1:4463301/93‑1 (MQ=255)
                                                            cGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCAttg                                <  2:3510837/93‑2 (MQ=255)
                                                            cGTCTCCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTa                                <  1:5904938/93‑1 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              >  2:3499400/1‑93 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              <  2:2599576/93‑1 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              <  1:1273240/93‑1 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              <  1:5722698/93‑1 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              <  1:6266784/93‑1 (MQ=255)
                                                              tctcCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtc                              <  2:3420862/93‑1 (MQ=255)
                                                               ctcCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtct                             >  2:3607986/1‑93 (MQ=255)
                                                                tcCTGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTAtctc                            <  1:1744222/93‑1 (MQ=255)
                                                                  cTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCaa                          <  2:3025715/93‑1 (MQ=255)
                                                                   tGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAAt                         <  1:2766243/93‑1 (MQ=255)
                                                                   tGGCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAAt                         <  2:4472093/93‑1 (MQ=255)
                                                                    ggCGCACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATg                        <  2:223654/93‑1 (MQ=255)
                                                                     gcgcACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGc                       <  1:4956786/93‑1 (MQ=255)
                                                                     gcgcACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGc                       >  1:1630438/1‑93 (MQ=255)
                                                                     gcgcACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGc                       <  2:6671263/93‑1 (MQ=255)
                                                                       gcACCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCt                     <  2:1044773/93‑1 (MQ=255)
                                                                         aCCTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTa                   >  2:6177757/1‑93 (MQ=255)
                                                                          ccTTTTTTTATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCttnn                  >  2:419883/1‑91 (MQ=255)
                                                                               ttttATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCagag             <  1:1000860/93‑1 (MQ=255)
                                                                               ttttATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCagag             <  1:1587407/93‑1 (MQ=255)
                                                                                 ttATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGAt           <  1:7237007/93‑1 (MQ=255)
                                                                                  tATGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGAtt          <  1:3919510/93‑1 (MQ=255)
                                                                                    tGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCt        >  2:3642675/1‑93 (MQ=255)
                                                                                    tGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCt        >  2:4603557/1‑93 (MQ=255)
                                                                                    tGCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCt        <  1:2468719/93‑1 (MQ=255)
                                                                                     gCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCtt       <  1:2570471/93‑1 (MQ=255)
                                                                                     gCTTCCCTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCtt       >  1:847514/1‑93 (MQ=255)
                                                                                      cttccctcCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCttt      <  1:4780930/93‑1 (MQ=255)
                                                                                       ttccttccTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCtttt     <  2:3271291/93‑1 (MQ=255)
                                                                                        tccctcctGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCTTTTc    <  2:6381021/93‑1 (MQ=255)
                                                                                        tccctcctGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCTTTTc    >  2:7190700/1‑93 (MQ=255)
                                                                                          cttcctGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCTTTTCCt  >  2:4130430/1‑93 (MQ=255)
                                                                                           |                                                                                           
AGCAGTGTAAAAGGCAAGGGGTAATTACGCCCCACAGTGCTGATTTTGCAACAACTGGTGCGTCTCCTGGCGCACCTTTTTTTATGCTTCCTTCCTGGGATATGAGCGATTTTTTATAGTAACTCACTTCTTCTTCACTAAGAATATCCATTATCTCAATGCCTTATCAGAGATTCTTTTCCT  >  NC_000913/4248030‑4248212

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