Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,207,981 G→A 58.8% A161T (GCC→ACC)  stfP → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,9810GA58.8% 34.6 / 102.8 90A161T (GCC→ACC) stfPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (22/31);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

GAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACG  >  NC_000913/1207882‑1208075
                                                                                                   |                                                                                              
gAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATg                                                                                                <  2:550058/100‑1 (MQ=255)
 aGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGAAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGc                                                                                               >  2:2386421/1‑100 (MQ=255)
  gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATg                                                                                                <  1:2753689/98‑1 (MQ=255)
  gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGcc                                                                                              >  2:2092683/1‑100 (MQ=255)
  gCAAGCAGTAATCATTCAGTGGGGGTGGAATGCCGCCAAAGCATCTGCCTCTGGGGGGGATGGAAAGACAGTTGTATTCCCGGTTGCGTTTAATAATGcc                                                                                              <  1:1994317/100‑1 (MQ=255)
    aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtg                                                                                            <  1:2521819/100‑1 (MQ=255)
    aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtg                                                                                            <  1:1542650/100‑1 (MQ=255)
     aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtgt                                                                                           >  2:306666/1‑100 (MQ=255)
      tCAGTAATCATTCAGTGGGGGTGGAAGGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtgtg                                                                                          <  2:1628695/100‑1 (MQ=255)
       cAGTAATCATTCAGTGGGGGTGGAATACCGCAAAAGCATCTGCCTCTGGGAGGGATGGAAATACAGATGTATTCCCGGTTGCGTTTAATAATACCTgtgt                                                                                         >  1:645529/1‑100 (MQ=255)
        aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt                                                                                        >  1:3292742/1‑100 (MQ=255)
        aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGtt                                                                                        <  1:2217339/100‑1 (MQ=255)
         gTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTg                                                                                       >  2:665264/1‑100 (MQ=255)
          tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                      >  2:786879/1‑100 (MQ=255)
          tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                      <  2:2999590/100‑1 (MQ=255)
          tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCAGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                      >  1:1804434/1‑100 (MQ=255)
            aTCATTCAGTGGGGGGGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGcc                                                                                     <  1:3024495/99‑1 (MQ=255)
              cATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGTGGAAATACAGTTGTATTCCCGGTTGCGTTTAATACTACATGTGTTGCCgtt                                                                                  >  1:1977426/1‑100 (MQ=255)
                ccgatctggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgttgt                                                                                <  1:775657/93‑1 (MQ=255)
                 tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttgtt                                                                               >  1:2602015/1‑100 (MQ=255)
                   cagtggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                             >  1:283192/5‑100 (MQ=255)
                   aGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCATGTGTTGCCGTTGTTGc                                                                             >  2:2612854/1‑100 (MQ=255)
                   aGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                             >  1:2437141/1‑100 (MQ=255)
                      gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaaa                                                                          >  1:2610085/1‑100 (MQ=255)
                      gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaa                                                                           >  1:1952220/1‑99 (MQ=255)
                      gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGGGTTGCCGTTGTTGCaaa                                                                          >  2:1451809/1‑100 (MQ=255)
                       ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgttgtagatcgg                                                                         >  2:775660/1‑93 (MQ=255)
                       ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                         <  1:542439/100‑1 (MQ=255)
                        gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAAtt                                                                        <  2:1712800/100‑1 (MQ=255)
                         ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTa                                                                       <  1:2478013/100‑1 (MQ=255)
                         ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTa                                                                       >  2:35202/1‑100 (MQ=255)
                           tGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATg                                                                     <  1:3057255/100‑1 (MQ=255)
                           tGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATg                                                                     <  2:2038324/100‑1 (MQ=255)
                            ggAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGa                                                                    <  1:2806394/100‑1 (MQ=255)
                             gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGAc                                                                   <  2:2743842/100‑1 (MQ=255)
                               aTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                                                                 <  2:983568/100‑1 (MQ=255)
                                 gCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                                               >  2:331619/1‑100 (MQ=255)
                                  ccGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGt                                                              >  1:2352553/1‑100 (MQ=255)
                                   cGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTc                                                             <  1:549504/100‑1 (MQ=255)
                                    gCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTc                                                             <  2:1231985/99‑1 (MQ=255)
                                      aaaaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                          <  1:569834/100‑1 (MQ=255)
                                      aaaaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc                                                          >  2:2787231/1‑100 (MQ=255)
                                       aaggCATCTGCGTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                                                         <  1:2623/97‑1 (MQ=255)
                                        aaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcgc                                                        <  2:2032038/100‑1 (MQ=255)
                                         aGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                                       <  1:1925080/100‑1 (MQ=255)
                                          gCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAc                                                      <  2:775060/100‑1 (MQ=255)
                                           cATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                     <  2:2602035/100‑1 (MQ=255)
                                           cATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                     <  1:2787208/100‑1 (MQ=255)
                                           cATCTGCCTCGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                     <  1:2010630/100‑1 (MQ=255)
                                             tCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTa                                                   <  1:288017/100‑1 (MQ=255)
                                             tCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTATAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTa                                                   >  1:1876848/1‑100 (MQ=255)
                                              cTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTAt                                                  >  2:2197327/1‑100 (MQ=255)
                                               tGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATc                                                 >  1:1691777/1‑100 (MQ=255)
                                               tGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTAtn                                                 <  1:1685060/100‑2 (MQ=255)
                                                   tctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGGGTGGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATg                                             >  2:631109/1‑100 (MQ=255)
                                                   tctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATg                                             <  2:3197061/100‑1 (MQ=255)
                                                     tGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                            <  1:2406270/99‑1 (MQ=255)
                                                        gggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTTATAATACATGTGATGCCGTTGTTGCAAATTATGACAATGTAAGCGCACCTATGAATGCAGTg                                        >  1:3206310/1‑100 (MQ=255)
                                                        gggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                                        >  1:998399/1‑100 (MQ=255)
                                                        gggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTg                                        >  2:1809309/1‑100 (MQ=255)
                                                         ggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTgg                                       <  1:974439/100‑1 (MQ=255)
                                                          gggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                                      <  1:3038801/100‑1 (MQ=255)
                                                           ggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                                      <  1:620354/99‑1 (MQ=255)
                                                              tGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACg                                  >  2:3055547/1‑100 (MQ=255)
                                                              tGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCaagg                                  >  1:2534021/1‑98 (MQ=255)
                                                               ggAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgg                                 >  2:992149/1‑100 (MQ=255)
                                                                  aaTACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                              <  2:2617313/100‑1 (MQ=255)
                                                                    tACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                            <  2:2931624/100‑1 (MQ=255)
                                                                    tACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                            <  2:2417351/100‑1 (MQ=255)
                                                                    tACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                            <  1:1038957/100‑1 (MQ=255)
                                                                      cAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata                          >  2:2500230/1‑100 (MQ=255)
                                                                        gTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                        <  2:2072638/100‑1 (MQ=255)
                                                                        gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                        >  1:3042467/1‑100 (MQ=255)
                                                                          tGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACa                      <  2:1052603/100‑1 (MQ=255)
                                                                            tATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAAc                    <  1:819772/100‑1 (MQ=255)
                                                                             aTTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAAcc                   >  2:1004312/1‑100 (MQ=255)
                                                                               tCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAAAGGGGGGATATACAACCAc                 >  2:1668727/1‑100 (MQ=255)
                                                                                cccGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACt                <  1:2606536/100‑1 (MQ=255)
                                                                                     ttGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGtt           <  2:2847113/100‑1 (MQ=255)
                                                                                        cGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt        >  1:511919/1‑100 (MQ=255)
                                                                                        cGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt        <  1:737124/100‑1 (MQ=255)
                                                                                         ntttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta       >  1:95965/2‑100 (MQ=255)
                                                                                         gTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGtttt         >  1:1841141/1‑98 (MQ=255)
                                                                                         gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta       <  2:60421/100‑1 (MQ=255)
                                                                                         gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta       <  1:3111814/100‑1 (MQ=255)
                                                                                          tttAATAATAGCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAACGCACCTATCAATGCAGTGGCAACGGTGGGCTATACGACCACTTCGtttttt       <  1:790000/99‑2 (MQ=255)
                                                                                          tttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttat      >  2:486025/1‑100 (MQ=255)
                                                                                          tttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttat      <  1:2930162/100‑1 (MQ=255)
                                                                                           ttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatt     >  1:1586917/1‑100 (MQ=255)
                                                                                           ttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatn     <  1:1400629/100‑2 (MQ=255)
                                                                                            taataaTGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatta    >  2:2950712/1‑100 (MQ=255)
                                                                                             aataatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTAc   <  1:310497/100‑1 (MQ=255)
                                                                                              ataataCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACg  >  1:728091/1‑100 (MQ=255)
                                                                                                   |                                                                                              
GAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACG  >  NC_000913/1207882‑1208075

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: