Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 661,369 T→C 63.5% N4D (AAT→GAT)  ybeF ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913661,3690TC63.5% 120.7 / 119.5 132N4D (AAT→GAT) ybeFLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base C (46/37);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 142 total aligned reads displayed.

AAATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCTCAC  >  NC_000913/661270‑661467
                                                                                                   |                                                                                                  
aaaTAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTGGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGattatt                                                                                                    <  2:3356925/100‑1 (MQ=255)
  aTAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCAc                                                                                                  <  1:2722611/100‑1 (MQ=255)
   tAGTCAGAAGGTTAAGATCAATATTTCGTAATGTTGTTAATATTTGTCGTTTATCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACt                                                                                                 <  1:2291596/100‑1 (MQ=255)
   tAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACt                                                                                                 <  1:2977432/100‑1 (MQ=255)
    aGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTa                                                                                                <  2:2866129/100‑1 (MQ=255)
     gTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTAt                                                                                               >  1:3206284/1‑100 (MQ=255)
     gTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACATTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTAt                                                                                               >  1:1396285/1‑100 (MQ=255)
     gTCAGAAGGTTAAGATCAATATGTCGTAAGGTTGTAGATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTAt                                                                                               >  1:625500/1‑100 (MQ=255)
       cAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACGTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATcc                                                                                             >  1:1086073/1‑100 (MQ=255)
       cAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATcc                                                                                             <  1:1612102/100‑1 (MQ=255)
         gAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCAc                                                                                           >  2:1280039/1‑100 (MQ=255)
          aaGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACg                                                                                          <  1:538064/100‑1 (MQ=255)
           aGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGc                                                                                         <  1:2373874/100‑1 (MQ=255)
            ggTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCa                                                                                        >  1:684604/1‑100 (MQ=255)
            ggTTAAAATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCa                                                                                        <  2:142464/100‑1 (MQ=255)
             gTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCAc                                                                                       >  2:3197872/1‑100 (MQ=255)
               tAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACtt                                                                                     >  2:2632281/1‑100 (MQ=255)
               tAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACtt                                                                                     <  1:3045333/100‑1 (MQ=255)
                 aGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTAc                                                                                   <  1:3164449/100‑1 (MQ=255)
                 aGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTAc                                                                                   >  1:2768224/1‑100 (MQ=255)
                  ggtCAATGTTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTAAATTTGATTATCACTATGCACGCACTTa                                                                                    >  1:757540/3‑98 (MQ=255)
                  gATCAATATTTCGTAAGGTTGTAAATATTTTTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGAGTTCAATTTGATTATCACTATCCACGCACt                                                                                      >  1:2068764/1‑96 (MQ=255)
                  gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                                  <  1:9657/100‑1 (MQ=255)
                   aTCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTTACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTc                                                                                 <  1:347627/100‑1 (MQ=255)
                   aTCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTc                                                                                 <  1:1600585/100‑1 (MQ=255)
                    tCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTcc                                                                                >  2:710076/1‑100 (MQ=255)
                       atatTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAt                                                                             >  1:1012338/1‑100 (MQ=255)
                         atTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAttt                                                                           <  2:66545/100‑1 (MQ=255)
                          tttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAtttt                                                                          <  2:294400/100‑1 (MQ=255)
                              gTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTc                                                                      <  2:2488734/100‑1 (MQ=255)
                                 aGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATg                                                                   <  1:1893640/100‑1 (MQ=255)
                                  ggTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATgg                                                                  >  2:65084/1‑100 (MQ=255)
                                  ggTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATgg                                                                  <  1:1049646/100‑1 (MQ=255)
                                   gTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTAATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGa                                                                 >  1:392689/1‑100 (MQ=255)
                                    ttGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGaa                                                                <  1:157509/100‑1 (MQ=255)
                                     tGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGaaa                                                               >  1:1217503/1‑100 (MQ=255)
                                      gTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGaa                                                                <  1:2840238/98‑1 (MQ=255)
                                        aaaTATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaat                                                            <  2:1644830/100‑1 (MQ=255)
                                        aaaTATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaat                                                            <  2:3309978/100‑1 (MQ=255)
                                        aaaTATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaat                                                            <  2:1937835/100‑1 (MQ=255)
                                          atatTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaataa                                                          <  2:1750238/100‑1 (MQ=255)
                                            atTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAataata                                                        <  1:2344572/100‑1 (MQ=255)
                                            atTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAataata                                                        <  2:2568873/100‑1 (MQ=255)
                                              ttGTGGTTTACCTTCTGTCGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAtt                                                      <  1:3348166/100‑1 (MQ=255)
                                              ttGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAtt                                                      <  1:2386445/100‑1 (MQ=255)
                                               tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                     >  2:552119/1‑100 (MQ=255)
                                                gtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTa                                                    >  2:2320683/1‑100 (MQ=255)
                                                gtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAGTTTATTCATGGAAAAATAATATTTa                                                    <  1:218118/100‑1 (MQ=255)
                                                      tACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                              >  2:930091/1‑100 (MQ=255)
                                                      tACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                              >  1:609774/1‑100 (MQ=255)
                                                      tACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                              <  2:718915/100‑1 (MQ=255)
                                                         cTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAc                                           <  2:3191807/100‑1 (MQ=255)
                                                         cTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAc                                           >  2:117414/1‑100 (MQ=255)
                                                          ttCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACa                                          <  1:1360286/100‑1 (MQ=255)
                                                          ttCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACa                                          >  2:1828479/1‑100 (MQ=255)
                                                          ttCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACa                                          >  2:170576/1‑100 (MQ=255)
                                                           tCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaa                                         >  2:2235498/1‑100 (MQ=255)
                                                           tCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaa                                         >  1:978118/1‑100 (MQ=255)
                                                                cgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCCCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATc                                    >  1:2901122/1‑100 (MQ=255)
                                                                 gacTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATACACGCACTTACTCCAATATTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCt                                   <  1:137525/100‑1 (MQ=255)
                                                                   cTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTa                                 >  2:2692698/1‑100 (MQ=255)
                                                                   cTTTCAGCTAAAACAGTCTTCAATTTGAATATCACTAAACACGCACTTAATCCAATTTCATTCATGGAAAAATAATATTTAAAAAATTACAATAATc                                    <  1:1450292/97‑1 (MQ=255)
                                                                    tttCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTAt                                <  2:1758527/100‑1 (MQ=255)
                                                                    tttCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTAt                                >  1:633962/1‑100 (MQ=255)
                                                                     ttCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATg                               <  2:325987/100‑1 (MQ=255)
                                                                       cTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTc                             <  1:710073/100‑1 (MQ=255)
                                                                        tGCTTAAACAGGGTTCAATTTGATTAGTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCt                            >  2:2056156/1‑100 (MQ=255)
                                                                         gCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                           >  2:2361007/1‑100 (MQ=255)
                                                                         gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGtata                           >  1:2092496/1‑100 (MQ=255)
                                                                         gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                           >  2:1765728/1‑100 (MQ=255)
                                                                         gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                           >  2:1229937/1‑100 (MQ=255)
                                                                          catAAACTAGGTTCAATTTGATTATCACTATAAAAGCACTTAATACAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATg                               <  1:1630762/93‑1 (MQ=255)
                                                                          cTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTaa                          >  2:1284790/1‑100 (MQ=255)
                                                                          cTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTaa                          >  1:2122769/1‑100 (MQ=255)
                                                                            tAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAAtt                        <  2:3031752/100‑1 (MQ=255)
                                                                             aaaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTg                       >  2:2453889/1‑100 (MQ=255)
                                                                             aaaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTg                       <  2:677878/100‑1 (MQ=255)
                                                                              aaCAGTGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                      <  1:288277/100‑1 (MQ=255)
                                                                              aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                      >  2:1895971/1‑100 (MQ=255)
                                                                                   ggTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaac                   <  1:2320666/98‑1 (MQ=255)
                                                                                   ggTTCAATTTGATTATCACTATCAACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacgg                 <  1:1644161/100‑1 (MQ=255)
                                                                                    gTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacgga                >  2:956909/1‑100 (MQ=255)
                                                                                    gTTCAATGTGATTATCACTATCCACGCACTTACTCCAATTTTATTCTTGGAAAAATAATATTTACACAATTACAATAATCTTATGTCTAATTggaac                   >  1:1666242/1‑97 (MQ=255)
                                                                                       cAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgaaacggaa               >  1:2391410/1‑98 (MQ=255)
                                                                                       cAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacggaac              >  1:2592355/1‑99 (MQ=255)
                                                                                        aaTTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGgaacggaacg             >  2:1888725/1‑99 (MQ=255)
                                                                                        aaTTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGc            >  2:1621406/1‑100 (MQ=255)
                                                                                            tGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtttt        >  2:2787437/1‑100 (MQ=255)
                                                                                             gATTCTCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTAAAAAAATTACAATAATCTTATATCTAATTggaacggaac              >  1:179507/1‑93 (MQ=255)
                                                                                             gATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATACAATTATGTATAATTGAAACGGAACACtttt        >  1:831460/1‑99 (MQ=255)
                                                                                             gATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTg       <  2:1959499/100‑1 (MQ=255)
                                                                                              attattACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc      <  2:2554661/100‑1 (MQ=255)
                                                                                              attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc      <  1:570848/100‑1 (MQ=255)
                                                                                              attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc      <  2:726980/100‑1 (MQ=255)
                                                                                              attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc      <  2:871910/100‑1 (MQ=255)
                                                                                              attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc      <  2:879860/100‑1 (MQ=255)
                                                                                                tatcaCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCTc    <  1:1130765/100‑1 (MQ=255)
                                                                                                  ttaCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAACATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTACTCAc  >  1:988144/1‑100 (MQ=255)
                                                                                                  ttaCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCTCAc  >  1:2546719/1‑100 (MQ=255)
                                                                                                  ttaCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCTCAc  <  2:2782545/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCTCAC  >  NC_000913/661270‑661467

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: