Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,498,550 (C)5→6 intergenic (‑101/‑278) yjgY ← / → yjgZ hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,498,5451.C98.0% 799.3 / 7.1 254intergenic (‑96/‑283)yjgY/yjgZhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base C (123/126);  total (124/130)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 284 total aligned reads displayed.

ACACTGGTACTGCAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTATGCATG  >  NC_000913/4498453‑4498632
                                                                                               |                                                                                         
   cTGGTACTGCAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCAccc                                                                                       <  1:662966/93‑1 (MQ=255)
   cTGGTACTGCAAGAGGGGGACTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCAccc                                                                                       >  2:6153853/1‑93 (MQ=255)
    tGGTACTGCAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCAcccc                                                                                      >  1:3077532/1‑93 (MQ=255)
        aCTGCAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCAcccccctc                                                                                 <  1:1069049/93‑4 (MQ=255)
            cAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGa                                                                              <  1:2890783/93‑1 (MQ=255)
              agagGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAgg                                                                            <  1:5682639/93‑1 (MQ=255)
               gagGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGc                                                                           >  2:4069384/1‑93 (MQ=255)
                 gggggCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCgt                                                                         >  2:6202348/1‑93 (MQ=255)
                 gggggCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCgt                                                                         >  2:5503880/1‑93 (MQ=255)
                 gggggCCTTCGCTCCTCTGACATTACTAAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCgt                                                                         <  1:2623706/93‑1 (MQ=255)
                  ggggCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCCCCCCCCTCGCGAGGCgtt                                                                        <  2:3469995/93‑1 (MQ=255)
                    ggCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCgttgt                                                                      <  2:3936162/93‑1 (MQ=255)
                      ccTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGttgttg                                                                    >  1:6147734/1‑93 (MQ=255)
                       cTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCTATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGc                                                                   <  2:1384350/93‑1 (MQ=255)
                       cTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGc                                                                   >  2:2439284/1‑93 (MQ=255)
                       cTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGc                                                                   <  2:3575050/93‑1 (MQ=255)
                       cTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGc                                                                   <  1:2321843/93‑1 (MQ=255)
                        ttCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGAGTTGTTGcc                                                                  >  2:3108773/1‑93 (MQ=255)
                            ctcctcTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTc                                                              >  1:4829340/1‑93 (MQ=255)
                             tcctcTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCa                                                             >  2:2602935/1‑93 (MQ=255)
                              cctcTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCAc                                                            <  2:5695792/93‑1 (MQ=255)
                              cctcTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCAc                                                            >  1:1611576/1‑93 (MQ=255)
                                tctGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCAccc                                                          >  2:4591408/1‑93 (MQ=255)
                                 ctGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCg                                                         >  2:4137331/1‑93 (MQ=255)
                                 ctGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACACg                                                         >  1:699016/1‑93 (MQ=255)
                                  tGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCgg                                                        >  2:2055986/1‑93 (MQ=255)
                                       ttACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGgacaga                                                   >  1:5519069/1‑93 (MQ=255)
                                        tACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGgacagac                                                  >  1:3170840/1‑93 (MQ=255)
                                        tACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGgacagac                                                  <  2:1520072/93‑1 (MQ=255)
                                          cTCAGTATTTTCACTACTACCCCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGacac                                                <  1:245266/93‑1 (MQ=255)
                                          cTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGacac                                                <  1:1374322/93‑1 (MQ=255)
                                           tCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACAcc                                               >  2:2885363/1‑93 (MQ=255)
                                             aGTATTTTCACTACTACCCCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGt                                             <  1:1782035/93‑1 (MQ=255)
                                             aGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGt                                             <  1:4724476/93‑1 (MQ=255)
                                              gTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGtt                                            <  2:5896378/93‑1 (MQ=255)
                                               tATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTa                                           <  2:3338676/93‑1 (MQ=255)
                                               tATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTa                                           <  1:4264938/93‑1 (MQ=255)
                                               tATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTa                                           <  1:2258013/93‑1 (MQ=255)
                                                aTTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCGCCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTAt                                          <  1:1734797/93‑1 (MQ=255)
                                                aTTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTAt                                          <  2:1457934/93‑1 (MQ=255)
                                                aTTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTAt                                          <  2:657817/93‑1 (MQ=255)
                                                      actactACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACtt                                    >  2:3765035/1‑93 (MQ=255)
                                                      actactACACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACtt                                    >  1:541313/1‑93 (MQ=255)
                                                        tactacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTtct                                  >  2:5788735/1‑93 (MQ=255)
                                                         cctacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTtctc                                 <  1:4418950/92‑1 (MQ=255)
                                                          ctacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTctct                                >  1:3578993/1‑93 (MQ=255)
                                                          ctacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTctct                                >  1:5189523/1‑93 (MQ=255)
                                                          ctacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTctct                                <  1:163380/93‑1 (MQ=255)
                                                           tacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCtt                               <  1:3216133/93‑1 (MQ=255)
                                                           tacACCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCtt                               <  1:5491459/93‑1 (MQ=255)
                                                            acacCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCttg                              >  1:1327943/1‑93 (MQ=255)
                                                            acacCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCttg                              >  1:2810809/1‑93 (MQ=255)
                                                            acacCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGGCACCCGGACAGACACCGTTATGGACTTCTCttt                              >  1:4588080/1‑92 (MQ=255)
                                                             cacCCCATCCGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCttgt                             <  2:5993707/93‑1 (MQ=255)
                                                             cacCCCAGCTGCCATCCATAACAACCGATTCACCACCCTCGCGAGGCGTTGTTGCCCGTCGCCCGGACAGACACCGTTATGGACTTCTCttgt                             <  1:3927668/93‑1 (MQ=255)
                                                              ctcccATTTGCCATCCACAACAACCGATTCACCCCGCTCGCGAGGCGTTGTTGCCCGTCCCCCGGACAGACACCGTTATGGACTTCTCttgtt                            <  1:2403204/91‑1 (MQ=255)
                                                              acCCCATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCttgtt                            >  1:150026/1‑93 (MQ=255)
                                                                 ccATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCa                         <  2:4836089/93‑1 (MQ=255)
                                                                  cATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAc                        <  2:5373334/93‑1 (MQ=255)
                                                                  cATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAc                        <  1:1064568/93‑1 (MQ=255)
                                                                  cATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAc                        <  2:4738934/93‑1 (MQ=255)
                                                                  cATCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAc                        <  2:3659456/93‑1 (MQ=255)
                                                                   aTCTGCCATCCATAACAACCGATTCGCCCCCCTCGCGAGGCGTTGTTGCCCGTCCCCCGGACAGACCCCGTTATGGACTTCTCTTGTTGCAcg                       <  1:1910894/93‑1 (MQ=255)
                                                                   aTCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcg                       <  2:4938366/93‑1 (MQ=255)
                                                                   aTCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcg                       >  1:4458175/1‑93 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgt                      >  1:947164/1‑92 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgt                      >  1:5075549/1‑92 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcga                      <  2:1890062/93‑1 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcga                      >  2:4491862/1‑93 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcga                      >  2:1784254/1‑93 (MQ=255)
                                                                    tCTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcga                      >  1:4756286/1‑93 (MQ=255)
                                                                     cTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgat                     <  1:3667993/93‑1 (MQ=255)
                                                                     cTGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgat                     <  2:1946255/93‑1 (MQ=255)
                                                                      tGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgatc                    <  1:4700245/93‑1 (MQ=255)
                                                                      tGCCATCCATAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgatc                    <  1:4956326/93‑1 (MQ=255)
                                                                         catccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACAGCGTTATGGACTTCTCTTGTTGCAcgatccat                 >  1:3450875/1‑93 (MQ=255)
                                                                         catccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcggtcgtt                 >  1:6243501/1‑91 (MQ=255)
                                                                         catccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgatcgat                 >  2:2620813/1‑93 (MQ=255)
                                                                         catccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgatcgat                 >  2:4687005/1‑93 (MQ=255)
                                                                         catccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCAcgatcgat                 <  2:5534189/93‑1 (MQ=255)
                                                                          atccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACgatcgatc                <  2:820922/93‑1 (MQ=255)
                                                                          atccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACgatcgatc                <  1:2282612/93‑1 (MQ=255)
                                                                          atccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACgatcgatc                >  1:713257/1‑93 (MQ=255)
                                                                           tccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCa               >  1:601144/1‑93 (MQ=255)
                                                                            ccatAACAGCCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCaa              <  2:4393459/93‑1 (MQ=255)
                                                                            ccatAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCaa              <  1:2177295/93‑1 (MQ=255)
                                                                             catAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAAt             <  2:2374665/93‑1 (MQ=255)
                                                                             catAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAAt             <  2:1410408/93‑1 (MQ=255)
                                                                             catAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAAt             >  1:4799283/1‑93 (MQ=255)
                                                                             catAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAAt             <  1:4143440/93‑1 (MQ=255)
                                                                              atAACAACCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAAtt            <  1:2000778/93‑1 (MQ=255)
                                                                                 acaacCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACCGACACCGTTATGGACTTCTCTTGTTGCACGAGCGATCAATTGtt         >  2:2801647/1‑93 (MQ=255)
                                                                                 acaacCGATTCACCCCCCTCGCGAGGCGGTGTTGCCCGTCACCCGGACAGACACCCCTATGGACTTCTCTTGTTGGACGATCGTTCCattatt         >  2:4299738/1‑93 (MQ=255)
                                                                                  caacCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTa        <  2:1705009/93‑1 (MQ=255)
                                                                                  caacCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTa        <  1:2595487/93‑1 (MQ=255)
                                                                                   aacCGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTat       <  1:118156/93‑1 (MQ=255)
                                                                                     ccGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTatgc     <  1:5839287/93‑1 (MQ=255)
                                                                                      cGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTatgca    <  2:1221142/93‑1 (MQ=255)
                                                                                      cGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTatgca    <  2:1159117/93‑1 (MQ=255)
                                                                                      cGATTCACCCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTatgca    <  2:2659326/93‑1 (MQ=255)
                                                                                               |                                                                                         
ACACTGGTACTGCAAGAGGGGGCCTTCGCTCCTCTGACATTACTCAGTATTTTCACTACTACACCCCATCTGCCATCCATAACAACCGATTCACCCCCTCGCGAGGCGTTGTTGCCCGTCACCCGGACAGACACCGTTATGGACTTCTCTTGTTGCACGATCGATCAATTGTTATGCATG  >  NC_000913/4498453‑4498632

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