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.

AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGT  >  NC_000913/1207886‑1208077
                                                                                                |                                                                                                
aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtg                                                                                              >  2:436236/1‑100 (MQ=255)
 aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCtgt                                                                                             >  1:1941589/1‑100 (MQ=255)
   cAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTgtgt                                                                                           >  1:3111572/1‑100 (MQ=255)
    aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGtt                                                                                          >  2:3194314/1‑100 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                        >  2:1037978/1‑100 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGc                                                                                        <  2:1710417/100‑1 (MQ=255)
      tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGAGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGc                                                                                        <  1:1185121/100‑1 (MQ=255)
        aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCg                                                                                      <  2:294802/100‑1 (MQ=255)
         tCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgt                                                                                     <  2:2675741/100‑1 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgtt                                                                                    >  2:2168898/1‑100 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttg                                                                                   >  2:1618341/1‑100 (MQ=255)
           aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttg                                                                                   >  1:1110503/1‑100 (MQ=255)
            ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTTGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCgttgn                                                                                  <  1:1897533/100‑2 (MQ=255)
             tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATATCTGTGTTGCCgttgt                                                                                  <  1:2394380/99‑1 (MQ=255)
              cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGttgttg                                                                                <  1:945800/100‑1 (MQ=255)
               aGTGGGGGTGGGATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGc                                                                               >  1:942911/1‑100 (MQ=255)
                 tGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaa                                                                             >  1:1377223/1‑100 (MQ=255)
                  gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaaa                                                                            >  2:3270392/1‑100 (MQ=255)
                  gggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCaaa                                                                            <  2:1337689/100‑1 (MQ=255)
                   ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                           >  1:1000470/1‑100 (MQ=255)
                    gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAtt                                                                          <  1:1654294/100‑1 (MQ=255)
                    gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAAt                                                                           >  2:1036189/1‑99 (MQ=255)
                     ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTa                                                                         <  1:633541/100‑1 (MQ=255)
                         gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGAc                                                                     >  2:212197/1‑100 (MQ=255)
                          gaTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                    <  2:942923/99‑1 (MQ=255)
                          aaTGCCGCAAGAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGa                                                                      <  1:3429213/98‑1 (MQ=255)
                          aaTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                    >  1:2096425/1‑100 (MQ=255)
                            tGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAAt                                                                  >  1:1297110/1‑100 (MQ=255)
                             gCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                                                 <  1:2758479/100‑1 (MQ=255)
                               cGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTc                                                               >  2:2139871/1‑100 (MQ=255)
                               cGCAAAAGCATCTGCCTCTGGGGGGAATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATg                                                                 <  1:2978749/98‑1 (MQ=255)
                                    aaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                                                           <  2:1688818/99‑1 (MQ=255)
                                     aGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTTCAAATTATGACAATGTCAGCGCa                                                         >  1:2116460/1‑100 (MQ=255)
                                     aGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                                                         >  1:861570/1‑100 (MQ=255)
                                      gCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAc                                                        >  1:2515689/1‑100 (MQ=255)
                                       cATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       >  2:2734001/1‑100 (MQ=255)
                                       cATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       <  1:2301029/100‑1 (MQ=255)
                                       cATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAcc                                                       <  1:775528/100‑1 (MQ=255)
                                           tGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATc                                                   <  2:1020071/100‑1 (MQ=255)
                                            gCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  <  1:1013325/100‑1 (MQ=255)
                                            gCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  >  1:1492077/1‑100 (MQ=255)
                                            gCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCa                                                  <  1:3260610/100‑1 (MQ=255)
                                              ctctGGGTGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAAt                                                <  2:1222691/100‑1 (MQ=255)
                                                ctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  2:1833952/100‑1 (MQ=255)
                                                ctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  1:60018/100‑1 (MQ=255)
                                                ctGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              >  2:1978082/1‑100 (MQ=255)
                                                 tGGGGGGGATGGGAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCa                                             >  2:2282931/1‑100 (MQ=255)
                                                 tGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATACCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGc                                              <  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)
                                                                                                |                                                                                                
AATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGT  >  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: