Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,979 +G coding (1511/1833 nt) recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,9791.G99.8% 1846.5 / ‑4.4 532coding (1511/1833 nt)recQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (284/247);  total (284/248)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

Only 100 of 589 total aligned reads displayed.

CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060
                                                                                             |                                                                                     
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      <  2:5458794/93‑3 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      >  1:1827769/1‑91 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      <  1:4318728/93‑3 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      >  1:3431206/1‑91 (MQ=255)
  tGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgc                                                                                      <  1:4526114/93‑3 (MQ=255)
   gATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcg                                                                                     <  2:3021543/93‑4 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    <  1:1369824/93‑5 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    <  2:1383600/93‑5 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    <  2:7375725/93‑5 (MQ=255)
    aTCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCtgcgc                                                                                    >  2:5718492/1‑89 (MQ=255)
         ccAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCg‑‑gcga                                                                             >  2:1180064/1‑90 (MQ=255)
                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     >  2:8289535/1‑93 (MQ=255)
                 ttCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCtt                                                                     >  2:2534418/1‑93 (MQ=255)
                    cccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCG‑‑GCTAATCCTCTTTGc                                                                  <  1:1653531/92‑1 (MQ=255)
                     ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCTCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCa                                                                 <  1:7060824/93‑1 (MQ=255)
                     ccTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCa                                                                 <  2:7949319/93‑1 (MQ=255)
                       tCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAAc                                                               <  1:5603052/93‑1 (MQ=255)
                        cGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACt                                                              <  2:4491023/93‑1 (MQ=255)
                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             >  2:1104208/1‑93 (MQ=255)
                         ggCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACtt                                                             >  2:5158956/1‑93 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            <  1:729670/93‑1 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            >  1:1560574/1‑93 (MQ=255)
                          gCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtg                                                            <  1:1443977/93‑1 (MQ=255)
                           ccTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgc                                                           <  1:3021587/93‑1 (MQ=255)
                           ccTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgc                                                           <  1:2234252/93‑1 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:3334941/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:7496727/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  1:2570745/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  2:85079/1‑93 (MQ=255)
                            cTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgcc                                                          >  2:1529053/1‑93 (MQ=255)
                              ggTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGAAACTtgccgt                                                        >  2:4087323/1‑93 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtg                                                       <  2:7457090/93‑1 (MQ=255)
                               gTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtg                                                       <  2:1755320/93‑1 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     >  2:3639594/1‑93 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     <  2:2699824/93‑1 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     >  2:2372079/1‑93 (MQ=255)
                                 gACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgcc                                                     >  1:2659269/1‑93 (MQ=255)
                                  aCGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTtgccgtgccg                                                    >  1:5678389/1‑93 (MQ=255)
                                   cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCggc                                                   >  1:8696630/1‑93 (MQ=255)
                                   cGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCcgc                                                   <  1:8600061/93‑1 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                >  1:5442823/1‑93 (MQ=255)
                                      aaaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTa                                                <  2:7163395/93‑1 (MQ=255)
                                       aaaTATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTAt                                               <  2:1676830/93‑1 (MQ=255)
                                        aaTATTGCTCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATc                                              <  1:4958194/93‑1 (MQ=255)
                                         atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             >  1:5855180/1‑93 (MQ=255)
                                         atatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCg                                             <  2:1899351/93‑1 (MQ=255)
                                          tatTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGt                                            >  1:3131750/1‑93 (MQ=255)
                                            ttGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgg                                          >  1:6219767/1‑93 (MQ=255)
                                             tGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgc                                         <  2:6614977/93‑1 (MQ=255)
                                              gCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTgggg                                        >  2:4244141/1‑91 (MQ=255)
                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  2:2190267/1‑93 (MQ=255)
                                               cccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGgcgc                                       >  1:6469654/1‑93 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  2:4860036/1‑93 (MQ=255)
                                                ccAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCt                                      >  1:4053286/1‑93 (MQ=255)
                                                 cAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTc                                     >  1:7682385/1‑93 (MQ=255)
                                                  aGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCa                                    >  2:6791277/1‑93 (MQ=255)
                                                   gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTTCCGTGCCGCGTATCGTGGGGCTCaa                                   >  2:8595274/1‑93 (MQ=255)
                                                   gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                                   >  2:6067850/1‑93 (MQ=255)
                                                   gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCaa                                   >  1:1203445/1‑93 (MQ=255)
                                                   gCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGACGCGTATCGTGGCGCTCaa                                   >  1:5212241/1‑93 (MQ=255)
                                                     aTTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAc                                 >  2:313185/1‑93 (MQ=255)
                                                      ttCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAAcc                                >  2:6220063/1‑93 (MQ=255)
                                                          gCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGaaa                            <  1:1961979/93‑1 (MQ=255)
                                                           cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:6392245/93‑1 (MQ=255)
                                                           cccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  1:785851/93‑1 (MQ=255)
                                                           accTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACGTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAg                           <  2:5946847/92‑1 (MQ=255)
                                                            ccTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGc                          >  1:4150884/1‑93 (MQ=255)
                                                             cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         >  2:8774501/1‑93 (MQ=255)
                                                             cTACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCg                         <  2:2862315/93‑1 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        <  1:2428346/93‑1 (MQ=255)
                                                              tACAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGa                        >  1:4849878/1‑93 (MQ=255)
                                                               aCAACTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGAt                       >  1:1061815/1‑93 (MQ=255)
                                                                 aaCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCGCAAACCGGAAGAGATGc                     >  2:3768114/1‑93 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    <  1:3949056/93‑1 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:5047780/1‑93 (MQ=255)
                                                                  aCTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCa                    >  1:957223/1‑93 (MQ=255)
                                                                   cTGACAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAg                   >  2:7619875/1‑93 (MQ=255)
                                                                      aCAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGaaa                <  2:7841230/93‑1 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               >  2:2468632/1‑93 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               >  2:2117000/1‑93 (MQ=255)
                                                                       cAGAGGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAAt               <  1:8266905/93‑1 (MQ=255)
                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  1:6038055/93‑1 (MQ=255)
                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              <  2:5329117/93‑1 (MQ=255)
                                                                        agagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATc              >  2:774023/1‑93 (MQ=255)
                                                                         gagGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCg             >  1:6827979/1‑93 (MQ=255)
                                                                          agGCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGt            <  2:7423489/93‑1 (MQ=255)
                                                                           ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:3542905/93‑1 (MQ=255)
                                                                           ggCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGtt           <  2:7029097/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  2:3908612/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:8433529/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:6908145/93‑1 (MQ=255)
                                                                            gCCGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTc          <  1:6656785/93‑1 (MQ=255)
                                                                             ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         <  2:5540436/93‑1 (MQ=255)
                                                                             ccGCGCGCCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcg         >  1:5348950/1‑93 (MQ=255)
                                                                               gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:4985054/1‑93 (MQ=255)
                                                                               gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggc       >  2:3300374/1‑93 (MQ=255)
                                                                               gcgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAACCGATGCAGAAATCGTTcggg       >  1:2018766/1‑92 (MQ=255)
                                                                                cgcgcCCGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTcggcg      <  1:8759819/93‑1 (MQ=255)
                                                                                    cccGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  <  1:369899/93‑1 (MQ=255)
                                                                                    cccGGTGCTGCGCGGCG‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCaa  >  2:5821927/1‑93 (MQ=255)
                                                                                             |                                                                                     
CGTGATCCGCCAGCTGATTCACCTCGGCCTGGTGACGCAAAATATTGCCCAGCATTCTGCCCTACAACTGACAGAGGCCGCGCGCCCGGTGCTCGCG‑‑G‑‑AATCCTCTTTGCAACTTGCCGTGCCGCGTATCGTGGCGCTCAAACCGAAAGCGATGCAGAAATCGTTCGGCGGCAA  >  NC_000913/4004887‑4005060

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