Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,908,486 G→C 7.7% A394G (GCC→GGC)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,908,4860GC 7.7% 642.7 / 32.3 183A394G (GCC→GGC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base G (74/94);  new base C (8/6);  total (83/100)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 189 total aligned reads displayed.

TGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTC  >  NC_000913/2908338‑2908634
                                                                                                                                                    |                                                                                                                                                    
tgCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTgg                                                                                                                                                      <  4:1536736/149‑1 (MQ=255)
 gcgcCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGc                                                                                                                                                     <  3:1727015/149‑1 (MQ=255)
      gagaCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACAt                                                                                                                                                >  1:1779390/1‑149 (MQ=255)
         aCGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGCCAACATgga                                                                                                                                             <  3:1605513/149‑2 (MQ=255)
           gCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           >  3:423643/1‑149 (MQ=255)
             aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGc                                                                                                                                         >  1:1065896/1‑149 (MQ=255)
             aTGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGc                                                                                                                                         >  1:70438/1‑149 (MQ=255)
              tGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCg                                                                                                                                        <  3:152742/149‑1 (MQ=255)
                 tACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTaa                                                                                                                                     <  2:1128394/149‑1 (MQ=255)
                 tACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTaa                                                                                                                                     <  1:716480/149‑1 (MQ=255)
                     gccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCaa                                                                                                                                 <  1:729068/149‑1 (MQ=255)
                      ccgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAt                                                                                                                                <  3:871404/149‑1 (MQ=255)
                       cgAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAAtt                                                                                                                               <  3:220552/149‑1 (MQ=255)
                          gaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                            <  4:1507323/149‑1 (MQ=255)
                          gaTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                            <  4:10856/149‑1 (MQ=255)
                           aTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           <  1:1145192/133‑1 (MQ=255)
                           aTCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           >  2:1145192/1‑133 (MQ=255)
                             cGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACa                                                                                                                                                 <  4:178667/125‑1 (MQ=255)
                               ctctTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAc                                                                                                                       <  2:1395351/149‑1 (MQ=255)
                                tcttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCAcc                                                                                                                      >  4:624885/1‑149 (MQ=255)
                                  ttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCaa                                                                                                                                                   <  4:128560/118‑1 (MQ=255)
                                  ttctCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCaa                                                                                                                                                   >  3:128560/1‑118 (MQ=255)
                                         tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGTCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                             >  4:1077970/1‑149 (MQ=255)
                                         tACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                             >  3:1106526/1‑149 (MQ=255)
                                           cGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAg                                                                                                           <  2:1065896/149‑1 (MQ=255)
                                           cGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAg                                                                                                           >  1:793783/1‑149 (MQ=255)
                                               cTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCaa                                                                                                       >  3:1557800/1‑149 (MQ=255)
                                                ttCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAAt                                                                                                                                    <  3:1613204/119‑1 (MQ=255)
                                                ttCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAAt                                                                                                                                    >  4:1613204/1‑119 (MQ=255)
                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAAt                                                                                                     >  1:1040832/1‑149 (MQ=255)
                                                 tCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAAt                                                                                                     >  1:920983/1‑149 (MQ=255)
                                                      gTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCa                                                                                                <  4:452956/149‑1 (MQ=255)
                                                         gCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                             <  2:793783/149‑1 (MQ=255)
                                                          ccGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGAt                                                                                             <  1:85362/148‑1 (MQ=255)
                                                            gTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCTATGTTGCCAACATGgagaa                                                                                                                                          <  3:1568600/101‑2 (MQ=255)
                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           >  1:396052/1‑99 (MQ=255)
                                                             ttttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           <  2:396052/99‑1 (MQ=255)
                                                              tttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGa                                                                                        <  3:954576/149‑1 (MQ=255)
                                                               ttCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGaa                                                                                       <  3:354749/149‑1 (MQ=255)
                                                                   tCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATAtt                                                                                   <  3:604730/149‑1 (MQ=255)
                                                                    cgcgCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTg                                                                                  <  1:453850/149‑1 (MQ=255)
                                                                         cACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGtt                                                                                >  1:700113/1‑146 (MQ=255)
                                                                         cACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGtt                                                                                <  2:700113/146‑1 (MQ=255)
                                                                          aCTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATcccc                                                                                                             >  3:1672821/1‑116 (MQ=255)
                                                                          aCTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATAccc                                                                                                             <  4:1672821/116‑1 (MQ=255)
                                                                          aCTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCa                                                                            >  2:1664345/1‑149 (MQ=255)
                                                                           cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCTAGGCAAATGCCCAGATAAg                                                                                          <  3:1119107/134‑1 (MQ=255)
                                                                           cTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCAc                                                                           >  4:1709116/1‑149 (MQ=255)
                                                                               gTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGa                                                                                                                                           <  2:859525/81‑1 (MQ=255)
                                                                               gTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACgcgc                                                                       >  1:362774/1‑149 (MQ=255)
                                                                                    cAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAg                                                                                                           >  2:360700/1‑108 (MQ=255)
                                                                                    cAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAc                                                                  <  4:1557800/149‑1 (MQ=255)
                                                                                     aGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                           <  4:27518/123‑1 (MQ=255)
                                                                                      gcgcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                            <  3:291879/89‑1 (MQ=255)
                                                                                      gcgcCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgc                                                                >  1:643217/1‑149 (MQ=255)
                                                                                         ccACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAAcgc                                                                <  3:1055560/146‑1 (MQ=255)
                                                                                          cacaACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCgg                                                            >  1:1867320/1‑149 (MQ=255)
                                                                                          cacaACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCgg                                                            >  4:1200729/1‑149 (MQ=255)
                                                                                            caACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGc                                                                                                   >  3:895163/1‑108 (MQ=255)
                                                                                            caACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGc                                                                                                   <  4:895163/108‑1 (MQ=255)
                                                                                            caACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTc                                                          <  4:423643/149‑1 (MQ=255)
                                                                                             aaCCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATAtt                                                                                   <  1:1303403/123‑1 (MQ=255)
                                                                                                cGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTaa                                                      <  3:533322/149‑1 (MQ=255)
                                                                                                    tACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCAt                                                  <  2:362774/149‑1 (MQ=255)
                                                                                                      cTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCa                                                                                                                        >  3:719458/1‑77 (MQ=255)
                                                                                                          cAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTc                                            <  2:643217/149‑1 (MQ=255)
                                                                                                             tCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTAc                                         <  3:1689204/149‑1 (MQ=255)
                                                                                                              cTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACg                                        <  2:1867320/149‑1 (MQ=255)
                                                                                                               tGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAAc                                                                                                                            >  1:1655921/1‑64 (MQ=255)
                                                                                                               tGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGc                                       >  4:1118061/1‑149 (MQ=255)
                                                                                                                 gCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCa                                     >  2:740642/1‑149 (MQ=255)
                                                                                                                   acaAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                          <  1:21012/94‑1 (MQ=255)
                                                                                                                    caAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGa                                  <  2:222508/149‑1 (MQ=255)
                                                                                                                     aaaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                           <  3:1378036/91‑1 (MQ=255)
                                                                                                                     aaaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATaa                                                                                           >  4:1378036/1‑91 (MQ=255)
                                                                                                                      aaTTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATa                                <  3:76712/149‑1 (MQ=255)
                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTaa                                                      >  2:1100494/1‑121 (MQ=255)
                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTaa                                                      <  1:1100494/121‑1 (MQ=255)
                                                                                                                            gTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGaa                          >  2:566355/1‑149 (MQ=255)
                                                                                                                                 gTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGc                              <  1:577380/140‑1 (MQ=255)
                                                                                                                                 gTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGc                              >  2:577380/1‑140 (MQ=255)
                                                                                                                                  ttGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGcc                    <  4:943496/149‑1 (MQ=255)
                                                                                                                                       gTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTa               >  4:1014223/1‑149 (MQ=255)
                                                                                                                                         tcgccATTTTGCCAACATGGAGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg             >  2:1033861/4‑149 (MQ=255)
                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTa                                          >  4:350840/4‑120 (MQ=255)
                                                                                                                                         tcgccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACg             >  1:1520149/4‑149 (MQ=255)
                                                                                                                                           tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCtt                                             <  4:572012/115‑1 (MQ=255)
                                                                                                                                           tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGa            <  2:384743/148‑1 (MQ=255)
                                                                                                                                           tcCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGa            >  1:384743/1‑148 (MQ=255)
                                                                                                                                           gccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGc                     >  2:804203/2‑139 (MQ=255)
                                                                                                                                            ccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATc                                                    <  3:806582/107‑1 (MQ=255)
                                                                                                                                            ccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATc                                                    >  4:806582/1‑107 (MQ=255)
                                                                                                                                            ccATGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAgg          >  1:339808/1‑149 (MQ=255)
                                                                                                                                              aTGTTGCCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGAt        <  2:1520149/149‑1 (MQ=255)
                                                                                                                                               tGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAg                                                                                          <  2:1246185/66‑1 (MQ=255)
                                                                                                                                               tGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACCCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGAtt       >  3:315361/1‑149 (MQ=255)
                                                                                                                                                 ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCa                   <  1:1472646/135‑1 (MQ=255)
                                                                                                                                                 ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc     <  4:733636/149‑1 (MQ=255)
                                                                                                                                                 ttGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGc     >  3:1769587/1‑149 (MQ=255)
                                                                                                                                                    gCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTc  <  4:315361/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGCCGAGACGCATGGTACCGCCGAGATCGCTCTTCTCGCTACGAACTTCAACGTTGCCGTTTTCATCGCGCCACTCGGTAATCAGCGCCACAACCGGGTACTTACAGTCTGGCACAAATTCCGTAGAGTTGGCGTTCTCCATGTTGGCAACATGGCGAGCGTAATCAATTAACGCCACCTGCATACCCAGGCAAATGCCCAGATAAGGAATATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTC  >  NC_000913/2908338‑2908634

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