Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,981 G→A 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,9810GA96.9% 408.5 / 0.9 130A161T (GCC→ACC) stfPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (65/61);  total (66/64)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 131 total aligned reads displayed.

AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGT  >  NC_000913/1207886‑1208077
                                                                                                |                                                                                                
aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtg                                                                                              >  2:436236/1‑100 (MQ=255)
 aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtgt                                                                                             >  1:1941589/1‑100 (MQ=255)
   cAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTgtgt                                                                                           >  1:3111572/1‑100 (MQ=255)
    aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGtt                                                                                          >  2:3194314/1‑100 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                        >  2:1037978/1‑100 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                        <  2:1710417/100‑1 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGAGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGc                                                                                        <  1:1185121/100‑1 (MQ=255)
        aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCg                                                                                      <  2:294802/100‑1 (MQ=255)
         tCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgt                                                                                     <  2:2675741/100‑1 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgtt                                                                                    >  2:2168898/1‑100 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttg                                                                                   >  2:1618341/1‑100 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttg                                                                                   >  1:1110503/1‑100 (MQ=255)
            ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑TGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttgn                                                                                  <  1:1897533/100‑2 (MQ=255)
             tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATATCTGTGTTGCCgttgt                                                                                  <  1:2394380/99‑1 (MQ=255)
              cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGttgttg                                                                                <  1:945800/100‑1 (MQ=255)
               aGTGGGGGTGGGATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                               >  1:942911/1‑100 (MQ=255)
                 tGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaa                                                                             >  1:1377223/1‑100 (MQ=255)
                  gggggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaaa                                                                            >  2:3270392/1‑100 (MQ=255)
                  gggggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaaa                                                                            <  2:1337689/100‑1 (MQ=255)
                   ggggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                           >  1:1000470/1‑100 (MQ=255)
                    gggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAtt                                                                          <  1:1654294/100‑1 (MQ=255)
                    gggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                           >  2:1036189/1‑99 (MQ=255)
                     ggTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTa                                                                         <  1:633541/100‑1 (MQ=255)
                         gAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGAc                                                                     >  2:212197/1‑100 (MQ=255)
                          gaTGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                    <  2:942923/99‑1 (MQ=255)
                          aaTGCCGCAAGAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGa                                                                      <  1:3429213/98‑1 (MQ=255)
                          aaTGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                    >  1:2096425/1‑100 (MQ=255)
                            tGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAAt                                                                  >  1:1297110/1‑100 (MQ=255)
                             gCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                                                 <  1:2758479/100‑1 (MQ=255)
                               cGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTc                                                               >  2:2139871/1‑100 (MQ=255)
                               cGCAAAAGCATCTGCCTCT‑GGGGGGAATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                                                 <  1:2978749/98‑1 (MQ=255)
                                    aaGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                                                           <  2:1688818/99‑1 (MQ=255)
                                     aGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTTCAAATTATGACAATGTCAGCGCa                                                         >  1:2116460/1‑100 (MQ=255)
                                     aGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                                         >  1:861570/1‑100 (MQ=255)
                                      gCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAc                                                        >  1:2515689/1‑100 (MQ=255)
                                       cATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       >  2:2734001/1‑100 (MQ=255)
                                       cATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       <  1:2301029/100‑1 (MQ=255)
                                       cATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       <  1:775528/100‑1 (MQ=255)
                                           tGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATc                                                   <  2:1020071/100‑1 (MQ=255)
                                            gCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  <  1:1013325/100‑1 (MQ=255)
                                            gCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  >  1:1492077/1‑100 (MQ=255)
                                            gCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  <  1:3260610/100‑1 (MQ=255)
                                              ctct‑GGGTGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAAt                                                <  2:1222691/100‑1 (MQ=255)
                                                ct‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  2:1833952/100‑1 (MQ=255)
                                                ct‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  1:60018/100‑1 (MQ=255)
                                                ct‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              >  2:1978082/1‑100 (MQ=255)
                                                 t‑GGGGGGGATGGGAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCa                                             >  2:2282931/1‑100 (MQ=255)
                                                 t‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  1:1501182/99‑1 (MQ=255)
                                                   gggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                            >  1:3226823/1‑100 (MQ=255)
                                                   gggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                            <  2:2642639/100‑1 (MQ=255)
                                                   gggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                            >  2:3454287/1‑100 (MQ=255)
                                                   gggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAg                                            >  1:2687355/1‑100 (MQ=255)
                                                    ggggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt                                           >  2:1252756/1‑100 (MQ=255)
                                                     gggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt                                           >  2:273367/1‑99 (MQ=255)
                                                      ggggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTgg                                         <  2:1169219/100‑1 (MQ=255)
                                                       gggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                                        <  2:2315425/100‑1 (MQ=255)
                                                        ggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCa                                       >  1:986829/1‑100 (MQ=255)
                                                        ggATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCa                                       >  2:1747476/1‑100 (MQ=255)
                                                             gAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggg                                  >  2:786715/1‑100 (MQ=255)
                                                              aaaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                                 >  2:1136948/1‑100 (MQ=255)
                                                              aaaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                                 <  1:3375999/100‑1 (MQ=255)
                                                              aaaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgg                                   <  1:3118412/98‑1 (MQ=255)
                                                               aaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                                <  1:2622121/100‑1 (MQ=255)
                                                               aaTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggggg                                <  2:1128812/100‑1 (MQ=255)
                                                                aTACAGTTGTATTCCGGGTTGTGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGAAACgggggg                               <  1:2725569/100‑1 (MQ=255)
                                                                aTACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggggg                               >  1:808374/1‑100 (MQ=255)
                                                                 tACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                              <  1:1221492/100‑1 (MQ=255)
                                                                 tACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                              >  2:2865940/1‑100 (MQ=255)
                                                                   cAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata                            <  2:779598/100‑1 (MQ=255)
                                                                   cAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGata                            <  2:2929239/100‑1 (MQ=255)
                                                                   cAGTTGTATTAACGATAGCGTTTAATAAAACTTGGGTTGGCGTTGTTGCAAATTATGACAATTTCAGGGCACCTATCAGTGCAGTGGCGAAGGGGGGcta                            >  1:1155620/1‑100 (MQ=255)
                                                                    aGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGatat                           >  2:2371128/1‑100 (MQ=255)
                                                                     gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                          >  2:434898/1‑100 (MQ=255)
                                                                     gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                          <  2:2020910/100‑1 (MQ=255)
                                                                     gTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                          >  2:590377/1‑100 (MQ=255)
                                                                      ttGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATAc                         <  2:2323517/100‑1 (MQ=255)
                                                                      ttGTATTCCCGGTTGCGTTTAATAATACCTGTATTGCCGTTGTTGCAAATTATGACAGTGGCAGCGCACCCATCAATGCAGTGGCAACAGGGGGat                             <  1:1180931/96‑1 (MQ=255)
                                                                        gTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACaa                       >  1:356815/1‑100 (MQ=255)
                                                                          aTTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAAcc                     >  1:724770/1‑100 (MQ=255)
                                                                          aTTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAAcc                     <  2:1519667/100‑1 (MQ=255)
                                                                          aTTCCCGGGTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAAcc                     <  2:3270679/100‑1 (MQ=255)
                                                                            tCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGTCAACGGGGGGATATACAACCAc                   <  1:2015371/100‑1 (MQ=255)
                                                                            tCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCAc                   >  1:2308926/1‑100 (MQ=255)
                                                                                ggTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCg               >  1:1365523/1‑100 (MQ=255)
                                                                                 gTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGt              <  2:3205217/100‑1 (MQ=255)
                                                                                  ttGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGtt             <  1:3242064/100‑1 (MQ=255)
                                                                                   tGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATCCAACCACTTCGttt            >  1:2369867/1‑100 (MQ=255)
                                                                                   tGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttt            <  2:1526440/100‑1 (MQ=255)
                                                                                     tccgatctaatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt          <  1:2690578/93‑1 (MQ=255)
                                                                                     cGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGttttt          >  2:831421/1‑100 (MQ=255)
                                                                                      gTTTAATAATACCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGAAGTGGCAACGGGGGGATATACAACCACTTCGttttt          >  1:2767363/1‑99 (MQ=255)
                                                                                      gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta         >  1:3248079/1‑100 (MQ=255)
                                                                                      gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTtta         <  2:604527/100‑1 (MQ=255)
                                                                                      gTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGATTtta         >  2:571154/1‑100 (MQ=255)
                                                                                        ttAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttatt       >  2:2787890/1‑100 (MQ=255)
                                                                                          aataatACCTGTGTTGCCGATGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTAc     >  1:608020/1‑100 (MQ=255)
                                                                                           ataataCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACg    >  2:779036/1‑100 (MQ=255)
                                                                                            taatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACgg   <  2:1289321/100‑1 (MQ=255)
                                                                                            taatACCGGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACgg   >  2:885508/1‑100 (MQ=255)
                                                                                             aatACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTttagatcgga  >  2:2690620/1‑93 (MQ=255)
                                                                                                |                                                                                                
AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGT  >  NC_000913/1207886‑1208077

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: ‑