Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,473 T→C 15.9% L61L (TTG→CTG) ‡ uspD → stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4730TC15.9% 490.5 / 79.1 218L61L (TTG→CTG) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (95/88);  new base C (20/15);  total (115/103)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-01

Only 100 of 229 total aligned reads displayed.

AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAA  >  NC_000913/4113374‑4113571
                                                                                                   |                                                                                                  
     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               <  2:717977/100‑1 (MQ=255)
     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               >  2:1878400/1‑100 (MQ=255)
     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               <  1:1487258/100‑1 (MQ=255)
     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               >  1:1125545/1‑100 (MQ=255)
      aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               <  2:1351995/99‑1 (MQ=255)
        tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                            <  1:733612/100‑1 (MQ=255)
           cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                         >  2:2786572/1‑100 (MQ=255)
           cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                         >  2:3015927/1‑100 (MQ=255)
           cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAAt                                                                                         <  1:2181824/100‑1 (MQ=255)
             cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACATCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATaa                                                                                       >  1:697973/1‑100 (MQ=255)
             cTCACCTGACATTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                       <  2:3959634/100‑2 (MQ=255)
              tctccTGTCGTTGATTCATTTTTATGTTGGTTTAAGCGAGTTGTATCCTGGCATCTACTTCCCTGCATCATTTTATATTTTTCAATTGTTGCAGAATa                                                                                        <  1:889831/95‑1 (MQ=255)
              tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCAGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                      <  2:2923078/100‑1 (MQ=255)
               caacTGACGTTAATTCATAGTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACGTCCTTTAAACAGAAGATATGCTTCAACAATTCAAGAATa                                                                                        <  1:2495134/94‑1 (MQ=255)
               cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                     <  2:3870617/100‑4 (MQ=38)
               cACCTGACGTTAATTCATATTGATGAAAGCATAAACTAGTTGTATCCGGGTATCTACTTACCTGCAACAGAAGCTATTCTTCAACAATTGAAGAATatc                                                                                      <  1:2488473/99‑3 (MQ=25)
                aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                    <  2:4047573/100‑1 (MQ=255)
                aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                    >  2:3589699/1‑100 (MQ=255)
                aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                     <  1:3842384/99‑4 (MQ=38)
                 ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                   <  1:611450/100‑1 (MQ=255)
                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                  <  2:3898129/100‑1 (MQ=255)
                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                  <  2:1017226/100‑1 (MQ=255)
                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                 >  1:816380/1‑100 (MQ=255)
                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                 <  1:2923596/100‑1 (MQ=255)
                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                 <  2:2272708/100‑8 (MQ=37)
                    gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                <  2:3801670/100‑1 (MQ=255)
                      cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                              >  2:4194421/1‑100 (MQ=255)
                        ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                            <  1:2171204/100‑1 (MQ=255)
                          aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                          >  2:347394/1‑100 (MQ=255)
                          aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATAGTTGAAGAATAAGTCGGATAACAAg                                                                          >  2:604228/1‑100 (MQ=255)
                              cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTAAAGAATAAGTCGGATAACAAGCTGt                                                                      >  1:383810/1‑100 (MQ=255)
                                tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                    >  2:54108/1‑100 (MQ=255)
                                 atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                   <  1:1942226/100‑1 (MQ=255)
                                  ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtatag                                                                  >  2:2378603/1‑99 (MQ=255)
                                   tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGGTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                   <  1:192963/98‑1 (MQ=255)
                                    gatgatGGCTTAAGCGAGTTGTATCGGGGTATCTACTTCCCTTCAACAGAAGATATTCGTCAATTGTTAAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                >  1:1651893/1‑100 (MQ=255)
                                    gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                >  2:655176/1‑100 (MQ=255)
                                    gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                >  2:1100691/1‑100 (MQ=255)
                                    gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCGGTATAAAc                                                                >  2:1869366/1‑100 (MQ=255)
                                     atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                               <  2:1988761/100‑1 (MQ=255)
                                      tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                              >  2:1430845/1‑100 (MQ=255)
                                      tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                              >  1:2399915/1‑100 (MQ=255)
                                      tgatgGCTTAAGCGAGTTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                              <  1:306043/100‑1 (MQ=255)
                                      tgatgGCTTAAGCAAGTTGTATCCGGGTATCGACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                              >  2:2973602/1‑100 (MQ=255)
                                        atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCGTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAc                                                            >  2:3106187/1‑100 (MQ=255)
                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                         >  1:2778934/1‑100 (MQ=255)
                                            cTTGAGCGAGTTGTGTTCGAGTATCCAGGTCCCGGCTACAGAAGATATTCTGCAATTGTTGATGAGGAAGTCGGATAACAAGCGGTATAAAATGACGaa                                                         >  1:2564403/1‑99 (MQ=255)
                                            cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                        <  1:474638/100‑1 (MQ=255)
                                             ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                       <  2:2562280/100‑1 (MQ=255)
                                               aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                     <  2:4177458/100‑1 (MQ=255)
                                                 gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACCAACTGTATAAACTGACGAAAAatat                                                   >  1:3798736/1‑100 (MQ=255)
                                                 gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                   >  2:3754071/1‑100 (MQ=255)
                                                  cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                  >  1:4169342/1‑100 (MQ=255)
                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGGATAATATTc                                                 <  1:3443081/100‑1 (MQ=255)
                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                 <  2:2427826/100‑1 (MQ=255)
                                                    aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                <  1:2701704/100‑1 (MQ=255)
                                                     gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                               >  2:1759192/1‑100 (MQ=255)
                                                       tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                             >  1:993171/1‑100 (MQ=255)
                                                         tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                           <  2:1217185/100‑1 (MQ=255)
                                                          aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                          <  1:1688034/100‑1 (MQ=255)
                                                           tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                         >  2:1492023/1‑100 (MQ=255)
                                                            gcGGGTATCTACTTCCCTGCACCCGAAGAGATTGTCCAATATTTTTTCAATCATTCGGATAACAAGCTGTAACAACTGACGAGAAGTATTCAAt                                              >  1:3694856/2‑94 (MQ=255)
                                                            ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAAAAAAAAAAAATATTCAATGGCCGa                                        <  1:2425424/100‑1 (MQ=255)
                                                              gggTATCTAGTTCCATGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                      <  1:353229/100‑1 (MQ=255)
                                                               ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                     >  2:3698419/1‑100 (MQ=255)
                                                                     tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                               <  1:1100687/100‑1 (MQ=255)
                                                                        ttCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                            <  1:2122986/100‑1 (MQ=255)
                                                                         tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                           >  1:3134699/1‑100 (MQ=255)
                                                                          cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                          >  2:3180778/1‑100 (MQ=255)
                                                                            cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                        >  2:2729010/1‑100 (MQ=255)
                                                                            cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                        >  1:2327610/1‑100 (MQ=255)
                                                                             tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                       <  2:1639914/100‑1 (MQ=255)
                                                                             tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGCATAAACAAAAGAAAAATATACAATGGCAGAAGACAAAACAGCg                           >  1:1515974/1‑96 (MQ=255)
                                                                              gCAACAGAAGATATTCTTGAATTGTTTAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                      <  1:1807791/100‑1 (MQ=255)
                                                                                   agaagaTATTCTTGAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                 <  1:1069173/100‑1 (MQ=255)
                                                                                   agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                 <  2:2650462/100‑1 (MQ=255)
                                                                                    gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                <  1:2584328/100‑1 (MQ=255)
                                                                                    cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                >  1:2005609/6‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                >  1:2043207/5‑99 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               >  1:1217439/5‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               >  1:3234989/5‑100 (MQ=38)
                                                                                      agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga              >  2:851494/1‑100 (MQ=255)
                                                                                      acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga              >  2:3311468/4‑100 (MQ=38)
                                                                                       cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag             >  1:3664195/3‑100 (MQ=255)
                                                                                       cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag             >  1:4199543/3‑100 (MQ=255)
                                                                                         tatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa           >  1:3879378/1‑100 (MQ=255)
                                                                                         tatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAAGCGGAGaa           >  1:4070321/1‑100 (MQ=255)
                                                                                         tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa           >  2:4256140/1‑100 (MQ=255)
                                                                                          atTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa          >  1:4357191/1‑100 (MQ=255)
                                                                                            tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg        <  1:1878394/100‑1 (MQ=255)
                                                                                            tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGAGTATTGAACGCGGAGAAATg        <  1:1525884/100‑1 (MQ=255)
                                                                                             cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg        >  2:476754/1‑99 (MQ=255)
                                                                                              ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc      <  2:1676117/100‑1 (MQ=255)
                                                                                               tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg     <  1:2820547/100‑1 (MQ=255)
                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg    >  1:4322099/1‑100 (MQ=255)
                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg    >  2:2096871/1‑100 (MQ=255)
                                                                                                 aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcagaa   >  1:867063/1‑96 (MQ=255)
                                                                                                 aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa   <  2:3372835/100‑1 (MQ=255)
                                                                                                 aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa   <  1:3410365/100‑1 (MQ=255)
                                                                                                  aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa  >  1:4244273/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAA  >  NC_000913/4113374‑4113571

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