Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,631,832 (T)6→5 coding (259/1467 nt) guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,631,8270T.99.3% 649.2 / ‑2.7 143coding (264/1467 nt)guaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (63/79);  total (63/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 164 total aligned reads displayed.

CAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCA  >  NC_000913/2631736‑2631914
                                                                                           |                                                                                       
cAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtt                                                                                        <  1:3660113/93‑1 (MQ=255)
  aaaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttt                                                                                      <  1:3881998/93‑1 (MQ=255)
   aaCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGttttt                                                                                     >  1:1249557/1‑93 (MQ=255)
    aCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttc                                                                                    <  1:4208745/93‑2 (MQ=255)
    aCCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttc                                                                                    <  1:965213/93‑2 (MQ=255)
      cGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcac                                                                                  <  1:4360919/93‑4 (MQ=255)
      cGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcac                                                                                  <  1:4713986/93‑4 (MQ=255)
      cGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcac                                                                                  >  1:1433021/1‑90 (MQ=255)
      cGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcac                                                                                  <  2:1664169/93‑4 (MQ=255)
       gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                 <  2:870137/93‑5 (MQ=255)
       gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                 >  2:342754/1‑89 (MQ=255)
       gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                 <  1:2694145/93‑5 (MQ=255)
       gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGtttttcaca                                                                                 <  2:1627391/93‑5 (MQ=255)
       gTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGAAGAACAGTCTGCGGTTCAGTCACCACACCAAATTCGTGtttttcaca                                                                                 >  2:1080872/1‑89 (MQ=255)
         tACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACg                                                                              >  2:4293727/1‑93 (MQ=255)
         tACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACg                                                                              <  1:2915136/93‑1 (MQ=255)
          aCGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACgg                                                                             <  1:696375/93‑1 (MQ=255)
            gCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCg                                                                           >  2:73553/1‑93 (MQ=255)
               cGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAAc                                                                        <  2:136923/93‑1 (MQ=255)
                ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAAct                                                                       <  2:1032646/93‑1 (MQ=255)
                ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAAct                                                                       >  1:4509223/1‑93 (MQ=255)
                ggTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAAct                                                                       <  2:4366138/93‑1 (MQ=255)
                 gTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActt                                                                      <  1:4071183/93‑1 (MQ=255)
                 gTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActt                                                                      >  2:654930/1‑93 (MQ=255)
                 gTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActt                                                                      <  2:4081605/93‑1 (MQ=255)
                  tCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtc                                                                     <  2:4442187/93‑1 (MQ=255)
                  tCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtc                                                                     <  1:1679170/93‑1 (MQ=255)
                  tCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtc                                                                     <  2:1460582/93‑1 (MQ=255)
                  tCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtc                                                                     <  2:824192/93‑1 (MQ=255)
                   cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtct                                                                    <  1:3259319/93‑1 (MQ=255)
                   cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtct                                                                    <  2:867196/93‑1 (MQ=255)
                   cAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtct                                                                    >  1:3990639/1‑93 (MQ=255)
                    aGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAActtctt                                                                   <  1:4544317/93‑1 (MQ=255)
                     gTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACttcttc                                                                  >  2:4085937/1‑93 (MQ=255)
                     gTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACAAGATTCGTGTTTTTCACACGGCGAACttcttc                                                                  >  1:1757121/1‑93 (MQ=255)
                      ttCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTtcttct                                                                 >  2:370354/1‑93 (MQ=255)
                       tCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTg                                                                <  1:378368/93‑1 (MQ=255)
                           tCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGCCGAACTTCTTCTGCCTg                                                            >  1:1258452/1‑93 (MQ=255)
                            cACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTgg                                                           <  1:2185498/93‑1 (MQ=255)
                              cTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCg                                                         >  1:1040703/1‑93 (MQ=255)
                               ttCGCGGAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGt                                                        >  2:4667428/1‑93 (MQ=255)
                                   cgcAGCGTCGTGGTTGGCAGAACAGTCTGCTGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCaa                                                    <  2:3209363/93‑1 (MQ=255)
                                   cgcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCaa                                                    <  1:575494/93‑1 (MQ=255)
                                    gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                   <  1:2407850/93‑1 (MQ=255)
                                    gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                   <  1:207721/93‑1 (MQ=255)
                                    gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                   <  1:1786904/93‑1 (MQ=255)
                                    gcAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAAt                                                   <  2:2720035/93‑1 (MQ=255)
                                     cAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATg                                                  >  1:3240905/1‑93 (MQ=255)
                                     cAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATg                                                  <  1:1917033/93‑1 (MQ=255)
                                      aGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATg                                                  <  1:3673283/93‑1 (MQ=255)
                                      aGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATgg                                                 >  2:3316166/1‑93 (MQ=255)
                                      aGCGTCGCGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATgg                                                 >  1:1721333/1‑93 (MQ=255)
                                        cgtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGAc                                               >  2:902485/1‑93 (MQ=255)
                                         gtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACa                                              >  1:3796351/1‑93 (MQ=255)
                                         gtcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACa                                              >  2:1521740/1‑93 (MQ=255)
                                          tcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACTGCGAACTTCTTCTGCCTGGCGTTCAATGGACAt                                             >  2:4188070/1‑93 (MQ=255)
                                          tcgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACAt                                             >  2:3583501/1‑93 (MQ=255)
                                           cgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATg                                            >  2:664359/1‑93 (MQ=255)
                                           cgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATg                                            >  1:237946/1‑93 (MQ=255)
                                           cgtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATg                                            <  1:1422385/93‑1 (MQ=255)
                                            gtGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGt                                           >  1:1421617/1‑93 (MQ=255)
                                               gtgagCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGtttt                                        <  2:2749478/89‑1 (MQ=255)
                                                   gCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTtgtg                                    <  1:907886/93‑1 (MQ=255)
                                                      gAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGAt                                 <  1:4574716/93‑1 (MQ=255)
                                                      gAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGAt                                 <  2:2464037/93‑1 (MQ=255)
                                                       aaCAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATa                                <  2:3753408/93‑1 (MQ=255)
                                                       aaCAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATa                                >  1:3980750/1‑93 (MQ=255)
                                                         cAGTCTGCGGATCAGTCACCACCCCAGATTCGTGTTTTTCACAGGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaaa                              <  2:2968958/93‑1 (MQ=255)
                                                         cAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATaaa                              <  1:2326164/93‑1 (MQ=255)
                                                            tCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGcc                           <  1:3434609/93‑1 (MQ=255)
                                                            tCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGcc                           <  1:1320682/93‑1 (MQ=255)
                                                             cTGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCg                          <  1:1511114/93‑1 (MQ=255)
                                                              tGCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGa                         >  2:3917799/1‑93 (MQ=255)
                                                               tccgATCAGTCACCACACCAGATTCGTGTTTTTCTCAGTGCAAGCTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGAt                        <  2:3169890/90‑1 (MQ=255)
                                                               gCGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGAt                        >  1:3081602/1‑93 (MQ=255)
                                                                cGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                       <  2:2151439/93‑1 (MQ=255)
                                                                cGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                       <  1:2861947/93‑1 (MQ=255)
                                                                cGGATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATa                       <  2:2392394/93‑1 (MQ=255)
                                                                  gATCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATAcc                     <  2:756082/93‑1 (MQ=255)
                                                                   aTCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTCCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCg                    >  1:4021071/1‑93 (MQ=255)
                                                                    tCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGc                   >  2:229873/1‑93 (MQ=255)
                                                                    tCAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGc                   >  1:3394431/1‑93 (MQ=255)
                                                                     cAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGcc                  <  1:4290450/93‑1 (MQ=255)
                                                                     cAGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGCTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGcc                  >  2:3300730/1‑93 (MQ=255)
                                                                      aGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCt                 <  1:53971/93‑1 (MQ=255)
                                                                      aGTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCt                 <  2:2251878/93‑1 (MQ=255)
                                                                       gTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCtt                >  2:2460758/1‑93 (MQ=255)
                                                                       gTCACCACACCAGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCtt                >  2:1332415/1‑93 (MQ=255)
                                                                         caccacaccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTcc              <  1:668262/93‑1 (MQ=255)
                                                                          accacaccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCt             <  2:2137023/93‑1 (MQ=255)
                                                                          accacaccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCt             <  1:737886/93‑1 (MQ=255)
                                                                            cacaccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGa           >  2:341151/1‑93 (MQ=255)
                                                                             acaccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAg          >  1:1154775/1‑93 (MQ=255)
                                                                               accaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGcc        <  2:2767284/93‑1 (MQ=255)
                                                                               accaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGcc        <  2:4281770/93‑1 (MQ=255)
                                                                                ccaGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCa       >  2:2020241/1‑93 (MQ=255)
                                                                                 caGATTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCag      >  1:3909044/1‑93 (MQ=255)
                                                                                    aTTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGc   <  2:784306/93‑1 (MQ=255)
                                                                                    aTTCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGc   >  2:2588087/1‑93 (MQ=255)
                                                                                     ttCGTGTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCa  >  1:3710171/1‑93 (MQ=255)
                                                                                           |                                                                                       
CAAAACCGTTACGCTCGGTCAGTTCTTTCACTTCGCGCAGCGTCGTGGTTGGCAGAACAGTCTGCGGATCAGTCACCACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCA  >  NC_000913/2631736‑2631914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: