Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,510,121 +G intergenic (‑358/+151) yhfZ ← / ← trpS hypothetical protein/tryptophan tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,1201.G96.0% 708.8 / 17.2 202intergenic (‑357/+152)yhfZ/trpShypothetical protein/tryptophan tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (7/1);  new base G (120/74);  total (127/75)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGG  >  NC_000913/3510028‑3510209
                                                                                              |                                                                                         
aaGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTcccc                                                                                            <  1:3325264/93‑1 (MQ=255)
 aGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCg                                                                                          >  2:351935/1‑93 (MQ=255)
 aGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCg                                                                                          >  2:1474470/1‑93 (MQ=255)
    tAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtt                                                                                       >  2:3654180/1‑93 (MQ=255)
     aaGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgttt                                                                                      >  1:3275555/1‑92 (MQ=255)
      aGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttc                                                                                     <  1:5031341/93‑3 (MQ=255)
      aGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttc                                                                                     >  1:346288/1‑91 (MQ=255)
         cgcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga                                                                                   <  1:3090971/93‑1 (MQ=255)
          gcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag                                                                                  <  1:1468798/93‑1 (MQ=255)
          gcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag                                                                                  <  2:3305831/93‑1 (MQ=255)
               aGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatg                                                                             >  1:3372686/1‑93 (MQ=255)
                gCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatga                                                                            >  2:715939/1‑93 (MQ=255)
                gCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatga                                                                            <  2:1251052/93‑1 (MQ=255)
                gCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatga                                                                            >  1:4872240/1‑93 (MQ=255)
                  gagaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATa                                                                          <  2:1045332/93‑1 (MQ=255)
                   agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGATTCCCCGGTTTCAGAGTGATGATaa                                                                         >  2:2394514/1‑93 (MQ=255)
                    gaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaa                                                                        >  1:586266/1‑93 (MQ=255)
                     aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGTTaaaa                                                                       <  2:4618556/93‑1 (MQ=255)
                      tgtgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAg                                                                      >  1:3121360/1‑93 (MQ=255)
                       gtgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGc                                                                     <  2:6015809/93‑1 (MQ=255)
                          gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaa                                                                  >  1:5162941/1‑93 (MQ=255)
                          gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaa                                                                  >  1:365515/1‑93 (MQ=255)
                           aaaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAATTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaaa                                                                 >  2:3136155/1‑93 (MQ=255)
                           aaaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaaaa                                                                 <  2:3359099/93‑1 (MQ=255)
                            aaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAt                                                                >  2:5002756/1‑93 (MQ=255)
                             aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAtt                                                               <  2:2895237/93‑1 (MQ=255)
                               ggCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGc                                                             <  2:5092354/93‑1 (MQ=255)
                               ggCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGc                                                             >  1:2025500/1‑93 (MQ=255)
                                gCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGcc                                                            <  2:3881497/93‑1 (MQ=255)
                                   gggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         <  2:4040490/92‑1 (MQ=255)
                                   aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         <  1:4230599/93‑1 (MQ=255)
                                   aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         >  2:5304497/1‑93 (MQ=255)
                                   aGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGa                                                         <  1:2793612/93‑1 (MQ=255)
                                    ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                                        >  2:5443668/1‑93 (MQ=255)
                                    ggCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGAt                                                        <  2:4434542/93‑1 (MQ=255)
                                       aCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg                                                     >  1:708985/1‑93 (MQ=255)
                                       aCTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcg                                                     >  2:5898575/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    >  1:5415574/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                    >  1:128133/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTGCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                   >  2:3408402/1‑93 (MQ=255)
                                          tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTa                                                  >  2:516384/1‑93 (MQ=255)
                                           cacaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                                                 >  2:1356576/1‑93 (MQ=255)
                                                tttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                                            <  2:4956942/93‑1 (MQ=255)
                                                tttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTa                                            <  2:476490/93‑1 (MQ=255)
                                                  tttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                                          >  2:5830468/1‑93 (MQ=255)
                                                   ttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCa                                         >  2:254812/1‑93 (MQ=255)
                                                     cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                       >  2:5591095/1‑93 (MQ=255)
                                                     cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                       <  2:561596/93‑1 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                      >  1:2488634/1‑93 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                      <  1:6042108/93‑1 (MQ=255)
                                                       gATTATTGCAAATGAGGATAAAAAACCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGcc                                     <  2:119834/93‑1 (MQ=255)
                                                         ttattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                   >  2:1733333/1‑93 (MQ=255)
                                                          tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                  >  2:2992739/1‑93 (MQ=255)
                                                            ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt                                <  2:2719405/93‑1 (MQ=255)
                                                            ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCATACAt                                >  2:1546501/1‑93 (MQ=255)
                                                             tGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt                               <  2:1705158/93‑1 (MQ=255)
                                                               cAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTc                             <  2:1426699/93‑1 (MQ=255)
                                                                aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc                            >  1:616592/1‑93 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           <  2:4901315/93‑1 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           >  2:2124833/1‑93 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           >  1:2949502/1‑93 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           >  1:3960102/1‑93 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                           <  1:3837072/93‑1 (MQ=255)
                                                                   tGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTg                         <  1:2854796/93‑1 (MQ=255)
                                                                      ggATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCaa                      <  2:5718615/93‑1 (MQ=255)
                                                                       gATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAat                     >  1:5203175/1‑93 (MQ=255)
                                                                       gATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAat                     <  2:5271821/93‑1 (MQ=255)
                                                                        aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAata                    <  1:4057218/93‑1 (MQ=255)
                                                                         tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAatat                   >  2:2070289/1‑93 (MQ=255)
                                                                         tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAatat                   >  1:5307761/1‑93 (MQ=255)
                                                                          aaaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtg                  <  2:3946140/93‑1 (MQ=255)
                                                                          aaaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtg                  <  2:1528754/93‑1 (MQ=255)
                                                                          aaaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtg                  >  1:2161170/1‑93 (MQ=255)
                                                                           aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtgt                 >  1:2686405/1‑93 (MQ=255)
                                                                            aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtgtg                >  1:235064/1‑93 (MQ=255)
                                                                            aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtgtg                >  2:3975571/1‑93 (MQ=255)
                                                                            aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtgtg                <  1:4416681/93‑1 (MQ=255)
                                                                             aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGc               <  1:1066660/93‑1 (MQ=255)
                                                                             aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGc               <  2:5011230/93‑1 (MQ=255)
                                                                             aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGc               >  2:3018918/1‑93 (MQ=255)
                                                                              aaaCCGGGTTTCCCCGGGTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCa              >  2:1686054/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt             >  1:1020927/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt             >  1:3028572/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt             <  1:37908/93‑1 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCGTATCAGGCCTACATTTCCTTGCAATATGTGCAt             >  2:5930783/1‑93 (MQ=255)
                                                                                aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt            >  1:467937/1‑93 (MQ=255)
                                                                                aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt            >  1:2572438/1‑93 (MQ=255)
                                                                                aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt            <  1:2092829/93‑1 (MQ=255)
                                                                                 ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTa           >  1:804204/1‑93 (MQ=255)
                                                                                   gggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACt         >  1:2258343/1‑93 (MQ=255)
                                                                                     gTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACttt       >  1:561648/1‑93 (MQ=255)
                                                                                     gTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACttt       <  2:5634508/93‑1 (MQ=255)
                                                                                      tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTg      <  2:2563592/93‑1 (MQ=255)
                                                                                       ttccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGt     <  2:1795158/90‑1 (MQ=255)
                                                                                       ttccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGt     <  2:2643919/90‑1 (MQ=255)
                                                                                        tccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTa    <  1:3563104/91‑1 (MQ=255)
                                                                                        tCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTa    >  2:4311815/1‑93 (MQ=255)
                                                                                          cccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAgg  >  1:462894/1‑93 (MQ=255)
                                                                                          cccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAgg  >  2:6117581/1‑93 (MQ=255)
                                                                                          cccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAgg  >  2:3517842/1‑93 (MQ=255)
                                                                                              |                                                                                         
AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGG  >  NC_000913/3510028‑3510209

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