Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 366,290 C→T D6N (GAT→AAT)  lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,2900CT99.4% 88939.8 / 90.6 26327D6N (GAT→AAT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (73/82);  new base T (12978/13188);  total (13054/13273)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 26333 total aligned reads displayed.

ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366192‑366388
                                                                                                  |                                                                                                  
  tgcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCg                                                                                                  <  2:793455‑M2/7‑1 (MQ=255)
   gcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGt                                                                                                 <  2:2304044‑M2/8‑1 (MQ=255)
    caccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTa                                                                                                >  2:1397124‑M2/91‑99 (MQ=255)
     accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                               <  1:3668622‑M2/10‑1 (MQ=255)
      ccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAt                                                                                              <  2:2900288‑M2/11‑1 (MQ=255)
        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  1:3252866‑M2/13‑1 (MQ=255)
        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  2:1184947‑M2/13‑1 (MQ=255)
         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           >  1:904326‑M2/86‑99 (MQ=255)
         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           >  2:2977050‑M2/86‑99 (MQ=255)
         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           <  1:1955781‑M2/14‑1 (MQ=255)
         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           <  1:1889583‑M2/14‑1 (MQ=255)
         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           >  1:2575621‑M2/86‑99 (MQ=255)
            ccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCatgtc                                                                                        >  1:970963‑M2/83‑97 (MQ=255)
               gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                    >  1:18193‑M2/80‑99 (MQ=255)
               gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                    >  2:1539784‑M2/80‑99 (MQ=255)
                cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                   <  2:3603085‑M2/21‑1 (MQ=255)
                  cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                 >  2:3252059‑M2/77‑99 (MQ=255)
                  cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                 >  2:2100144‑M2/77‑99 (MQ=255)
                   gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                                >  2:2437981‑M2/76‑99 (MQ=255)
                    agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                               <  2:2838161‑M2/25‑1 (MQ=255)
                    agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                               <  2:721816‑M2/25‑1 (MQ=255)
                     gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                              <  1:2919152‑M2/26‑1 (MQ=255)
                      cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGttt                                                                             <  1:2780781‑M2/27‑1 (MQ=255)
                       gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTc                                                                            <  1:3529976‑M2/28‑1 (MQ=255)
                       gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTc                                                                            <  2:1749565‑M2/28‑1 (MQ=255)
                         gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCt                                                                          <  2:3675182‑M2/30‑1 (MQ=255)
                          ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                         >  1:2436678‑M2/69‑99 (MQ=255)
                           tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                        <  1:2996354‑M2/32‑1 (MQ=255)
                             cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                      >  2:3121386‑M2/66‑99 (MQ=255)
                             cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                      <  1:768067‑M2/34‑1 (MQ=255)
                              gtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgtg                                                                     >  1:2301477‑M2/65‑99 (MQ=255)
                                gccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaa                                                                   <  1:1131227‑M2/37‑1 (MQ=255)
                                  cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                 <  1:2160987‑M2/39‑1 (MQ=255)
                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              <  1:2716268‑M2/42‑1 (MQ=255)
                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              <  1:3323810‑M2/42‑1 (MQ=255)
                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              >  2:1330444‑M2/58‑99 (MQ=255)
                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              >  1:5725‑M2/58‑99 (MQ=255)
                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              <  1:127034‑M2/42‑1 (MQ=255)
                                      tgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGtt                                                             <  2:1040151‑M2/43‑1 (MQ=255)
                                      tAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              <  1:1481311/99‑1 (MQ=255)
                                        tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                           >  2:3322093‑M2/55‑99 (MQ=255)
                                         gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                          <  1:2090711‑M2/46‑1 (MQ=255)
                                          aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                         >  1:83878‑M2/53‑99 (MQ=255)
                                           agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                        >  2:3745788‑M2/52‑99 (MQ=255)
                                           agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                        >  1:443334‑M2/52‑99 (MQ=255)
                                            gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                       <  2:859802‑M2/49‑1 (MQ=255)
                                            gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                       >  1:2224364‑M2/51‑99 (MQ=255)
                                            gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                       <  2:516727‑M2/49‑1 (MQ=255)
                                              aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                     <  1:3190229‑M2/51‑1 (MQ=255)
                                               aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                    <  1:509768‑M2/52‑1 (MQ=255)
                                               aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                    >  1:640113‑M2/48‑99 (MQ=255)
                                                aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTNAAATTGTTATCCGCTcac                                                   >  1:2368544‑M2/47‑99 (MQ=255)
                                                 agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                  >  2:223532‑M2/46‑99 (MQ=255)
                                                   tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                <  2:2368839‑M2/56‑1 (MQ=255)
                                                     aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                              <  2:3458655‑M2/58‑1 (MQ=255)
                                                     aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                              >  1:1754608‑M2/42‑99 (MQ=255)
                                                      atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCATAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                             <  2:1253702‑M2/59‑1 (MQ=255)
                                                       tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                            >  1:3067321‑M2/40‑99 (MQ=255)
                                                         atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                          <  2:583210‑M2/62‑1 (MQ=255)
                                                         atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                          >  1:1196106‑M2/38‑99 (MQ=255)
                                                          tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                         >  1:1689086‑M2/37‑99 (MQ=255)
                                                             agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                      >  2:3530938‑M2/34‑99 (MQ=255)
                                                             agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                      <  1:1413702‑M2/66‑1 (MQ=255)
                                                              gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                     >  1:3811584‑M2/33‑99 (MQ=255)
                                                              gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                     >  1:1348356‑M2/33‑99 (MQ=255)
                                                                 gtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                  <  1:2021999‑M2/70‑1 (MQ=255)
                                                                  tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  2:660505‑M2/29‑99 (MQ=255)
                                                                  tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  2:928204‑M2/29‑99 (MQ=255)
                                                                  tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 <  1:1546483‑M2/71‑1 (MQ=255)
                                                                   caagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGa                                <  2:1962349‑M2/72‑1 (MQ=255)
                                                                    aagttgtcatgattgttaccgagctcGAATTCGTAATGATGTCATTGCTGTTTCCTGTGTGAATTTGTTATCCGCTCACAATTCCACACAACATACGAg                               <  2:3390296‑M2/73‑1 (MQ=255)
                                                                    aagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                               >  2:2034661‑M2/27‑99 (MQ=255)
                                                                    aagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                               <  2:447943‑M2/73‑1 (MQ=255)
                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              <  1:3598879‑M2/74‑1 (MQ=255)
                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              <  2:3044238‑M2/74‑1 (MQ=255)
                                                                       ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                            >  1:699880‑M2/24‑99 (MQ=255)
                                                                       ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                            <  1:1617814‑M2/76‑1 (MQ=255)
                                                                       ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                            >  2:2701879‑M2/24‑99 (MQ=255)
                                                                         gtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGa                          >  1:3393302‑M2/22‑99 (MQ=255)
                                                                          tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                         >  1:1002747‑M2/21‑99 (MQ=255)
                                                                          tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                         >  2:2501718‑M2/21‑99 (MQ=255)
                                                                           cataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                        <  2:3815602‑M2/80‑1 (MQ=255)
                                                                            ataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGc                       >  2:217080‑M2/19‑99 (MQ=255)
                                                                             taattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                      >  2:2634210‑M2/18‑99 (MQ=255)
                                                                               attggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    <  2:1466162‑M2/84‑1 (MQ=255)
                                                                                 tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                  >  1:3459653‑M2/14‑99 (MQ=255)
                                                                                  ggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                 <  1:3129546‑M2/87‑1 (MQ=255)
                                                                                   gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                <  1:2637753‑M2/88‑1 (MQ=255)
                                                                                   gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                <  2:1682478‑M2/88‑1 (MQ=255)
                                                                                    taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg               >  2:1822515‑M2/11‑99 (MQ=255)
                                                                                     accgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt              >  2:189359‑M2/10‑99 (MQ=255)
                                                                                      ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa             >  1:2504849‑M2/9‑99 (MQ=255)
                                                                                      ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa             <  1:1066847‑M2/91‑1 (MQ=255)
                                                                                       cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa            >  2:3185678‑M2/8‑99 (MQ=255)
                                                                                       cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa            <  2:1610558‑M2/92‑1 (MQ=255)
                                                                                        gagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa           <  2:1114929‑M2/93‑1 (MQ=255)
                                                                                        gagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa           <  2:2571263‑M2/93‑1 (MQ=255)
                                                                                           ctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc        <  2:2772876‑M2/96‑1 (MQ=255)
                                                                                               aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg    >  1:3740866/1‑99 (MQ=255)
                                                                                               aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg    >  1:284558/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366192‑366388

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