Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,966 T→A 5.6% intergenic (‑39/+18) 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,9660TA 5.6% 198.2 / 4.5 72intergenic (‑39/+18)ralR/recTRac prophage; endodeoxyribonuclease toxin RalR/Rac prophage; recombinase, DNA renaturation
Reads supporting (aligned to +/- strand):  ref base T (31/37);  new base A (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-02

TTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGC  >  NC_000913/1413826‑1414107
                                                                                                                                             |                                                                                                                                              
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                                                >  2:67966/1‑143 (MQ=255)
                                                                                              atatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc                                                <  1:702605/143‑1 (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                                   <  1:770782/143‑1 (MQ=255)
                                                                                                           tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCAtt                                   <  2:810291/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)
                                                                                                                                          gtaTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACgg    <  1:378500/141‑1 (MQ=255)
                                                                                                                                            aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGc  <  1:822850/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTCATCCATTGATACTGCACGCTGGATCTCAATTGATACGGGC  >  NC_000913/1413826‑1414107

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: