Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,221,645 T→G E281D (GAA→GAC arpA ← ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,221,6450TG97.9% 841.6 / 1.6 241E281D (GAA→GACarpAankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (120/115);  total (124/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

CAGATTGATTAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGC  >  NC_000913/4221548‑4221741
                                                                                                 |                                                                                                
cAGATTGATTAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTc                                                                                                <  2:1726548/100‑1 (MQ=255)
  gattgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGc                                                                                              >  1:518461/1‑100 (MQ=255)
   attgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATACAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCt                                                                                             >  1:3837915/1‑100 (MQ=255)
   attgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCt                                                                                             <  1:1077355/100‑1 (MQ=255)
   attgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCt                                                                                             <  1:3612034/100‑1 (MQ=255)
   attgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCt                                                                                             <  2:333972/100‑1 (MQ=255)
     tgattAATAACATGATCATTTGACCATTTCAGGAGTCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCa                                                                                           <  2:5418565/100‑1 (MQ=255)
     tgattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATGCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCa                                                                                           >  2:5372114/1‑100 (MQ=255)
      gattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCTCTCAGAGAATAGTTTTTAATATTCTCTTTGCTGTCGCTCAt                                                                                          >  2:4685729/1‑100 (MQ=255)
      gattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCAt                                                                                          <  2:2958391/100‑1 (MQ=255)
      gattAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCAt                                                                                          <  1:589552/100‑1 (MQ=255)
       attAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATa                                                                                         <  2:2362487/100‑1 (MQ=255)
       attAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATa                                                                                         >  2:3021301/1‑100 (MQ=255)
       attAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATa                                                                                         >  2:3926/1‑100 (MQ=255)
         taataaCATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAAt                                                                                       >  1:227334/1‑100 (MQ=255)
           ataaCATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTc                                                                                     <  2:2310371/100‑1 (MQ=255)
            taaCATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTcn                                                                                    <  1:1359820/100‑2 (MQ=255)
            taaCATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCt                                                                                    <  1:725082/100‑1 (MQ=255)
            taaCATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCt                                                                                    >  1:3716044/1‑100 (MQ=255)
                aagATCATATGACCATTTCAGGAGGCTATAATCTAGATTCATAAAATATAACAGAAAATAGTGTTTAATATTCTATCTGCTGTCGCTCATACTTCTTAtt                                                                                >  2:2020424/3‑100 (MQ=255)
                aTGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTAtt                                                                                <  2:4250125/100‑1 (MQ=255)
                 tGATCATTTGACCATTTCAGGAGGCTATAATCTACATTTATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTAttt                                                                               <  1:5321266/100‑1 (MQ=255)
                  gATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTAtttt                                                                              >  1:4206006/1‑100 (MQ=255)
                  gATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTAtttt                                                                              >  2:4347447/1‑100 (MQ=255)
                   aTCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTa                                                                             <  2:2568789/100‑1 (MQ=255)
                      aTTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCg                                                                          <  2:3723936/100‑1 (MQ=255)
                      aTTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCg                                                                          <  2:1476127/100‑1 (MQ=255)
                        ttGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCgtt                                                                        <  1:4275459/100‑1 (MQ=255)
                           aCCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCgttagt                                                                     >  1:3474061/1‑100 (MQ=255)
                           aCCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCgttagt                                                                     <  2:1308172/100‑1 (MQ=255)
                            ccATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTTCTGTCGCTCATAATTCTTATTTTATCtttagtt                                                                    >  1:1651565/1‑100 (MQ=255)
                             cATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCgttagtta                                                                   <  1:4250729/100‑1 (MQ=255)
                             cATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCgttagtta                                                                   >  1:4558441/1‑100 (MQ=255)
                              aTTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTAc                                                                  <  2:5112836/100‑1 (MQ=255)
                               tttCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACa                                                                 <  2:4143603/100‑1 (MQ=255)
                                ttCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACAt                                                                <  1:2880860/100‑1 (MQ=255)
                                 tCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACAtt                                                               >  2:455579/1‑100 (MQ=255)
                                  cAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTa                                                              <  1:4848118/100‑1 (MQ=255)
                                  cAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTa                                                              <  2:3175099/100‑1 (MQ=255)
                                   aggaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAg                                                             >  1:4135891/1‑100 (MQ=255)
                                    ggaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAgg                                                            >  1:3373353/1‑100 (MQ=255)
                                     gaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGa                                                           <  2:4715105/100‑1 (MQ=255)
                                     gaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGa                                                           <  1:543024/100‑1 (MQ=255)
                                     gaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGa                                                           <  1:2119846/100‑1 (MQ=255)
                                     gaggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGa                                                           >  1:5128930/1‑100 (MQ=255)
                                      aggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGAt                                                          <  2:1463785/100‑1 (MQ=255)
                                      aggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGAt                                                          <  2:36223/100‑1 (MQ=255)
                                      aggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGAt                                                          <  1:4022708/100‑1 (MQ=255)
                                       ggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATc                                                         >  1:3039111/1‑100 (MQ=255)
                                       ggCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATCCTTATTTTATCGTTAGTTACATTAGGATc                                                         >  2:2287642/1‑100 (MQ=255)
                                         cTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATTGTTAGTTACATTAGGATCAt                                                       >  2:2261594/1‑100 (MQ=255)
                                         cTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCAt                                                       >  2:1055960/1‑100 (MQ=255)
                                           ataATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATaa                                                     >  1:577237/1‑100 (MQ=255)
                                           ataATCTACATTCATAAAATCCATAAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAAGTCTTATTTTATCGTTAGTTACATTAGGATCATaa                                                     <  2:1258356/100‑1 (MQ=255)
                                               tCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCg                                                 >  2:1514845/1‑100 (MQ=255)
                                               tCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCg                                                 <  1:398925/100‑1 (MQ=255)
                                                   cATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAAc                                             <  2:871907/100‑1 (MQ=255)
                                                   cATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAAc                                             <  1:169119/100‑1 (MQ=255)
                                                    aTTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAAcc                                            <  2:1110173/100‑1 (MQ=255)
                                                     ttCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCa                                           <  1:1089695/100‑1 (MQ=255)
                                                      tCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCGTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCAc                                          <  1:4740434/100‑1 (MQ=255)
                                                        aTAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGt                                        <  1:454951/100‑1 (MQ=255)
                                                         tAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTa                                       >  2:4925045/1‑100 (MQ=255)
                                                           aaaTCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAAt                                     <  2:2111342/100‑1 (MQ=255)
                                                           aaaTCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAAt                                     <  2:706977/100‑1 (MQ=255)
                                                           aaaTCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAAt                                     >  2:1776907/1‑100 (MQ=255)
                                                               ccATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGt                                 >  1:5437308/1‑100 (MQ=255)
                                                               ccATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGt                                 <  1:1828313/100‑1 (MQ=255)
                                                               ccATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGt                                 <  1:2333507/100‑1 (MQ=255)
                                                               ccATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGt                                 <  1:3330653/100‑1 (MQ=255)
                                                                cATCAGAGAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTa                                <  2:5031916/100‑1 (MQ=255)
                                                                    agagAATAGTGTTTAATATTCTGTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATc                            >  1:4666810/1‑100 (MQ=255)
                                                                    agagAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATc                            >  2:840483/1‑100 (MQ=255)
                                                                    agagAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATc                            <  2:609405/100‑1 (MQ=255)
                                                                    agagAATAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATc                            >  2:3062225/1‑100 (MQ=255)
                                                                        aaTAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTc                        <  2:4334175/100‑1 (MQ=255)
                                                                        aaTAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTc                        <  1:1534066/100‑1 (MQ=255)
                                                                          tAGTGTTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAg                      <  2:1829411/100‑1 (MQ=255)
                                                                            gtgtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGtt                    >  2:2446868/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  1:4094000/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  1:3409987/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  2:350935/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  2:4626871/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  1:1985145/1‑100 (MQ=255)
                                                                              gtTTAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTtgt                  >  2:4806181/1‑100 (MQ=255)
                                                                                ttAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTgtgt                <  1:2404495/100‑1 (MQ=255)
                                                                                ttAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTgtgt                >  2:3771301/1‑100 (MQ=255)
                                                                                 tAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGtt               >  2:4005258/1‑100 (MQ=255)
                                                                                 tAATATTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGtt               >  2:3227159/1‑100 (MQ=255)
                                                                                   atatTCTCTTTGGTGTAGCTCATAATTCGTATTGTATCGTTAGTTACATTAGTATCATAATCCGAAACCGCGTAATGGGTACATCCTTCAGTTGTGTTCt             <  2:4219364/100‑1 (MQ=255)
                                                                                   atatTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCt             >  1:4061128/1‑100 (MQ=255)
                                                                                    tatTCTCTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCtt            >  2:4477595/1‑100 (MQ=255)
                                                                                       tctcTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTAt         >  1:739779/1‑100 (MQ=255)
                                                                                       tctcTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTAt         >  1:1553681/1‑100 (MQ=255)
                                                                                       tctcTTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTAt         <  2:1631847/100‑1 (MQ=255)
                                                                                         tctTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTc       >  2:5293282/1‑100 (MQ=255)
                                                                                         tctTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTc       <  1:933596/100‑1 (MQ=255)
                                                                                         tctTTGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTc       >  2:2650860/1‑100 (MQ=255)
                                                                                             tGCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAg   <  2:4883275/100‑1 (MQ=255)
                                                                                              gCTGTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGc  <  2:1856763/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CAGATTGATTAATAACATGATCATTTGACCATTTCAGGAGGCTATAATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGC  >  NC_000913/4221548‑4221741

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