Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 740,004 G→A M166I (ATG→ATA)  phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913740,0040GA99.1% 371.1 / ‑6.1 108M166I (ATG→ATA) phrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (52/55);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATA  >  NC_000913/739909‑740098
                                                                                               |                                                                                              
cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCgg                                                                                            <  1:2132703/100‑1 (MQ=255)
 tGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGa                                                                                           >  2:1505529/1‑100 (MQ=255)
  gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGa                                                                                           >  2:2361899/1‑99 (MQ=255)
  gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAg                                                                                          <  2:2380423/100‑1 (MQ=255)
  gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAg                                                                                          <  2:2317917/100‑1 (MQ=255)
    cgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTg                                                                                        <  2:2273937/100‑1 (MQ=255)
    cgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTg                                                                                        <  1:1199141/100‑1 (MQ=255)
        tGGCGCGGTGATGACCGGTAATCACGAGATGTGCAAAGTTTTAACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTc                                                                                    <  2:2143539/100‑1 (MQ=255)
           cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCt                                                                                 <  1:2108865/100‑1 (MQ=255)
           cgcgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCt                                                                                 <  2:2324/100‑1 (MQ=255)
               gTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTgcgc                                                                             >  1:108482/1‑100 (MQ=255)
                 gatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCa                                                                           >  1:93944/1‑100 (MQ=255)
                  atgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaa                                                                          <  1:1095235/100‑1 (MQ=255)
                   tgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCgaa                                                                         >  2:2032285/1‑100 (MQ=255)
                     aCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAg                                                                       <  2:373095/100‑1 (MQ=255)
                       cGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGtt                                                                     >  1:2757010/1‑100 (MQ=255)
                        ggTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                    >  2:2784649/1‑100 (MQ=255)
                         gTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCg                                                                   <  2:1503403/100‑1 (MQ=255)
                         gTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCg                                                                   <  2:1024881/100‑1 (MQ=255)
                          tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGACACCGGAGTGCGTCGCTGTGCCAAAAGTTCgt                                                                  <  1:2054661/100‑1 (MQ=255)
                          tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCGGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                  <  1:2529899/100‑1 (MQ=255)
                           aaTCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTACGCCAAAAGTTCgta                                                                 >  1:3011937/1‑100 (MQ=255)
                             tCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagt                                                               <  1:554804/100‑1 (MQ=255)
                              cACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCTCCAAAAGTTCgtagta                                                              <  2:1091251/100‑1 (MQ=255)
                               aCGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGtagtag                                                             <  2:1015553/100‑1 (MQ=255)
                                cGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGc                                                            <  1:679872/100‑1 (MQ=255)
                                 gagaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCg                                                           <  2:2529704/100‑1 (MQ=255)
                                   gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                         >  1:1938314/1‑100 (MQ=255)
                                   gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                         >  2:741290/1‑100 (MQ=255)
                                    aTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGAt                                                        >  2:3022080/1‑100 (MQ=255)
                                     tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAATGGATACAGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                         >  2:629585/1‑98 (MQ=255)
                                      gTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCg                                                      >  1:2528187/1‑100 (MQ=255)
                                      gTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCg                                                      >  2:1925239/1‑100 (MQ=255)
                                       tACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                                     >  2:1213799/1‑100 (MQ=255)
                                        aCAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAt                                                    <  1:2236628/100‑1 (MQ=255)
                                         cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                   >  1:476224/1‑100 (MQ=255)
                                         cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                   <  1:1343890/100‑1 (MQ=255)
                                         cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                   <  1:255054/100‑1 (MQ=255)
                                          aaaGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATag                                                  <  2:12943/100‑1 (MQ=255)
                                          aaaGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATag                                                  <  1:2017242/100‑1 (MQ=255)
                                           aaGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATaga                                                 >  2:1650254/1‑100 (MQ=255)
                                            aGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATagag                                                <  1:2317373/100‑1 (MQ=255)
                                             gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                               >  1:1876119/1‑100 (MQ=255)
                                             gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                               <  1:519566/100‑1 (MQ=255)
                                               cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                             >  2:166298/1‑100 (MQ=255)
                                                  tACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCg                                          <  1:1259416/100‑1 (MQ=255)
                                                  tACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCg                                          <  2:1256714/100‑1 (MQ=255)
                                                  tACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCg                                          <  1:1925262/100‑1 (MQ=255)
                                                   aCGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAGAGAGCCCTCGc                                         >  2:449015/1‑100 (MQ=255)
                                                     gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                       >  1:2708000/1‑100 (MQ=255)
                                                     gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                       <  1:1063568/100‑1 (MQ=255)
                                                      ccTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccat                                      >  2:1182686/1‑100 (MQ=255)
                                                      ccTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccat                                      >  2:2686140/1‑100 (MQ=255)
                                                       cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                     <  1:2855229/100‑1 (MQ=255)
                                                       cTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatc                                     <  2:1555788/100‑1 (MQ=255)
                                                        ttttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                    >  1:113545/1‑100 (MQ=255)
                                                         tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                   <  1:1280141/100‑1 (MQ=255)
                                                         tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                   <  2:2709471/100‑1 (MQ=255)
                                                         tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                   <  1:2792012/100‑1 (MQ=255)
                                                          ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                  <  1:848705/100‑1 (MQ=255)
                                                          ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCACTCGccatccat                                  >  1:2422611/1‑100 (MQ=255)
                                                           tAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCAtt                                 <  1:2784688/100‑1 (MQ=255)
                                                             aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                               <  2:1207613/100‑1 (MQ=255)
                                                             aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                               <  1:2748158/100‑1 (MQ=255)
                                                              gAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACg                              >  1:1599846/1‑100 (MQ=255)
                                                              gAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACg                              >  2:2154568/1‑100 (MQ=255)
                                                               aaTGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGc                             <  2:1327149/100‑1 (MQ=255)
                                                               aaTGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGc                             <  2:61480/100‑1 (MQ=255)
                                                                  gCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTgn                          <  2:2972374/100‑2 (MQ=255)
                                                                  gCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGa                          >  1:1249163/1‑100 (MQ=255)
                                                                  gCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGa                          >  1:2118420/1‑100 (MQ=255)
                                                                   ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGaa                         <  2:777633/100‑1 (MQ=255)
                                                                   ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGaa                         >  2:1105335/1‑100 (MQ=255)
                                                                   ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGaa                         >  1:2838634/1‑100 (MQ=255)
                                                                    cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                        <  1:80873/100‑1 (MQ=255)
                                                                      ggCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTa                      >  1:1905457/1‑100 (MQ=255)
                                                                       gCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTAt                     >  2:44722/1‑100 (MQ=255)
                                                                        cTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                    >  1:74640/1‑100 (MQ=255)
                                                                         tGAAACGGCTGCGGGAAGGGATACCGGTGTGCGTCGCTGCTCCAATAGTTCGTAGTAGCGGATCGATAGAGCCCCCGCCATCTATTACGCTGAATTATcc                   <  2:256948/100‑1 (MQ=255)
                                                                         tGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGGATTAt                     >  2:2899531/1‑98 (MQ=255)
                                                                         tGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc                   >  1:1290466/1‑100 (MQ=255)
                                                                         tGAAACGCCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATcc                   >  1:1608894/1‑100 (MQ=255)
                                                                          gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                  >  1:694948/1‑100 (MQ=255)
                                                                          gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                  >  1:898475/1‑100 (MQ=255)
                                                                          gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                  >  1:1770104/1‑100 (MQ=255)
                                                                          gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                  <  2:1422767/100‑1 (MQ=255)
                                                                              cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc              <  1:1182698/100‑1 (MQ=255)
                                                                              cGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTc              >  2:1943344/1‑100 (MQ=255)
                                                                               ggCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCa             >  2:2240778/1‑100 (MQ=255)
                                                                                 cTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGt           <  1:166305/100‑1 (MQ=255)
                                                                                  tGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc          <  1:1241225/100‑1 (MQ=255)
                                                                                  tGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc          >  2:2819150/1‑100 (MQ=255)
                                                                                   gCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt         <  1:1844547/100‑1 (MQ=255)
                                                                                     gggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCttt       >  2:1917671/1‑100 (MQ=255)
                                                                                      ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc      <  1:2324475/100‑1 (MQ=255)
                                                                                      ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc      <  1:3022121/100‑1 (MQ=255)
                                                                                      ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc      >  1:785216/1‑100 (MQ=255)
                                                                                       gAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCg     <  2:1290454/100‑1 (MQ=255)
                                                                                          gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa  >  1:1549725/1‑100 (MQ=255)
                                                                                          gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa  >  1:1141479/1‑100 (MQ=255)
                                                                                               |                                                                                              
CTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATA  >  NC_000913/739909‑740098

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