Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,595 G→T P18T (CCT→ACT)  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,5950GT98.4% 189.3 / NA 62P18T (CCT→ACT) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (36/25);  minor base A (1/0);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAT  >  NC_000913/710459‑710733
                                                                                                                                        |                                                                                                                                          
tcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGt                                                                                                                                      >  1:712049/1‑143 (MQ=255)
tcAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGt                                                                                                                                      >  2:572579/1‑143 (MQ=255)
 cAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGCTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGtt                                                                                                                                     >  2:1063535/1‑143 (MQ=255)
 cAAACTCGTTCATCACGCGATATACCGTAGCCAGACCAATTTCTTCCCCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGtt                                                                                                                                     <  2:1113860/143‑1 (MQ=255)
  aaaCTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGATTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACATCAAGGATTTTTAAAGGAGTAAGCGGTa                                                                                                                                    >  2:439832/1‑143 (MQ=255)
      tGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACttt                                                                                                                                >  1:830710/1‑143 (MQ=255)
      tGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACttt                                                                                                                                >  1:313793/1‑143 (MQ=255)
       tgTTCAGCATTCGATCTACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAACCTTCCGCACTTACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTc                                                                                                                               <  1:1119372/142‑1 (MQ=255)
         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGCTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTCTTACACGAGAACGCGCTACTTTCAg                                                                                                                             >  2:914290/1‑143 (MQ=255)
         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTAAg                                                                                                                             >  1:695618/1‑143 (MQ=255)
         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTCTCGTGATGATTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTAAg                                                                                                                             >  2:933126/1‑143 (MQ=255)
         ttCAGTACGCGAGATACCGTAGCCACACCAATTTCTTCACCCATATCGAGCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAg                                                                                                                             <  1:454672/143‑1 (MQ=255)
          tCAGTACGCGATAAACCGTAACCAGACAAATTACTTCACCCATAACAATCAGACATATGTATAAATCTTCCGCACTGACGTGACAGTTGTCCGGATCATGAAGAAATTCCAGGATCTTTAATCGACTAACCGTTAcataccag                                                                                                                            >  2:158744/1‑136 (MQ=255)
              tACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAg                                                                                                                        <  2:702044/143‑1 (MQ=255)
                  cgATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGcttt                                                                                                                    >  1:404519/1‑143 (MQ=255)
                         ccGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCTAGATTTCTTTAgg                                                                                                             >  2:497686/1‑143 (MQ=255)
                               ccAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTa                                                                                                       >  1:277250/1‑143 (MQ=255)
                                      aaTTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTAt                                                                                                <  1:886138/143‑1 (MQ=255)
                                      aaTTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTAt                                                                                                <  1:232700/143‑1 (MQ=255)
                                      aaTTTCTGCACCCATATCGATGAGGCGTTTGTATAAGTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTAt                                                                                                <  1:655455/143‑1 (MQ=255)
                                        tttCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCa                                                                                              <  1:572579/143‑1 (MQ=255)
                                         ttcttcACCCATATCGATCAGACTTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCTGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAg                                                                                             <  2:716574/143‑1 (MQ=255)
                                         ttcttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAg                                                                                             <  2:717202/143‑1 (MQ=255)
                                            ttcACCCATAACGAACTGATGTCTGTATACATCGTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCa                                                                                          <  2:861675/143‑1 (MQ=255)
                                                 ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCgg                                                                                     >  1:340443/1‑143 (MQ=255)
                                                         aTCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTc                                                                             >  1:641991/1‑143 (MQ=255)
                                                          tCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTcc                                                                            >  1:985314/1‑143 (MQ=255)
                                                          tCAGACGGTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGCACTTCCAGGATTTTTAAAGGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTGATCAGTCATGCGGAATCTGTcc                                                                            >  2:864629/1‑143 (MQ=255)
                                                                    gTAGAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCGAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTATGGCGGTATTGTTATCAGTCATGGGGAATCTGTTCTGTTACTAAg                                                                  <  2:1078732/143‑1 (MQ=255)
                                                                        aaaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGt                                                              >  1:638849/1‑143 (MQ=255)
                                                                           tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCa                                                           <  1:1063535/143‑1 (MQ=255)
                                                                              tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGACCTACCAAGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACtt                                                        >  2:945208/1‑143 (MQ=255)
                                                                              tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTTAGTGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACtt                                                        <  1:252729/143‑1 (MQ=255)
                                                                              tCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACtt                                                        <  1:522982/143‑1 (MQ=255)
                                                                                cGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCa                                                      <  2:277250/143‑1 (MQ=255)
                                                                                 gCACTGACGTGATGGTTGTCCGGCTCCTGAAGTACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAt                                                     <  1:1120368/143‑1 (MQ=255)
                                                                                  cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                    >  1:1016201/1‑143 (MQ=255)
                                                                                   aCTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACt                                                         <  1:777811/137‑1 (MQ=255)
                                                                                   aCTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACt                                                         >  2:777811/1‑137 (MQ=255)
                                                                                     tGACGTGATGGTTGTCCGGCTCCTGAAGAACTGCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTaga                                                 <  1:741622/143‑1 (MQ=255)
                                                                                      gACGTAATTTTTGTCCGTATACTGGCGAACTGCCATGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaa                                                <  1:96308/143‑1 (MQ=255)
                                                                                      ctcGTGCTGGTGGTCCGGCTCCTGAAGAACTGCCAGGATTTTTAAACGAGTCAGCGTTACTTTCAGGCCAGCTTTCTCTAGCGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTagaa                                                <  2:581612/141‑1 (MQ=255)
                                                                                         gTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAgaagaa                                             >  2:1155469/1‑143 (MQ=255)
                                                                                           gATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGCTATCAGTCATGCGGAATCTGTCCTGTTACTCAGCGGTTCACTTCATTAGAAGAAGt                                           >  2:254825/1‑143 (MQ=255)
                                                                                             tGGTTGTCCGGCTCCTGAAGGGCCTCCCGGATGTTTGAACGAGTGAGCGTTACTTTCCGGCCAGCTTTCTTTAGGGCGGTATTTTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAGTAGAAGAAGTgc                                         <  2:906628/143‑2 (MQ=255)
                                                                                                          ccTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGAAAGAACGTGctct                            >  2:60/1‑143 (MQ=255)
                                                                                                            tGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTg                          >  2:1024025/1‑143 (MQ=255)
                                                                                                             gaagaaCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCCCTTCATTAGAAGAAGTGACAGAATTTGCTCTTga                         >  1:869598/1‑143 (MQ=255)
                                                                                                             gaagaaCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTga                         >  1:64049/1‑143 (MQ=255)
                                                                                                             gaagaaCTGCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTga                         <  2:165605/143‑1 (MQ=255)
                                                                                                                 aaCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTCGACGAAGTGACGGAATTTGCTCTTGAGATa                     >  1:553653/1‑143 (MQ=255)
                                                                                                                    ttCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGCCAGAATCTGCTCTTGAGATAAta                  >  2:999312/1‑142 (MQ=255)
                                                                                                                    ttCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATg                  >  1:1003905/1‑143 (MQ=255)
                                                                                                                       cAGGATTTTTATACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGt               <  2:595379/143‑1 (MQ=255)
                                                                                                                       cAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCgc               >  2:342904/1‑142 (MQ=255)
                                                                                                                       cAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGt               >  1:1137089/1‑143 (MQ=255)
                                                                                                                        aGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGATCTTGAGGTAATGCGTa              >  1:696919/1‑143 (MQ=255)
                                                                                                                         ggATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTAATTAGAAGAAGTGACAGAATTAGCTCTTGAGATAATGCGAAt             >  2:1065956/1‑143 (MQ=255)
                                                                                                                           aTTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCTGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCa           <  2:833868/143‑1 (MQ=255)
                                                                                                                             ttttAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATTTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCAtt         >  2:52359/1‑143 (MQ=255)
                                                                                                                               ttAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtat       >  1:492481/1‑143 (MQ=255)
                                                                                                                                    cGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAt  >  1:1087629/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAT  >  NC_000913/710459‑710733

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