Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,278 A→G intergenic (+87/‑28) talB → / → mog transaldolase B/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,2780AG89.8% 810.6 / 70.6 303intergenic (+87/‑28)talB/mogtransaldolase B/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/19);  new base G (112/160);  total (124/179)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 304 total aligned reads displayed.

CTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCA  >  NC_000913/9186‑9367
                                                                                            |                                                                                         
   tAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCga                                                                                        >  1:7010643/1‑93 (MQ=255)
   tAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACga                                                                                        <  2:7082601/93‑1 (MQ=255)
     aTCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCgaga                                                                                      <  2:3397383/93‑1 (MQ=255)
      tCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGAc                                                                                     >  1:5827602/1‑93 (MQ=255)
      tCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGAc                                                                                     >  2:5343327/1‑93 (MQ=255)
       cATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACt                                                                                    >  2:5200129/1‑93 (MQ=255)
       cATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACt                                                                                    >  1:5091115/1‑93 (MQ=255)
        aTTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTg                                                                                   <  1:1397901/93‑1 (MQ=255)
         ttCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGt                                                                                  <  2:1300615/93‑1 (MQ=255)
            ttAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTa                                                                               <  1:744479/93‑1 (MQ=255)
              aGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTaaa                                                                             <  1:681519/93‑1 (MQ=255)
              aGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTaaa                                                                             <  1:1002689/93‑1 (MQ=255)
               gCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAAc                                                                            <  2:2994009/93‑1 (MQ=255)
                cGTGACCGGGAAGTCGGCCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACg                                                                           <  1:6557497/93‑1 (MQ=255)
                 gTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACgg                                                                          >  2:4441625/1‑93 (MQ=255)
                 gTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACgg                                                                          <  2:4702580/93‑1 (MQ=255)
                   gACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGaa                                                                        >  1:4426497/1‑93 (MQ=255)
                    aCCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGaaa                                                                       <  1:5972504/93‑1 (MQ=255)
                      cGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGaaaaa                                                                     <  2:6216277/93‑1 (MQ=255)
                        ggAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATc                                                                   >  2:6459850/1‑93 (MQ=255)
                        ggAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATc                                                                   >  2:1640681/1‑93 (MQ=255)
                         gAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCt                                                                  >  1:3877298/1‑93 (MQ=255)
                         gAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCt                                                                  >  2:2750898/1‑93 (MQ=255)
                          aaGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCtt                                                                 >  1:541820/1‑93 (MQ=255)
                          aaGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCtt                                                                 <  1:4153370/93‑1 (MQ=255)
                           aGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTtg                                                                <  1:5308632/93‑1 (MQ=255)
                           aGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTtg                                                                >  1:4354533/1‑93 (MQ=255)
                            gTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTtga                                                               <  2:5037407/93‑1 (MQ=255)
                              cGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTtgatg                                                             <  1:814188/93‑1 (MQ=255)
                               ggTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTtgatga                                                            <  1:3677587/93‑1 (MQ=255)
                                gTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGaa                                                           >  2:4331342/1‑93 (MQ=255)
                                 tCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAAt                                                          <  1:4879207/93‑1 (MQ=255)
                                 tCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAAt                                                          <  1:1697836/93‑1 (MQ=255)
                                  cACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATa                                                         <  2:4214901/93‑1 (MQ=255)
                                   aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATAc                                                        >  2:2728806/1‑93 (MQ=255)
                                   aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATAc                                                        >  2:3491617/1‑93 (MQ=255)
                                   aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATAc                                                        <  2:6853703/93‑1 (MQ=255)
                                   aCGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATAc                                                        <  1:2291731/93‑1 (MQ=255)
                                    cGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACt                                                       <  2:1047597/93‑1 (MQ=255)
                                     gCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACtt                                                      <  2:3699601/93‑1 (MQ=255)
                                     gCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACtt                                                      <  1:1476343/93‑1 (MQ=255)
                                     gCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACtt                                                      <  1:2456990/93‑1 (MQ=255)
                                     gCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACtt                                                      <  1:29382/93‑1 (MQ=255)
                                       tACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTa                                                    <  2:2882854/93‑1 (MQ=255)
                                        aCCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTAc                                                   <  1:2669035/93‑1 (MQ=255)
                                           tcttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTa                                                <  2:4624166/93‑1 (MQ=255)
                                           tcttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTa                                                >  1:7260544/1‑93 (MQ=255)
                                             ttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTAtt                                              <  1:7109136/93‑1 (MQ=255)
                                             ttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTAtt                                              >  2:5541943/1‑93 (MQ=255)
                                             ttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTAtt                                              <  2:7253116/93‑1 (MQ=255)
                                             ttctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTAtt                                              <  2:1453513/93‑1 (MQ=255)
                                              tctGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTg                                             <  1:4089657/93‑1 (MQ=255)
                                                tGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGc                                           <  2:4923033/93‑1 (MQ=255)
                                                tGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGc                                           <  2:3554811/93‑1 (MQ=255)
                                                 gAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCt                                          <  1:1651429/93‑1 (MQ=255)
                                                    gCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAg                                       <  2:3613832/93‑1 (MQ=255)
                                                    gCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAg                                       <  1:1124718/93‑1 (MQ=255)
                                                      ctgtctgtCCCTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGtt                                     <  1:3238224/93‑1 (MQ=255)
                                                          ctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCa                                 <  1:922842/93‑1 (MQ=255)
                                                          ctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCa                                 <  1:3388523/93‑1 (MQ=255)
                                                          ctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCa                                 >  1:4754040/1‑93 (MQ=255)
                                                           tgtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCAt                                >  2:5968405/1‑93 (MQ=255)
                                                           tgtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCAt                                <  1:5165076/93‑1 (MQ=255)
                                                           tgtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCAt                                >  1:5737969/1‑93 (MQ=255)
                                                            gtcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCAtc                               <  1:4276055/93‑1 (MQ=255)
                                                             tcACTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCAtct                              >  2:853838/1‑93 (MQ=255)
                                                                cTCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTg                           <  2:2333682/93‑1 (MQ=255)
                                                                 tCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGa                          <  2:2516234/93‑1 (MQ=255)
                                                                 tCCCTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGa                          <  1:2997469/93‑1 (MQ=255)
                                                                  cccTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGAt                         <  1:6434909/93‑1 (MQ=255)
                                                                  cccTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGAt                         <  2:6543241/93‑1 (MQ=255)
                                                                   ccTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATc                        <  1:1945414/93‑1 (MQ=255)
                                                                    cTTCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATcg                       <  1:3914871/93‑1 (MQ=255)
                                                                     ttCGCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAAGGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATcgc                      >  1:735345/1‑93 (MQ=255)
                                                                     ttCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATcgc                      >  2:1836913/1‑93 (MQ=255)
                                                                        gCAGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCa                   >  2:2084265/1‑93 (MQ=255)
                                                                          aGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATc                 <  2:5625420/93‑1 (MQ=255)
                                                                          aGTGTATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATc                 <  1:6709807/93‑1 (MQ=255)
                                                                           gtgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATcc                >  2:613229/1‑93 (MQ=255)
                                                                           gtgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATcc                <  1:6564494/93‑1 (MQ=255)
                                                                           gtgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATcc                <  2:1940439/93‑1 (MQ=255)
                                                                           gtgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATcc                >  1:1843988/1‑93 (MQ=255)
                                                                           gtgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATcc                >  1:454065/1‑93 (MQ=255)
                                                                            tgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCa               >  2:2242209/1‑93 (MQ=255)
                                                                            tgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCa               <  2:4390203/93‑1 (MQ=255)
                                                                            tgtATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCa               <  2:3880470/93‑1 (MQ=255)
                                                                              tATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAgc             >  2:4809153/1‑93 (MQ=255)
                                                                              tATCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAgc             >  1:457958/1‑93 (MQ=255)
                                                                                tCATTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAgcgg           <  1:120497/93‑1 (MQ=255)
                                                                                tCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAgcgg           <  2:3076243/93‑1 (MQ=255)
                                                                                  aTTCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAgcggcg         >  2:4086869/1‑93 (MQ=255)
                                                                                   ttCTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGt        <  2:5696737/93‑1 (MQ=255)
                                                                                     cTGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGttt      <  2:7314179/93‑1 (MQ=255)
                                                                                      tGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTa     >  1:4999811/1‑93 (MQ=255)
                                                                                      tGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTa     <  2:1739783/93‑1 (MQ=255)
                                                                                      tGTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTa     >  1:1512281/1‑93 (MQ=255)
                                                                                       gTTTAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTAt    >  1:2101323/1‑93 (MQ=255)
                                                                                        tttAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATc   <  1:629469/93‑1 (MQ=255)
                                                                                        tttAGCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATc   >  2:7485229/1‑93 (MQ=255)
                                                                                         ttAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCa  <  2:3181153/93‑1 (MQ=255)
                                                                                            |                                                                                         
CTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCA  >  NC_000913/9186‑9367

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