Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,123,637 T→A 5.8% E295V (GAG→GTG)  cpsB ← mannose‑1‑phosphate guanyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,123,6370TA 5.8% 255.5 / 2.8 88E295V (GAG→GTG) cpsBmannose‑1‑phosphate guanyltransferase
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base A (3/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCA  >  NC_000913/2123493‑2123786
                                                                                                                                                |                                                                                                                                                     
ttACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGACAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcg                                                                                                                                                  >  1:1126293/1‑150 (MQ=255)
  aCGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGgt                                                                                                                                                <  1:359085/150‑1 (MQ=255)
   cGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGGTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGgtg                                                                                                                                               >  1:405911/1‑150 (MQ=255)
    ggTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCAAGATTCTTTTTTTTTTTTTTTGTTTTGTTTTTTTTGTTTATTTACATCGCCGTGGCAAAGGTTGCCCTCGGCGgtgt                                                                                                                                              <  2:41486/150‑1 (MQ=255)
     gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGtgtg                                                                                                                                             >  2:909141/1‑150 (MQ=255)
     gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGtgtg                                                                                                                                             >  1:770473/1‑150 (MQ=255)
     gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGtgtg                                                                                                                                             <  2:291655/150‑1 (MQ=255)
     gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGAGGGCAACGTTTGCCACGGCGAtgtg                                                                                                                                             <  2:293920/150‑1 (MQ=255)
     gTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtg                                                                                                                                             <  1:1123303/150‑1 (MQ=38)
         gCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGgcg                                                                                                                                         <  2:1693829/150‑1 (MQ=255)
           aaTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                       >  1:695536/1‑150 (MQ=255)
           aaTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACACCGCCGAGGGCAACGTTTGCCACGGCGAtgtgattaat                                                                                                                                       <  1:1459286/150‑7 (MQ=37)
           aaTCAGCACCGCATCTTTGGGCTCCACCACAACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCt                                                                                                                                       <  2:299602/150‑1 (MQ=255)
            aTCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTg                                                                                                                                      <  2:104625/150‑1 (MQ=255)
             tCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGa                                                                                                                                     <  1:44099/150‑1 (MQ=255)
              cAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGTCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAt                                                                                                                                    >  2:1194823/1‑150 (MQ=255)
                 cACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGAtctc                                                                                                                                 <  2:60550/150‑1 (MQ=255)
                   ccGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                               >  1:1470890/1‑150 (MQ=255)
                   ccGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTccc                                                                                                                               <  2:750052/150‑1 (MQ=255)
                            tGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAg                                                                                                                      <  1:947170/150‑1 (MQ=255)
                               tCTGCACCACTACCAGATCTATCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGAATCATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAcc                                                                                                                   <  2:204037/150‑1 (MQ=255)
                               gcTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGAcc                                                                                                                   <  2:1331107/149‑1 (MQ=255)
                                cTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGCCCAGGCCAGATTCAGCATACCCATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGCAGAccc                                                                                                                  >  1:875090/1‑149 (MQ=255)
                                  gCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATATCCCATAATGAAGACAAgg                                                                                                                >  2:1409423/1‑150 (MQ=255)
                                    accacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAg                                                                                                              >  2:413659/1‑150 (MQ=255)
                                     ccacTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGc                                                                                                             >  1:360035/1‑150 (MQ=255)
                                        cTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcgg                                                                                                                                                 >  1:10782/1‑111 (MQ=255)
                                        cTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTcggcgg                                                                                                                                                 <  2:10782/111‑1 (MQ=255)
                                           ccAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGCCCAGGAGCCAACaa                                                                                                       >  2:1621071/1‑149 (MQ=255)
                                                tCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTc                                                                                                  >  1:1315451/1‑150 (MQ=255)
                                                 cTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGATTTCAGTTTTGTGATTAATCACATCGCCGTTGCAAACTTTGCCCTCGGCGGTGTGGGAGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTcc                                                                                                 >  2:247447/1‑150 (MQ=255)
                                                  tttCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCa                                                                                                <  1:453196/150‑1 (MQ=255)
                                                   ttCCCCCCGACGGTGGTGAACAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAg                                                                                               <  2:1063996/150‑1 (MQ=255)
                                                   ttCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAg                                                                                               <  2:545427/150‑1 (MQ=255)
                                                    tCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGc                                                                                              <  2:653301/150‑1 (MQ=255)
                                                     caccacGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATAAATCACATCGCCGTGTCAAACGTTGACCTCGGCGGTGTGGGCGCTGATCTACCATAATCAAGACCACGATCCAACATAGCTCCATcc                                                                                             >  2:349385/1‑150 (MQ=255)
                                                       ccccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCg                                                                                           <  1:1083925/150‑1 (MQ=255)
                                                       ccccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCg                                                                                           >  2:1083925/1‑150 (MQ=255)
                                                       ccccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCGGTTTACAGATTTGTGATTAATCACATAGCCGTGGCAAACGTTGCCCTCGGCGCTGTGGGCGCTGATCTCACAGAATGAAGACCAGCAGCCACCATCGCACCAGCCCg                                                                                           >  1:1548650/1‑150 (MQ=255)
                                                         ccGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGcc                                                                                         >  2:1731831/1‑149 (MQ=255)
                                                           gacgtgggtgAGCAGGCCAGATTCAGCATACAGATAGCTGTTTTCATTTTTGTGTTTAATCAAATCGCCGTGGCAAACGTGGCCCTCGGCGGTGTGGGAGCTGATCTCCCATAATGAAGACCAGGAGCCTACATCGCTCCAGCCCGcatc                                                                                       <  2:1548650/144‑1 (MQ=255)
                                                           gACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGCTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCATCGTCGTTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGcatc                                                                                       <  1:679356/150‑1 (MQ=255)
                                                                  gtgACCAGGCCAGATTCAGCTTACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGc                                                                                <  1:653733/150‑1 (MQ=255)
                                                                      ccAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGcacca                                                                            >  1:1677583/1‑150 (MQ=255)
                                                                       cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATGAGATCGCCGTGGCATACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGTAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGcacgac                                                                           <  2:751392/150‑1 (MQ=255)
                                                                       cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCATCATCGCTCCAGCCCGCATCCATCGGcaccac                                                                           >  1:1540006/1‑150 (MQ=255)
                                                                       cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGcaccac                                                                           <  2:36281/150‑1 (MQ=255)
                                                                       cAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGcaccac                                                                           <  1:1011313/150‑1 (MQ=255)
                                                                        aGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACcaca                                                                          <  2:1374094/150‑1 (MQ=255)
                                                                        aGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACcaca                                                                          <  1:271445/150‑1 (MQ=255)
                                                                          gCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCacaac                                                                        >  1:1300872/1‑150 (MQ=255)
                                                                            cAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCAGCGGCCCCacaacac                                                                      >  1:118923/1‑149 (MQ=255)
                                                                                  cAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTg                                                                <  2:1470890/150‑1 (MQ=255)
                                                                                  cAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCAGAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACCGCATCTg                                                                >  2:683855/1‑150 (MQ=255)
                                                                                    gCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGACACCACAACAGCATCTGcc                                                              <  2:696995/150‑1 (MQ=255)
                                                                                               gCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACCTCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                                   >  2:1559737/1‑150 (MQ=255)
                                                                                               gCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCAATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCAt                                                   >  1:1124900/1‑150 (MQ=255)
                                                                                                cTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGACAACATCGCTCCAGCCCGCATCCATCGGCACAACAACAGCATCGGCCGTACGTTCCATg                                                  >  1:683409/1‑150 (MQ=255)
                                                                                                 tGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACCACCGCAACTGCCGTACGTTCCATGa                                                 >  2:250712/1‑150 (MQ=255)
                                                                                                    tttCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACcg                                              >  2:1262424/1‑150 (MQ=255)
                                                                                                       cAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGt                                           <  1:366888/150‑1 (MQ=255)
                                                                                                        aGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCTTCGGCAGCACAACAGCATCTGCCGTACGTTCCATGACCGCGTa                                          >  1:97771/1‑150 (MQ=255)
                                                                                                         gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCCGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACACCAGCATCTGCCGTACGGGCCATGCCCGCGTaa                                         >  2:1313051/1‑150 (MQ=255)
                                                                                                         gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCCTCTGCCGTACGTAGCATGACCGCGTaa                                         >  2:1014018/1‑150 (MQ=255)
                                                                                                         gTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTaa                                         >  2:671189/1‑150 (MQ=255)
                                                                                                           tttGTGATTAATCCCATCGCCGTGGCAAACCTTGCCCACGGCGGTGTGGGCGCAGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCCTCCATCTGCACCACAACAGCATCTGCCGTACGTTCCATGACCGAGTAATc                                       >  1:146834/1‑150 (MQ=255)
                                                                                                           tttGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATc                                       <  1:413659/150‑1 (MQ=255)
                                                                                                             tgtgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAACATCTGCCGTACGTTCCATGACCGCGTAATCCa                                     <  2:1432733/150‑1 (MQ=255)
                                                                                                              gtgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCCCCCCAACAGCATCTGCCGTACGTTCCATGACCGCGTAATcccc                                    >  1:1119564/1‑148 (MQ=255)
                                                                                                              gtgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAc                                    >  1:3012/1‑150 (MQ=255)
                                                                                                               tgATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCAcc                                   <  1:1194823/150‑1 (MQ=255)
                                                                                                                    aaTCACATCTCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGAGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGACGTACGTTCCATGACCGCGTAATCCCCCGActc                              >  2:199576/1‑150 (MQ=255)
                                                                                                                     aTCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTTATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGGTCCAGGCCCGCGTAATCCCCCGActct                             >  1:1063517/1‑150 (MQ=255)
                                                                                                                         caTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTCCGTTCCATGACCGCGTAATCCACCGACTCTTCCg                         >  1:927432/1‑150 (MQ=255)
                                                                                                                          aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCCTCCATCGGCACCACAACAGCATCTGCCGTACGTTCCTTGACCGCGTAATCCCCCGACTCTTCCgg                        >  1:1017079/1‑150 (MQ=255)
                                                                                                                          aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAAt                                        >  2:1522601/1‑134 (MQ=255)
                                                                                                                          aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAAt                                        <  1:1522601/134‑1 (MQ=255)
                                                                                                                          aTCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCAGCGGCACCACAACCGCATCTGCCGTACGTTCCATGCCCGCGTAATCCACCGACTCTTCCgg                        >  1:1226875/1‑150 (MQ=255)
                                                                                                                           ccGCCGAGGGCAACGTTTGCCACGGCGATGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCggg                       >  1:138643/2‑150 (MQ=255)
                                                                                                                              ccGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCCCCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAg                    >  2:1621278/1‑150 (MQ=255)
                                                                                                                                   gCAAACGTTGCCCACGGCGGTGTGGGCGCTGAACACCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCCCAACCGCCACTGCCGTACGTACCATGCCCGCGTATACCACCGCCTCAGCCGGGCCGGCgag               >  1:1126886/1‑150 (MQ=255)
                                                                                                                                        cGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACg          >  2:645442/1‑150 (MQ=255)
                                                                                                                                            gCCCTCTGCGGTGTGGGCGCTGATCTCCCATCATGACGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCCCGTAATCCACCGACTCTTCCGGCCAGGCGAGAAACGcttc      <  1:349385/150‑1 (MQ=255)
                                                                                                                                            gCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGcttc      <  1:183469/150‑1 (MQ=255)
                                                                                                                                             cccTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGcttct     <  2:770473/150‑1 (MQ=255)
                                                                                                                                             cccTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGcttct     <  2:1126293/150‑1 (MQ=255)
                                                                                                                                              ccTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACAGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGcttctt    <  1:247447/150‑1 (MQ=255)
                                                                                                                                               cTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCttcttc   <  2:695536/150‑1 (MQ=255)
                                                                                                                                                tCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCATAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCa  <  1:1715291/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTACGGTCGGCAATCAGCACCGCATCTTTGGTCTGCACCACTACCAGATCTTTCACCCCGACGGTGGTGACCAGGCCAGATTCAGCATACACATAGCTGTTTTCAGTTTTGTGATTAATCACATCGCCGTGGCAAACGTTGCCCTCGGCGGTGTGGGCGCTGATCTCCCATAATGAAGACCAGGAGCCAACATCGCTCCAGCCCGCATCCATCGGCACCACAACAGCATCTGCCGTACGTTCCATGACCGCGTAATCCACCGACTCTTCCGGGCAGGCGAGAAACGCTTCTTCA  >  NC_000913/2123493‑2123786

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