Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,964 C→T 7.2% intergenic (‑37/+20) ralR ← / ← recT Rac prophage; endodeoxyribonuclease toxin RalR/Rac prophage; recombinase, DNA renaturation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,9640CT 7.2% 180.3 / 5.3 70intergenic (‑37/+20)ralR/recTRac prophage; endodeoxyribonuclease toxin RalR/Rac prophage; recombinase, DNA renaturation
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-02

TTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATA  >  NC_000913/1413826‑1414102
                                                                                                                                           |                                                                                                                                           
ttCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                                                         <  1:812544/143‑1 (MQ=255)
ttCGCATCCGTTACACTGCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACTACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                                                         <  1:558505/143‑1 (MQ=255)
   gCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCCTAAAAAATATTTCAAGTTGACAGAGCATTacac                                                                                                                                      >  1:680824/1‑143 (MQ=255)
        cGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAATATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCa                                                                                                                                 <  1:697404/143‑1 (MQ=255)
        cGTTACACTTCGCTCGGTAATATCCTGAAATGACTTTCACCGTTACTGATGGACAACCACAAAAGGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAACACTATTTCAAGTTCGCGGTGCATTACACCGccg                                                                                                                                 >  1:1019089/1‑142 (MQ=255)
         gTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaa                                                                                                                                <  1:887979/143‑2 (MQ=255)
          ttACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAgg                                                                                                                               >  2:278564/1‑143 (MQ=255)
               cTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                                                          <  1:1006604/143‑1 (MQ=255)
                    cTCGGCAATATTCTGAAATGGCTTTCACCGTGACGTACAGACCAATTACATAAGCGCATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                                                     <  2:1016114/143‑1 (MQ=255)
                     tCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatta                                                                                                                    >  2:770782/1‑142 (MQ=255)
                       ggTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcgt                                                                                                                  >  1:76961/1‑141 (MQ=255)
                          aaTTTCCTCAAATGCCTGTCACCGTTCCTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAGTTTTTCTCATAAAAAATATGTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                                                               <  2:628299/143‑1 (MQ=255)
                            tatCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGaa                                                                                                             <  1:119960/143‑1 (MQ=255)
                               ccTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTa                                                                                                          <  2:164018/143‑1 (MQ=255)
                                cTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATATCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCATCTGAATTAt                                                                                                         >  2:113154/1‑143 (MQ=255)
                                 tGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                                                        <  2:575271/143‑1 (MQ=255)
                                  gAAATGGCTTTCACCGTTAATGATGGACAACCACAAAATGGACATGGTTTAACATTGTAATAACTCATAATTTATCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTATGAATTaaag                                                                                                       >  2:552113/1‑140 (MQ=255)
                                   aaaTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATCAGACATATATCAAGTTGGCAGTACATTACACCGACAGGCTGAATAATTCCTCTGAATTATCGa                                                                                                      >  1:630412/1‑143 (MQ=255)
                                   aaaTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGa                                                                                                      >  1:614655/1‑143 (MQ=255)
                                    gaTGGCTGTCACCGGTACTGATGGACAACCACAAAATGGACATGGGTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAt                                                                                                     <  2:1017150/142‑1 (MQ=255)
                                    caTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAt                                                                                                     <  2:454445/142‑1 (MQ=255)
                                    aagtgCTTTCACCCTTACTGATGGACAACCACAAAATGGACATGGTGTAACATTGTCATATCTCCTAATTTTTCTCATCAAAAATATTTCACGTTGGCGCTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAt                                                                                                     <  1:485548/139‑1 (MQ=255)
                                    aaTGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTATACATTAttggg                                                                                                     <  2:452006/143‑5 (MQ=255)
                                          tttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTGCTCTGAATTATCGATTACACt                                                                                               >  1:519526/1‑143 (MQ=255)
                                                  ttACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTTAATTATCGTTTACAATGtcgcagc                                                                                       >  1:683676/1‑137 (MQ=255)
                                                  ttACTGATCGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTccc                                                                                       <  1:83842/143‑1 (MQ=255)
                                                      tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                                                                   >  2:71107/1‑143 (MQ=255)
                                                      tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATAGAGTACACGGTATTCCCCggg                                                                                   >  1:197999/1‑142 (MQ=255)
                                                         tGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGGTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTaa                                                                                >  1:338810/1‑143 (MQ=255)
                                                           gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATa                                                                              <  2:217858/143‑1 (MQ=255)
                                                                aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATCAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                                                          <  2:568927/143‑1 (MQ=255)
                                                                aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                                                          <  1:576775/143‑1 (MQ=255)
                                                                 ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                                                         <  2:1052493/143‑1 (MQ=255)
                                                                 ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCACATAAAAAATATTTCAAATTCGCGTTGCATTACACCGCCAGCCTGAATTATTCCTCGGAATTATCGATTACACCATATTCCTCTGTTGATCCagaa                                                                         >  2:213767/1‑142 (MQ=255)
                                                                   acaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGa                                                                       >  2:822850/1‑143 (MQ=255)
                                                                       aaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCCGTTAATACAGAGGAATCt                                                                   >  2:680863/1‑143 (MQ=255)
                                                                         tGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATCCAGAGGAATCTGc                                                                 >  2:396384/1‑143 (MQ=255)
                                                                         tGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                                                                 <  2:783094/143‑1 (MQ=255)
                                                                         tGGACATGGTTTAACATTGTCAGATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                                                                 <  1:452598/143‑1 (MQ=255)
                                                                               tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                                                          <  1:235732/143‑1 (MQ=255)
                                                                                ggTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGCATCTTCAGGATct                                                         >  2:787733/1‑142 (MQ=255)
                                                                                  tttAACATTGTCATATCTCATAATTTTCCTCATAAAAAATATTTCAAGTTGGCCGTGCGTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGAt                                                       <  2:662935/143‑1 (MQ=255)
                                                                                      aCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTATACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc                                                   >  2:551260/1‑143 (MQ=255)
                                                                                      aCATCGTCACACCTCACAAGTTTTCGCATAAAAAATATTTCAAGTTGGCGGTGCATTTCACCGCCAGGCTAAATTATGCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc                                                   <  1:300893/143‑1 (MQ=255)
                                                                                         ttGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCTGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt                                                <  1:713822/143‑1 (MQ=255)
                                                                                              atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTCTTCCCCGGTTCATCCAGCGGAACCTGCAGGATCGATTGTCCGTGGGTc                                           >  2:299641/1‑143 (MQ=255)
                                                                                              atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc                                           <  1:702605/143‑1 (MQ=255)
                                                                                              atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc                                           >  2:67966/1‑143 (MQ=255)
                                                                                               tatCTCATAATTTTTCTCATAAAAACTATCCATTTTTGTAGTTTCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTAAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTcc                                          >  2:342493/1‑143 (MQ=255)
                                                                                                 tctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGCATTTTTAACCTTAATTATCGACTACACTGTAGTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt                                        <  1:514524/143‑1 (MQ=255)
                                                                                                   tcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAAAGAGGAATCTGCAGGATCGATTGGCAGTGGTTCCtttt                                      >  1:384031/1‑143 (MQ=255)
                                                                                                     aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCa                                    <  1:278564/143‑1 (MQ=255)
                                                                                                     aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTAAACGGCCAGGATGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGGTAATACAGAGGAATCTGCAGGATCGATTGTCAGGGGGTGCTTTTCa                                    >  1:994610/1‑143 (MQ=255)
                                                                                                       aaTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGCATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGACAGAGGTTACTTTTCATc                                  >  1:53650/1‑143 (MQ=255)
                                                                                                           tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCAtt                              <  2:810291/143‑1 (MQ=255)
                                                                                                           tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCAtt                              <  1:770782/143‑1 (MQ=255)
                                                                                                               tcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctc                                                                                                                 <  2:263175/57‑1 (MQ=255)
                                                                                                               tcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctc                                                                                                                 >  1:263175/1‑57 (MQ=255)
                                                                                                               tcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCCGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATa                          >  2:112330/1‑143 (MQ=255)
                                                                                                                 ataataTATATTTCAAGTTTCCTGTGCATTACACCGCCATGCTGACTTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACt                        <  2:33160/143‑1 (MQ=255)
                                                                                                                 aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACt                        <  1:884036/143‑1 (MQ=255)
                                                                                                                    aaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGGGGATCCTTTTCATACATTGGTACTGCa                     >  1:443004/1‑143 (MQ=255)
                                                                                                                           tttCAAGTTGGCGGTGTATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGa              >  2:378500/1‑143 (MQ=255)
                                                                                                                            ttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGCTTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAt             <  2:972653/143‑1 (MQ=255)
                                                                                                                            ttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGAt             <  1:136297/143‑1 (MQ=255)
                                                                                                                                 gTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTACACGGCGAACAAAGAGGAATATGCAGAGATGTCTTTCGGTGGTTCCTTTTCATCCATTGAAACTGCCCGCTGGATCTCaa        >  2:784964/1‑143 (MQ=255)
                                                                                                                                  ctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAAt       >  1:837509/2‑143 (MQ=255)
                                                                                                                                   tgtgggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAAtt      <  2:974518/139‑1 (MQ=255)
                                                                                                                                   tGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCACtt      >  1:488992/1‑143 (MQ=255)
                                                                                                                                       gttGCATTACTCCTCGAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATa  <  2:479258/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATA  >  NC_000913/1413826‑1414102

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