Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,663,466 G→T 7.4% V9L (GTG→TTG)  suhB → inositol monophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,663,4660GT 7.4% 382.9 / 17.5 108V9L (GTG→TTG) suhBinositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base T (3/5);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

Only 100 of 108 total aligned reads displayed.

AGCATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAA  >  NC_000913/2663319‑2663612
                                                                                                                                                   |                                                                                                                                                  
agcaTTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                   <  2:1404021/149‑1 (MQ=255)
agcaTTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                   <  2:1779918/149‑1 (MQ=255)
       tctttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                            >  3:1116925/1‑149 (MQ=255)
         tttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagc                                                                                                                                          >  3:1263951/1‑149 (MQ=255)
         tttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagc                                                                                                                                          <  3:205513/149‑1 (MQ=255)
                  cTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                 <  3:497208/149‑1 (MQ=255)
                   tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgca                                                                                                                                     >  2:1421226/1‑144 (MQ=255)
                   tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcagcgcgca                                                                                                                                     <  1:1421226/144‑1 (MQ=255)
                   tAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCg                                                                                                                                >  1:1019742/1‑149 (MQ=255)
                    aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               <  2:1019742/149‑1 (MQ=255)
                    aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               <  1:367942/149‑1 (MQ=255)
                    aaaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               <  4:1263951/149‑1 (MQ=255)
                     aaGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGCgat                                                                                                                              <  1:1724388/149‑3 (MQ=255)
                       gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgc                                                                                                                                                 <  1:1324906/128‑1 (MQ=255)
                              cACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                     <  4:643304/149‑1 (MQ=255)
                               aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                           >  2:495983/1‑126 (MQ=255)
                               aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgcag                                                                                                                                           <  1:495983/126‑1 (MQ=255)
                                       tATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                            <  2:34608/149‑1 (MQ=255)
                                           aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACgg                                                                                                                                  >  3:1596978/1‑123 (MQ=255)
                                           aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACgg                                                                                                                                  <  4:1596978/123‑1 (MQ=255)
                                              aaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                            >  3:1875334/1‑110 (MQ=255)
                                              aaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                            <  4:1875334/110‑1 (MQ=255)
                                              aaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                    <  4:1432455/134‑1 (MQ=255)
                                              aaaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAACACTATGaa                                                                                                     >  1:702314/1‑149 (MQ=255)
                                               aaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                         >  4:631341/1‑144 (MQ=255)
                                               aaCGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                         <  3:631341/144‑1 (MQ=255)
                                                  gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgc                                                                                                                                                 <  1:1251712/101‑1 (MQ=255)
                                                  gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTgc                                                                                                                                                 >  2:1251712/1‑101 (MQ=255)
                                                  gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                 >  1:1813849/1‑149 (MQ=255)
                                                  gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                 >  1:1621442/1‑149 (MQ=255)
                                                         gATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               <  3:919967/112‑1 (MQ=255)
                                                         gATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               >  4:919967/1‑112 (MQ=255)
                                                         gATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               <  1:1054618/112‑1 (MQ=255)
                                                         gATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCgg                                                                                                                               >  2:1054618/1‑112 (MQ=255)
                                                               gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGa                                                                                    >  3:873308/1‑149 (MQ=255)
                                                                aTCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGaa                                                                                   >  2:716257/1‑149 (MQ=255)
                                                                  ctcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGc                                                                                 <  4:1116925/149‑1 (MQ=255)
                                                                   tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                     <  1:1736933/112‑1 (MQ=255)
                                                                   tcGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                     >  2:1736933/1‑112 (MQ=255)
                                                                    cGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                               <  2:702314/149‑1 (MQ=255)
                                                                     gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAg                                                                                  >  4:1065212/1‑145 (MQ=255)
                                                                     gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAg                                                                                  <  3:1065212/145‑1 (MQ=255)
                                                                     gctgctATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                              >  2:815266/1‑149 (MQ=255)
                                                                          tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAg                                                                              >  3:979122/1‑144 (MQ=255)
                                                                          tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGa                                                                         >  4:1896208/1‑149 (MQ=255)
                                                                              ctctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                     >  1:7012/1‑149 (MQ=255)
                                                                               tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                    >  3:1316771/1‑149 (MQ=255)
                                                                                ctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                   <  3:1896208/149‑1 (MQ=255)
                                                                                    gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTaa                                                               <  2:1813849/149‑1 (MQ=255)
                                                                                        ttttCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAt                                                           >  2:1272685/1‑149 (MQ=255)
                                                                                             ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTc                                                        <  1:1080518/147‑1 (MQ=255)
                                                                                             ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTc                                                        >  2:1080518/1‑147 (MQ=255)
                                                                                             ccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGt                                                      <  2:1621442/149‑1 (MQ=255)
                                                                                              cGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg                                                     <  4:873308/149‑1 (MQ=255)
                                                                                              cGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg                                                     >  2:1213393/1‑149 (MQ=255)
                                                                                              cGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTg                                                     >  2:122609/1‑149 (MQ=255)
                                                                                               gTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                    >  3:1728731/1‑149 (MQ=255)
                                                                                               gTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                    >  3:256775/1‑149 (MQ=255)
                                                                                               gTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                    <  1:815266/149‑1 (MQ=255)
                                                                                                 tCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCaaaaa                                                                                                            <  2:723013/91‑1 (MQ=255)
                                                                                                   tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAc                                                   <  3:1498705/146‑1 (MQ=255)
                                                                                                   tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGAc                                                   >  4:1498705/1‑146 (MQ=255)
                                                                                                   tttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCaa                                                >  4:251558/1‑149 (MQ=255)
                                                                                                    ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTgcgcgctgcgcgcacgg                                                                                                                                  >  2:626856/1‑63 (MQ=37)
                                                                                                    ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACgg                                                                                                                                  <  1:626856/66‑1 (MQ=38)
                                                                                                    ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAAc                                               <  1:716257/149‑1 (MQ=255)
                                                                                                    ttAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAAc                                               >  3:299768/1‑149 (MQ=255)
                                                                                                       aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                                                                       <  1:1552650/106‑1 (MQ=255)
                                                                                                       aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGt                                                                                       >  2:1552650/1‑106 (MQ=255)
                                                                                                       aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                            <  4:1728731/149‑1 (MQ=255)
                                                                                                       aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                            <  4:1316771/149‑1 (MQ=255)
                                                                                                       aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                            <  4:256775/149‑1 (MQ=255)
                                                                                                               tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGa                                          <  2:1482199/143‑1 (MQ=255)
                                                                                                               tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                    >  2:1087281/1‑149 (MQ=255)
                                                                                                               tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                    <  1:1087281/149‑1 (MQ=255)
                                                                                                               tGAGAGAGACCGATGCATCCGATGCTAAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                    >  2:1449098/1‑149 (MQ=255)
                                                                                                                       aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                        >  1:1099283/1‑73 (MQ=255)
                                                                                                                       aCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                        <  2:1099283/73‑1 (MQ=255)
                                                                                                                         cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAg                                                                                  >  1:1067242/1‑93 (MQ=255)
                                                                                                                         cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAg                                                                                  <  2:1067242/93‑1 (MQ=255)
                                                                                                                           aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGa                                                                                             >  4:599208/1‑80 (MQ=255)
                                                                                                                           aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                        <  1:1272685/149‑1 (MQ=255)
                                                                                                                             gCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGa                      <  1:1449098/149‑1 (MQ=255)
                                                                                                                              cATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTg                                  <  1:595496/136‑1 (MQ=255)
                                                                                                                                tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                   >  3:1823634/1‑149 (MQ=255)
                                                                                                                                  cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                 <  4:1823634/149‑1 (MQ=255)
                                                                                                                                  cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                 >  2:182697/1‑149 (MQ=255)
                                                                                                                                   gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTg                       >  2:1152829/1‑142 (MQ=255)
                                                                                                                                   gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTg                       <  1:1152829/142‑1 (MQ=255)
                                                                                                                                   gATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCg                >  2:881783/1‑149 (MQ=255)
                                                                                                                                    aTGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGa               <  2:7012/149‑1 (MQ=255)
                                                                                                                                      gttcagcATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                 <  3:1892928/139‑1 (MQ=255)
                                                                                                                                      gttcagcATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                 >  4:1892928/7‑145 (MQ=255)
                                                                                                                                        tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACg           >  2:619670/1‑149 (MQ=255)
                                                                                                                                          aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGc                           >  4:645576/1‑131 (MQ=255)
                                                                                                                                          aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGc                           <  3:645576/131‑1 (MQ=255)
                                                                                                                                          aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAt         >  3:1217494/1‑149 (MQ=255)
                                                                                                                                                gCCTTGCGCGCTGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                   <  3:181392/117‑1 (MQ=255)
                                                                                                                                                 ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTaaa  <  1:881783/149‑1 (MQ=255)
                                                                                                                                                 ccGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTaaa  <  1:122609/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAA  >  NC_000913/2663319‑2663612

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