Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,098,512 C→A R259R (CGC→CGA rhaS → transcriptional activator of rhaBAD and rhaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,098,5120CA97.7% 293.8 / ‑0.5 85R259R (CGC→CGArhaStranscriptional activator of rhaBAD and rhaT
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (44/39);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATAT  >  NC_000913/4098415‑4098605
                                                                                                 |                                                                                             
aaGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACg                                                                                             <  1:394824/100‑1 (MQ=255)
 aGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTcgacga                                                                                            >  1:1097559/1‑100 (MQ=255)
  gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                           >  2:286351/1‑100 (MQ=255)
  gCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgag                                                                                           >  2:681631/1‑100 (MQ=255)
   cccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACgaga                                                                                          >  1:1550328/1‑100 (MQ=255)
    ccGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGagag                                                                                         >  2:1160343/1‑100 (MQ=255)
         aTCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTaa                                                                                    >  2:144325/1‑100 (MQ=255)
          tCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAAc                                                                                   >  2:1565708/1‑100 (MQ=255)
             gCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTgg                                                                                >  2:1724127/1‑100 (MQ=255)
              cTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGt                                                                               >  1:1517350/1‑100 (MQ=255)
              cTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGt                                                                               >  2:757425/1‑100 (MQ=255)
              cTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGt                                                                               >  1:106832/1‑100 (MQ=255)
               tACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTc                                                                              <  1:1242024/100‑1 (MQ=255)
               tACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTc                                                                              <  2:1141262/100‑1 (MQ=255)
               tACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTc                                                                              >  1:986335/1‑100 (MQ=255)
                aCGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCa                                                                             >  1:997852/1‑100 (MQ=255)
                 cGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCAc                                                                            >  2:32152/1‑100 (MQ=255)
                  gCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCAcc                                                                           >  2:1144590/1‑100 (MQ=255)
                   ccACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGCGTTTAACTGGTCACcg                                                                          >  1:714464/1‑100 (MQ=255)
                   ccACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcg                                                                          <  2:5785/100‑1 (MQ=255)
                    cacaGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACcgc                                                                         >  1:374226/1‑100 (MQ=255)
                       aGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGtn                                                                      <  2:523151/100‑2 (MQ=255)
                       aGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                      <  1:2982/100‑1 (MQ=255)
                       aGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTg                                                                      <  2:1326563/100‑1 (MQ=255)
                          gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGata                                                                   <  2:1453365/100‑1 (MQ=255)
                          gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGata                                                                   <  2:593685/100‑1 (MQ=255)
                          gAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGata                                                                   <  1:848452/100‑1 (MQ=255)
                            ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                 >  1:1161561/1‑100 (MQ=255)
                            ggCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATAtt                                                                 >  1:1444542/1‑100 (MQ=255)
                              ccAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCg                                                               >  2:435824/1‑100 (MQ=255)
                               cAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGc                                                              >  2:1690595/1‑100 (MQ=255)
                                aGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                             <  2:1833176/100‑1 (MQ=255)
                                 gCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGcc                                                             <  1:1945306/99‑1 (MQ=255)
                                   gTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAgg                                                          <  2:1191181/100‑1 (MQ=255)
                                    ttACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAgg                                                          <  2:1097560/99‑1 (MQ=255)
                                     tACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAggg                                                         >  2:965895/1‑99 (MQ=255)
                                       cTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACg                                                      >  1:533317/1‑100 (MQ=255)
                                       cTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACg                                                      >  2:1570763/1‑100 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                    <  2:1780283/100‑1 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                    <  1:1549986/100‑1 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                    >  2:1491337/1‑100 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                    >  2:1248179/1‑100 (MQ=255)
                                         gACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                    >  2:731906/1‑100 (MQ=255)
                                           cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                  >  1:603363/1‑100 (MQ=255)
                                           cATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGAt                                                  >  1:475579/1‑100 (MQ=255)
                                            aTCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                 <  1:1469900/100‑1 (MQ=255)
                                            aTCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                 <  2:975560/100‑1 (MQ=255)
                                             tCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATgg                                                <  1:1807779/100‑1 (MQ=255)
                                              cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                               >  1:301721/1‑100 (MQ=255)
                                              cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                               >  2:507264/1‑100 (MQ=255)
                                              cGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGc                                               <  2:1284616/100‑1 (MQ=255)
                                               gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCt                                              >  2:1693673/1‑100 (MQ=255)
                                               gCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGCGATGGCt                                              >  2:87195/1‑100 (MQ=255)
                                                 ntATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCttt                                            >  1:1369721/2‑100 (MQ=255)
                                                    tCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                         <  2:1595492/100‑1 (MQ=255)
                                                    tCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                         <  2:450575/100‑1 (MQ=255)
                                                     cGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                        <  2:374228/100‑1 (MQ=255)
                                                      gCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCa                                       <  1:1724116/100‑1 (MQ=255)
                                                      gCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCa                                       >  1:705086/1‑100 (MQ=255)
                                                         gtgGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATa                                    >  2:1891576/1‑100 (MQ=255)
                                                              ttCAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAAcgcg                               <  1:805561/100‑1 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                             <  2:106832/100‑1 (MQ=255)
                                                                cAGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGGTTAACTGGTCACCGCGTGATATTCGCCAAGGACGGGATGGCTTTATACAATAACGCGaa                             >  1:847597/1‑100 (MQ=255)
                                                                atctgACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                             <  2:1011934/96‑1 (MQ=255)
                                                                 aGCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAt                            >  1:108341/1‑100 (MQ=255)
                                                                  gCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                           <  2:1161563/100‑1 (MQ=255)
                                                                  gCGACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtc                           <  2:1590481/100‑1 (MQ=255)
                                                                    gACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaaagat                         >  1:1011933/1‑96 (MQ=255)
                                                                    gACAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtctt                         >  2:12165/1‑100 (MQ=255)
                                                                     aCAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttc                        <  2:1550338/100‑1 (MQ=255)
                                                                     aCAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttc                        <  1:1147287/100‑1 (MQ=255)
                                                                     aCAGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtcttc                        >  1:1101581/1‑100 (MQ=255)
                                                                       aGTAACCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTtctc                      <  2:533319/100‑1 (MQ=255)
                                                                          aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                   >  2:944901/1‑100 (MQ=255)
                                                                          aaCCACTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAc                   <  2:986336/100‑1 (MQ=255)
                                                                              aCTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAt               >  2:1558132/1‑100 (MQ=255)
                                                                               cTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt              <  1:681630/100‑1 (MQ=255)
                                                                               cTTTTCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt              <  1:1565697/100‑1 (MQ=255)
                                                                                ttttCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTAtt              <  1:688652/99‑1 (MQ=255)
                                                                                 tttCGACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGACATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTg            <  1:1248175/100‑1 (MQ=255)
                                                                                     gACGCTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACg        >  2:1286904/1‑100 (MQ=255)
                                                                                         cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat    >  2:1274489/1‑100 (MQ=255)
                                                                                         cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat    <  2:1232866/100‑1 (MQ=255)
                                                                                         cTTTTTCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat    <  1:1144589/100‑1 (MQ=255)
                                                                                           ttttCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCatat  >  2:980955/1‑100 (MQ=255)
                                                                                           ttttCGACGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCat    <  1:106413/98‑1 (MQ=255)
                                                                                                 |                                                                                             
AAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATAT  >  NC_000913/4098415‑4098605

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