Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,066,300 T→C intergenic (‑957/‑330) trs56 ← / → b1995 IS5 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,066,3000TC99.7% 1196.0 / ‑7.2 362intergenic (‑957/‑330)trs56/b1995IS5 transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (129/232);  total (130/232)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 365 total aligned reads displayed.

TCTGGCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGAC  >  NC_000913/2066210‑2066389
                                                                                          |                                                                                         
tctGGCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCt                                                                                         <  2:3131719/93‑1 (MQ=255)
 ctGGCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCtt                                                                                        >  2:1179147/1‑93 (MQ=255)
   ggCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAc                                                                                      <  1:727689/93‑1 (MQ=255)
   ggCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAc                                                                                      >  1:4927378/1‑93 (MQ=255)
   ggCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAc                                                                                      <  1:1853158/93‑1 (MQ=255)
   ggCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAc                                                                                      <  1:3872205/93‑1 (MQ=255)
    gCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAcc                                                                                     >  2:9116384/1‑93 (MQ=255)
    gCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTAcc                                                                                     >  2:2879220/1‑93 (MQ=255)
      cATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGc                                                                                   <  2:3595042/93‑1 (MQ=255)
       aTTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCa                                                                                  <  2:4146446/93‑1 (MQ=255)
        ttAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCat                                                                                 >  1:7017533/1‑93 (MQ=255)
          agagTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCatat                                                                               <  2:4904237/93‑1 (MQ=255)
           gagTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCAtata                                                                              <  2:6868432/93‑1 (MQ=255)
           gagTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCAtata                                                                              <  1:8799796/93‑1 (MQ=255)
           gagTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCAtata                                                                              >  1:4957310/1‑93 (MQ=255)
           gagTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCAtata                                                                              >  1:9350151/1‑93 (MQ=255)
              tAACCTTCTGCATTCATCCTTTTGTAAAAAGGTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACgg                                                                           >  1:5505003/1‑93 (MQ=255)
               aaCCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGt                                                                          <  2:9896432/93‑1 (MQ=255)
                aCCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTc                                                                         >  2:2702409/1‑93 (MQ=255)
                 ccTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGGCCTGATACTGATATGAGCCTTACCGCATATACGGTCa                                                                        <  1:8306247/93‑1 (MQ=255)
                   ttCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAAt                                                                      >  2:5285241/1‑93 (MQ=255)
                   ttCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAAt                                                                      >  2:6388823/1‑93 (MQ=255)
                    tCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAAtt                                                                     <  2:7190306/93‑1 (MQ=255)
                    tCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAAtt                                                                     >  1:6106227/1‑93 (MQ=255)
                    tCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAAtt                                                                     <  1:7646194/93‑1 (MQ=255)
                      tGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTc                                                                   >  2:9595934/1‑93 (MQ=255)
                      tGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTc                                                                   <  2:6501391/93‑1 (MQ=255)
                       gCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCa                                                                  <  2:2267625/93‑1 (MQ=255)
                        cATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAg                                                                 >  2:9534411/1‑93 (MQ=255)
                         aTTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGc                                                                >  1:6637233/1‑93 (MQ=255)
                         aTTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGc                                                                >  1:1699123/1‑93 (MQ=255)
                         aTTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGc                                                                <  1:4526260/93‑1 (MQ=255)
                          ttCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCa                                                               >  1:1113697/1‑93 (MQ=255)
                          ttCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCa                                                               <  1:4344212/93‑1 (MQ=255)
                            cATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaat                                                             <  2:2514756/93‑1 (MQ=255)
                            cATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaat                                                             <  2:4389601/93‑1 (MQ=255)
                             aTCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaatt                                                            <  1:5654222/93‑1 (MQ=255)
                              tCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaatta                                                           <  2:7780910/93‑1 (MQ=255)
                                 ttttGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaattaatt                                                        <  2:467987/93‑1 (MQ=255)
                                 ttttGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCaattaatt                                                        >  2:3372528/1‑93 (MQ=255)
                                  tttGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAattaatta                                                       <  1:148139/93‑1 (MQ=255)
                                  tttGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAattaatta                                                       >  2:2962522/1‑93 (MQ=255)
                                    tGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACa                                                     <  2:511971/93‑1 (MQ=255)
                                    tGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACa                                                     >  2:1680002/1‑93 (MQ=255)
                                    tGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACa                                                     <  1:7358966/93‑1 (MQ=255)
                                     gTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACAt                                                    <  2:34798/93‑1 (MQ=255)
                                      tAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACAtt                                                   <  1:3425549/93‑1 (MQ=255)
                                       aaaaaGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTa                                                  >  2:9413341/1‑93 (MQ=255)
                                        aaaaGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTAt                                                 <  2:3949738/93‑1 (MQ=255)
                                           aGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCa                                              <  2:9474902/93‑1 (MQ=255)
                                           aGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCa                                              <  1:7994057/93‑1 (MQ=255)
                                            gTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCAc                                             <  1:6323135/93‑1 (MQ=255)
                                              ttATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGc                                           >  2:2150879/1‑93 (MQ=255)
                                              ttATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGc                                           <  1:3721619/93‑1 (MQ=255)
                                              ttATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGc                                           <  2:6246860/93‑1 (MQ=255)
                                              ttATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGc                                           <  2:7844222/93‑1 (MQ=255)
                                              ttATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGc                                           <  2:811781/93‑1 (MQ=255)
                                               tataTTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGcc                                          >  2:5781843/1‑93 (MQ=255)
                                                 tatTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCaa                                        >  2:7827737/1‑93 (MQ=255)
                                                 tatTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCaa                                        <  2:8500388/93‑1 (MQ=255)
                                                 tatTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCaa                                        >  2:4088367/1‑93 (MQ=255)
                                                  atTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCaaa                                       >  2:6609573/1‑93 (MQ=255)
                                                   ttAGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAg                                      <  1:2526548/93‑1 (MQ=255)
                                                     aGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTa                                    <  1:8081989/93‑1 (MQ=255)
                                                     aGTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTa                                    <  1:3323320/93‑1 (MQ=255)
                                                      gTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAt                                   <  1:8247045/93‑1 (MQ=255)
                                                      gTATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAt                                   <  1:10008955/93‑1 (MQ=255)
                                                       tATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAtt                                  <  2:3286912/93‑1 (MQ=255)
                                                       tATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAtt                                  <  2:7861608/93‑1 (MQ=255)
                                                       tATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAtt                                  <  2:2385661/93‑1 (MQ=255)
                                                       tATCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAtt                                  <  1:2288927/93‑1 (MQ=255)
                                                        aTCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTAttt                                 <  1:5902013/93‑1 (MQ=255)
                                                         tCAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTg                                <  2:9018630/93‑1 (MQ=255)
                                                          cAGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGt                               <  1:8598514/93‑1 (MQ=255)
                                                           aGCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGtc                              >  2:10027673/1‑93 (MQ=255)
                                                            gCAATTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGtca                             <  1:2680792/93‑1 (MQ=255)
                                                              aaTTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGtcatc                           <  2:3853262/93‑1 (MQ=255)
                                                               aTTAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGtcatca                          <  2:5557351/93‑1 (MQ=255)
                                                                ttAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATcac                         <  1:6200803/93‑1 (MQ=255)
                                                                 tAACCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATcaca                        <  2:4005097/93‑1 (MQ=255)
                                                                  aaCCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACaa                       <  2:364510/93‑1 (MQ=255)
                                                                  aaCCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACaa                       <  2:6145404/93‑1 (MQ=255)
                                                                   aCCGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAat                      <  2:8539677/93‑1 (MQ=255)
                                                                    ccGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAatg                     <  1:6814480/93‑1 (MQ=255)
                                                                    ccGGACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAatg                     >  1:7425908/1‑93 (MQ=255)
                                                                     cGGACCTGATACTGATATGATCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAatga                    <  1:1551237/93‑1 (MQ=255)
                                                                      ggACCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAatgat                   <  2:5036534/93‑1 (MQ=255)
                                                                        aCCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATgg                 >  2:7883289/1‑93 (MQ=255)
                                                                        aCCTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATgg                 <  1:5210535/93‑1 (MQ=255)
                                                                         ccTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGt                <  1:8999780/93‑1 (MQ=255)
                                                                         ccTGATACTGATATGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGt                <  1:8480426/93‑1 (MQ=255)
                                                                           tgatactgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTAc              <  2:6771413/93‑1 (MQ=255)
                                                                            gatactgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTAcc             <  2:7581641/93‑1 (MQ=255)
                                                                            gatactgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTAcc             <  1:3192660/93‑1 (MQ=255)
                                                                             atactgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACct            >  1:9623626/1‑93 (MQ=255)
                                                                               actgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACcttc          <  2:1938668/93‑1 (MQ=255)
                                                                                 tgataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACcttctt        <  2:7423754/93‑1 (MQ=255)
                                                                                  gataTGAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCttt       <  1:4078867/93‑1 (MQ=255)
                                                                                       gAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGac  <  2:1086282/93‑1 (MQ=255)
                                                                                       gAGCCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGac  <  1:2011534/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCTGGCCATTAGAGTAACCTTCTGCATTCATCCTTTTGTAAAAAGTTTATATTAGTATCAGCAATTAACCGGACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGAC  >  NC_000913/2066210‑2066389

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