Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,522 C→T R42H (CGT→CAT)  fur ← DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,5220CT98.7% 580.7 / ‑2.3 159R42H (CGT→CAT) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (72/85);  total (72/87)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 159 total aligned reads displayed.

CAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGC  >  NC_000913/710427‑710616
                                                                                               |                                                                                              
cAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGttt                                                                                            >  1:3950833/1‑100 (MQ=255)
cAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGttt                                                                                            >  1:377698/1‑100 (MQ=255)
 aaaaTTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAGTTTCTTCACCCATATCGATCAGATGTTTg                                                                                           <  1:4748541/100‑1 (MQ=255)
   aaTTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGta                                                                                         <  1:1481342/100‑1 (MQ=255)
       gtgGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGGAGCCAGCCCAATTTCTTCACCCATATCGATCAGATGTTTGGATaaa                                                                                     >  1:4844019/1‑100 (MQ=255)
           cGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGttt                                                                                            <  2:4309251/89‑1 (MQ=255)
           cGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGttt                                                                                            >  2:4588141/1‑89 (MQ=255)
           cGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGttt                                                                                            >  2:4703475/1‑89 (MQ=255)
            gggTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTg                                                                                           <  2:2317420/89‑1 (MQ=255)
            gggTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTc                                                                                <  1:2718548/100‑1 (MQ=255)
             ggTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGt                                                                                          <  2:3805687/89‑1 (MQ=255)
              gTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGta                                                                                         <  2:92592/89‑1 (MQ=255)
              gTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTCCACCCATATCGATCAGATGTTTGTATAAATCTTCCg                                                                              <  1:133711/100‑1 (MQ=255)
               tGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGc                                                                             <  1:4055395/100‑1 (MQ=255)
                  cGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATaaa                                                                                     >  2:3154092/1‑89 (MQ=255)
                  cGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACt                                                                          >  1:1483986/1‑100 (MQ=255)
                     tACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGAc                                                                       >  1:2935613/1‑100 (MQ=255)
                      aCCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCtt                                                                                 >  2:1816197/1‑89 (MQ=255)
                       ccAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGt                                                                     >  1:3460407/1‑100 (MQ=255)
                         aGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCg                                                                              <  2:1383341/89‑1 (MQ=255)
                         aGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGa                                                                   >  1:170612/1‑100 (MQ=255)
                          gCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGc                                                                             >  2:4586888/1‑89 (MQ=255)
                          gCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGAt                                                                  >  1:2170327/1‑100 (MQ=255)
                            gtcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATgg                                                                <  1:321513/100‑1 (MQ=255)
                             tcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGt                                                               <  1:1643298/100‑1 (MQ=255)
                             tcgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGt                                                               <  1:190339/100‑1 (MQ=255)
                              cgtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTg                                                                         <  2:304094/89‑1 (MQ=255)
                               gtcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGa                                                                        >  2:582991/1‑89 (MQ=255)
                                  aaaCTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGt                                                                     >  2:3079603/1‑89 (MQ=255)
                                       ggTTCAGTACGCGATATACCGTAGCTAGACCAATAACAAAACAAAAAAAGAAAAGATGTTTGTATAAATCTTCCGCACTGACGTGATgg                                                                <  2:2107416/89‑1 (MQ=255)
                                       ggTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATgg                                                                >  2:1539277/1‑89 (MQ=255)
                                        gTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGt                                                               >  2:2868037/1‑89 (MQ=255)
                                         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCt                                                   >  1:2789844/1‑100 (MQ=255)
                                          tCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTg                                                  <  1:742755/100‑1 (MQ=255)
                                           cAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTga                                                 >  1:245047/1‑100 (MQ=255)
                                            aGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTc                                                           >  2:518576/1‑89 (MQ=255)
                                               aCGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                             >  1:4559194/1‑100 (MQ=255)
                                               aCGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                             >  1:95598/1‑100 (MQ=255)
                                                cgcgATATACCGTAGCTAGACCAATTTTTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGc                                                       <  2:314278/89‑1 (MQ=255)
                                                 gcgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCt                                                      <  2:4708620/89‑1 (MQ=255)
                                                 gcgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACt                                           <  1:4767285/100‑1 (MQ=255)
                                                 gcgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACt                                           >  1:4006368/1‑100 (MQ=255)
                                                  cgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACtt                                          >  1:900270/1‑100 (MQ=255)
                                                   gATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTcc                                                    <  2:3954262/89‑1 (MQ=255)
                                                     tataCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTg                                                  <  2:4158229/89‑1 (MQ=255)
                                                      ataCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAg                                      <  1:4227174/100‑1 (MQ=255)
                                                        aCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                               <  2:2832802/89‑1 (MQ=255)
                                                        aCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                               <  2:1904230/89‑1 (MQ=255)
                                                        aCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaag                                               <  2:674021/89‑1 (MQ=255)
                                                          cGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                             >  2:92700/1‑89 (MQ=255)
                                                          cGTAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa                                             >  2:4858352/1‑89 (MQ=255)
                                                            tAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACt                                           >  2:2236869/1‑89 (MQ=255)
                                                            tAGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAtttt                                >  1:626854/1‑100 (MQ=255)
                                                             aGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACtt                                          <  2:1753560/89‑1 (MQ=255)
                                                             aGCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAttttt                               >  1:3019725/1‑100 (MQ=255)
                                                              gCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTc                                         >  2:3426778/1‑89 (MQ=255)
                                                              gCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTa                              <  1:3756786/100‑1 (MQ=255)
                                                              gCCAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTa                              <  1:1967426/100‑1 (MQ=255)
                                                               ccAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTcc                                        >  2:4524280/1‑89 (MQ=255)
                                                                cAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCa                                       >  2:4534883/1‑89 (MQ=255)
                                                                cAGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCa                                       >  2:2140839/1‑89 (MQ=255)
                                                                cAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCa                                       <  2:1492077/89‑1 (MQ=255)
                                                                 aGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAg                                      <  2:1338317/89‑1 (MQ=255)
                                                                 aGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAg                                      <  2:3398032/89‑1 (MQ=255)
                                                                 aGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAg                                      <  2:1083607/89‑1 (MQ=255)
                                                                 aGACCAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAAc                           <  1:1044329/100‑1 (MQ=255)
                                                                    ccAATTTCTTCACCCATATCGATCAGATGTTTGTATAAGTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAg                        >  1:4862372/1‑100 (MQ=255)
                                                                    ccAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAt                                   >  2:2725961/1‑89 (MQ=255)
                                                                    ccAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAg                        >  1:2432534/1‑100 (MQ=255)
                                                                    ccAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAg                        >  1:509414/1‑100 (MQ=255)
                                                                     cAATTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAtt                                  >  2:1966804/1‑89 (MQ=255)
                                                                       aTTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAtttt                                <  2:1066081/89‑1 (MQ=255)
                                                                       aTTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGAtttt                                <  2:3906015/89‑1 (MQ=255)
                                                                       aTTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGaa                     >  1:2167/1‑100 (MQ=255)
                                                                       aTTTCTTCACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGaa                     >  1:4056841/1‑100 (MQ=255)
                                                                         ttcttcACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTa                              <  2:1193513/89‑1 (MQ=255)
                                                                            ttcACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAAc                           <  2:2531794/89‑1 (MQ=255)
                                                                            ttcACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAAc                           <  2:1167011/89‑1 (MQ=255)
                                                                            ttcACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGtt                <  1:4435419/100‑1 (MQ=255)
                                                                             tcACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACg                          >  2:2884121/1‑89 (MQ=255)
                                                                              cACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGa                         >  2:953060/1‑89 (MQ=255)
                                                                              cACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTAc              <  1:1042333/100‑1 (MQ=255)
                                                                              cACCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTAc              >  1:1328769/1‑100 (MQ=255)
                                                                               aCCCATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACt             <  1:3202170/100‑1 (MQ=255)
                                                                                cccATATCGATCAGATGTTTGTATAAATCTTCTGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGACGCGTTACtt            <  1:2396145/100‑1 (MQ=255)
                                                                                cccATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAgg                       <  2:1665318/89‑1 (MQ=255)
                                                                                cccATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACtt            >  1:4222302/1‑100 (MQ=255)
                                                                                cccATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACtt            <  1:649603/100‑1 (MQ=255)
                                                                                 ccATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGa                      <  2:4150704/89‑1 (MQ=255)
                                                                                 ccATATCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGa                      <  2:2749574/89‑1 (MQ=255)
                                                                                   atatCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAg                    >  2:2897552/1‑89 (MQ=255)
                                                                                    tatCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGc                   <  2:673112/89‑1 (MQ=255)
                                                                                     atCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCg                  >  2:4790308/1‑89 (MQ=255)
                                                                                     atCGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAgg       <  1:1813547/100‑1 (MQ=255)
                                                                                       cGATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGtt                <  2:3898013/89‑1 (MQ=255)
                                                                                        gATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCa    >  1:2099805/1‑100 (MQ=255)
                                                                                        gATCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCa    >  1:4708979/1‑100 (MQ=255)
                                                                                          tCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACt             <  2:1614169/89‑1 (MQ=255)
                                                                                          tCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACt             >  2:3948459/1‑89 (MQ=255)
                                                                                          tCAGATGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGc  <  1:3999908/100‑1 (MQ=255)
                                                                                               |                                                                                              
CAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGC  >  NC_000913/710427‑710616

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: