Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,471,443 A→G intergenic (‑177/‑181) yfdO ← / → b2359 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,471,4430AG85.8% 618.4 / 89.4 260intergenic (‑177/‑181)yfdO/b2359hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base G (116/107);  total (137/123)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

Only 100 of 260 total aligned reads displayed.

TGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTC  >  NC_000913/2471352‑2471533
                                                                                           |                                                                                          
tGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGGTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGAc                                                                                           >  2:1330150/1‑93 (MQ=255)
 gCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCg                                                                                          <  1:3504966/93‑1 (MQ=255)
    tgtCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATTTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGt                                                                                       <  1:5458927/93‑1 (MQ=255)
    tgtCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCgggg                                                                                       >  1:1713957/1‑92 (MQ=255)
        cGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGaaa                                                                                   >  2:6710284/1‑93 (MQ=255)
          gTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAAc                                                                                 >  2:7252479/1‑93 (MQ=255)
          gTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAAc                                                                                 >  2:1725897/1‑93 (MQ=255)
              cGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGa                                                                             <  1:2139441/93‑1 (MQ=255)
              cGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGa                                                                             <  1:7317658/93‑1 (MQ=255)
               gAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAt                                                                            <  1:3619648/93‑1 (MQ=255)
               gAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAt                                                                            >  1:3602875/1‑93 (MQ=255)
                aaGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAtt                                                                           <  2:5213745/93‑1 (MQ=255)
                aaGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAtt                                                                           <  1:3423594/93‑1 (MQ=255)
                 aGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAttt                                                                          >  1:4525660/1‑93 (MQ=255)
                 aGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGAttt                                                                          >  1:2833555/1‑93 (MQ=255)
                   cGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGAttttt                                                                        <  2:3137392/93‑1 (MQ=255)
                      ccGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGc                                                                     >  1:4313500/1‑93 (MQ=255)
                        ggACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCa                                                                   >  2:1218210/1‑93 (MQ=255)
                         gACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAg                                                                  >  2:215121/1‑93 (MQ=255)
                          aCGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGa                                                                 <  1:6197077/93‑1 (MQ=255)
                          aCGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGa                                                                 >  1:6379994/1‑93 (MQ=255)
                          aCGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGa                                                                 >  2:4212769/1‑93 (MQ=255)
                           cGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGAt                                                                >  1:6578074/1‑93 (MQ=255)
                           cGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGAt                                                                <  2:4925825/93‑1 (MQ=255)
                            gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGAtt                                                               >  2:4849060/1‑93 (MQ=255)
                            gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGAtt                                                               <  2:6600312/93‑1 (MQ=255)
                            gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGAtt                                                               >  1:5780954/1‑93 (MQ=255)
                            gagaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGGGAAAACCCGATTTTTTGCCCGAtt                                                               >  2:2855884/1‑93 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTc                                                             <  1:896408/93‑1 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTc                                                             >  1:4453240/1‑93 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTc                                                             >  1:7275140/1‑93 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGGGAAAACCCGATTTTTTGGCAGATTTc                                                             >  1:52393/1‑93 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTc                                                             <  2:2419752/93‑1 (MQ=255)
                              gaAGGTTTTTTATCTGACGGATCATGTCTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTc                                                             >  2:7049455/1‑93 (MQ=255)
                               aaGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCg                                                            >  1:6503229/1‑93 (MQ=255)
                                 ggTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGAc                                                          >  1:5294358/1‑93 (MQ=255)
                                  gTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACg                                                         <  2:3801751/93‑1 (MQ=255)
                                   ttttttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGc                                                        >  2:6033387/1‑93 (MQ=255)
                                   ttttttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGc                                                        >  2:6404014/1‑93 (MQ=255)
                                      tttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATc                                                     <  2:1374707/93‑1 (MQ=255)
                                      tttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATc                                                     <  1:694581/93‑1 (MQ=255)
                                      tttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATc                                                     <  1:4747269/93‑1 (MQ=255)
                                       ttATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGACTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCa                                                    >  2:6816003/1‑93 (MQ=255)
                                         aTCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCaaa                                                  >  2:439112/1‑93 (MQ=255)
                                         aTCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCaaa                                                  <  1:7123720/93‑1 (MQ=255)
                                         aTCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCaaa                                                  <  2:5757774/93‑1 (MQ=255)
                                         aTCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCaaa                                                  >  1:2243985/1‑93 (MQ=255)
                                           cTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAAtt                                                <  1:3124945/93‑1 (MQ=255)
                                               cGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGa                                            <  2:4596255/93‑1 (MQ=255)
                                               cGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGa                                            <  1:3806688/93‑1 (MQ=255)
                                               cGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGa                                            <  1:1884500/93‑1 (MQ=255)
                                                   tCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACggg                                        <  1:4258247/93‑1 (MQ=255)
                                                   tCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACggg                                        >  2:2886407/1‑93 (MQ=255)
                                                    cATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGt                                       >  1:3975806/1‑93 (MQ=255)
                                                     aTGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTc                                      >  2:3381487/1‑93 (MQ=255)
                                                     aTGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTc                                      <  2:861122/93‑1 (MQ=255)
                                                     aTGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTc                                      <  1:1768072/93‑1 (MQ=255)
                                                      tGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCa                                     <  1:3171649/93‑1 (MQ=255)
                                                      tGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCa                                     >  2:6410230/1‑93 (MQ=255)
                                                      tGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCa                                     >  1:5934334/1‑93 (MQ=255)
                                                      tGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCa                                     >  2:5849603/1‑93 (MQ=255)
                                                       gTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAg                                    >  2:3887812/1‑93 (MQ=255)
                                                        ttttGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGa                                   <  2:3219344/93‑1 (MQ=255)
                                                          ttGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGAtt                                 >  2:4810718/1‑93 (MQ=255)
                                                           tGATTTTACTGACGGATCCCCGCCAGATTCTGGCGGGGGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGAttt                                >  1:6427088/1‑93 (MQ=255)
                                                             aTTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTg                              <  1:7483335/93‑1 (MQ=255)
                                                             aTTTTACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTg                              >  2:2800299/1‑93 (MQ=255)
                                                              ttttACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGa                             <  2:6928897/93‑1 (MQ=255)
                                                               tttACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGAt                            >  2:2488334/1‑93 (MQ=255)
                                                               tttACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGAt                            <  2:7500259/93‑1 (MQ=255)
                                                                ttACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATg                           <  2:5227071/93‑1 (MQ=255)
                                                                 tACTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGc                          <  1:1579719/93‑1 (MQ=255)
                                                                 tACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGc                          >  2:184477/1‑93 (MQ=255)
                                                                   cTGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGt                        >  1:1135155/1‑93 (MQ=255)
                                                                   cTGACGGATCCCCGCCAGATTCTGGCGGGCGGAAACCCGATTTTTTGCCAGACTTCGAAGCATCAAATTTTTACGGGGCAGCTTTTGAtttgt                        >  1:146514/1‑89 (MQ=255)
                                                                    tGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTc                       <  1:6241873/93‑1 (MQ=255)
                                                                    tGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTc                       >  2:5368992/1‑93 (MQ=255)
                                                                    tGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTc                       <  1:4577/93‑1 (MQ=255)
                                                                    tGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTc                       <  2:3531458/93‑1 (MQ=255)
                                                                    tGACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTc                       >  1:6050178/1‑93 (MQ=255)
                                                                     gACGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCa                      >  1:5122503/1‑93 (MQ=255)
                                                                     aaCGGATCCCCGCCAGATTCTGGCGGGGGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGGCAGATTTTGATGCGTCa                      >  1:6688800/2‑93 (MQ=255)
                                                                       cGGATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGa                    <  2:555386/93‑1 (MQ=255)
                                                                       cGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGa                    <  1:1061929/93‑1 (MQ=255)
                                                                        ggATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGAt                   <  1:6876041/93‑1 (MQ=255)
                                                                         gATCCCCGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGAtt                  >  2:5801515/1‑93 (MQ=255)
                                                                             cccGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa              <  1:3689522/93‑1 (MQ=255)
                                                                             cccGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa              <  1:3258931/93‑1 (MQ=255)
                                                                             cccGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa              >  1:4908010/1‑93 (MQ=255)
                                                                             cccGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGa              <  2:7465800/93‑1 (MQ=255)
                                                                              ccGCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGAc             >  2:5004672/1‑93 (MQ=255)
                                                                               cGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACg            <  1:5619697/93‑1 (MQ=255)
                                                                                gCCAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACgg           >  2:386674/1‑93 (MQ=255)
                                                                                gCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACgg           <  1:49692/93‑1 (MQ=255)
                                                                                  cAGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGt         >  2:6184577/1‑93 (MQ=255)
                                                                                   aGATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTc        >  1:481454/1‑93 (MQ=255)
                                                                                    gATTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCa       >  2:4342765/1‑93 (MQ=255)
                                                                                     aTTCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACTGGTCag      <  2:6510725/93‑1 (MQ=255)
                                                                                      ttCTGGCGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCaga     <  2:6233274/93‑1 (MQ=255)
                                                                                         tGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTc  <  2:7181901/93‑1 (MQ=255)
                                                                                           |                                                                                          
TGCGTGTCCGGTTGCGAAGCGTCCGGACGAGAAGGTTTTTTATCTGACGGATCATGTTTTGATTTTACTGACGGATCCCCGCCAGATTCTGACGGGTGAAAACCCGATTTTTTGCCAGATTTCGACGCATCAAATTTTGACGGGTCAGATTTTGATGCGTCAGATTTTGACGGGTCAGAGTC  >  NC_000913/2471352‑2471533

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