Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,307,298 G→A 21.9% P168S (CCC→TCC)  alsK ← D‑allose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,307,2980GA21.9% 116.0 / 28.0 60P168S (CCC→TCC) alsKD‑allose kinase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base A (6/7);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

AGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCA  >  NC_000913/4307149‑4307447
                                                                                                                                                     |                                                                                                                                                     
aGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAgggg                                                                                                                                                       <  2:622228/150‑1 (MQ=255)
 ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTTGGTCATATCTCCCAggggg                                                                                                                                                      <  1:1313929/150‑1 (MQ=255)
 ggggCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAggggg                                                                                                                                                      <  1:293139/150‑1 (MQ=255)
      gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACTCGCGGAGTGTTTGGTCATATCTCCCAGGGGGATATg                                                                                                                                                 <  2:606423/150‑1 (MQ=255)
      gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATg                                                                                                                                                 >  2:948268/1‑150 (MQ=255)
        tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAc                                                                                                                                               <  1:403115/150‑1 (MQ=255)
        tttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAc                                                                                                                                               <  1:383632/150‑1 (MQ=255)
         ttCCGCATGGACGAAAAGATCGCTCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGAGAGATATGAcc                                                                                                                                              <  1:155084/150‑1 (MQ=255)
         ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGAcc                                                                                                                                              <  1:543295/150‑1 (MQ=255)
          tCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGAccc                                                                                                                                             <  2:886292/150‑1 (MQ=255)
           ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGGCCCa                                                                                                                                            <  1:293365/150‑1 (MQ=255)
           ccGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCACTTGTAGATATGAACCa                                                                                                                                            <  2:1304714/150‑1 (MQ=255)
           ccGCATGGACGAAAAGATCGCGCAATGGGTAAATTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCa                                                                                                                                            <  1:971850/150‑1 (MQ=255)
              cATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGtt                                                                                                                                         >  1:1028775/1‑150 (MQ=255)
                tGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATGTGACCCAGTTCg                                                                                                                                       >  2:972698/1‑150 (MQ=255)
                  gACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGcc                                                                                                                                     <  1:965580/150‑1 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGTCCCAGTTCGCCTg                                                                                                                                   >  1:641458/1‑150 (MQ=255)
                    cGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTg                                                                                                                                   <  1:626796/150‑1 (MQ=255)
                              gcgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCgtg                                                                                                                         <  1:514994/150‑1 (MQ=255)
                                gcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGAAGTGTTGGGTCATATCTCCCAGGGTGATATGACCCAGTTCGCCTGCCACACCGtgtg                                                                                                                       <  2:1503056/150‑1 (MQ=255)
                                  aaTGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCa                                                                                                                     <  1:1114305/150‑1 (MQ=255)
                                    tGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAcc                                                                                                                   <  2:1481457/150‑1 (MQ=255)
                                    tGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAcc                                                                                                                   <  2:1321595/150‑1 (MQ=255)
                                           tttCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCAc                                                                                                            >  1:796409/1‑150 (MQ=255)
                                                   cTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCACGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAcc                                                                                                    <  1:818053/150‑1 (MQ=255)
                                                       tgttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTc                                                                                                >  1:1525188/1‑150 (MQ=255)
                                                        gttCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATGTGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                                               <  1:972698/150‑1 (MQ=255)
                                                          tCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATc                                                                                             >  2:1403670/1‑150 (MQ=255)
                                                          tCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATc                                                                                             <  1:616889/150‑1 (MQ=255)
                                                                 aGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTg                                                                                      >  2:1102378/1‑150 (MQ=255)
                                                                        tAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                               <  1:546649/150‑1 (MQ=255)
                                                                        tAGCGCCATTCCAGAGCAATTGGATTCAAGGCACCCAGGATTGCCACACGCGCAGTGTAGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                                               >  2:607722/1‑150 (MQ=255)
                                                                          gcgcCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCa                                                                             >  1:389229/1‑150 (MQ=255)
                                                                             ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcc                                                                          <  1:485668/150‑1 (MQ=255)
                                                                                 tCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCCTCCCCGt                                                                      >  1:1437880/1‑150 (MQ=255)
                                                                                   cAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCTGGGGGATATGACCCAGTTCGCCTGCCACATTGTGTGATACTGTTAAGGGCTCACCGTTGATCTACACTGCGAACCCCATCCCCGTAc                                                                    <  2:1015996/150‑1 (MQ=255)
                                                                                      agCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGTCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCga                                                                 <  2:641458/150‑1 (MQ=255)
                                                                                       gCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCgag                                                                <  1:812330/150‑1 (MQ=255)
                                                                                         aaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCTCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGAt                                                              <  1:509519/150‑1 (MQ=255)
                                                                                         aaTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGAt                                                              <  2:1525188/150‑1 (MQ=255)
                                                                                                  ccAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCa                                                     >  1:429756/1‑150 (MQ=255)
                                                                                                     ggCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACCCCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCATCCCCGTCCCGAGATAGGCACCCAgca                                                  >  2:1228058/1‑148 (MQ=255)
                                                                                                          ccAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCAACCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgt                                             <  2:1061082/150‑1 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                     >  1:1151127/1‑150 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGt                                     <  1:860588/150‑1 (MQ=255)
                                                                                                                  gCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCCTTCCTCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGACCCAGGTGTTGCGt                                     >  2:972224/1‑150 (MQ=255)
                                                                                                                    cacaCGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACAGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTaa                                   >  1:1406264/1‑150 (MQ=255)
                                                                                                                      cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                 >  2:893883/1‑150 (MQ=255)
                                                                                                                      cacGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATACCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg                                 >  2:250388/1‑150 (MQ=255)
                                                                                                                              gtgtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc                         >  1:47762/1‑150 (MQ=255)
                                                                                                                                gtTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTcta                       <  1:262787/150‑1 (MQ=255)
                                                                                                                                 ttGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                      <  1:948268/150‑1 (MQ=255)
                                                                                                                                 ctGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctac                      >  1:1315339/2‑150 (MQ=255)
                                                                                                                                  tGGGTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCTCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctact                     >  1:1130700/1‑150 (MQ=255)
                                                                                                                                   gggTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                    <  2:429756/150‑1 (MQ=255)
                                                                                                                                   gggTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTctacta                    >  1:2303/1‑150 (MQ=255)
                                                                                                                                     gTCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACg                  >  2:1087378/1‑150 (MQ=255)
                                                                                                                                      tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCTCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGt                 <  2:1130700/150‑1 (MQ=255)
                                                                                                                                      tCATATCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGt                 >  2:943882/1‑150 (MQ=255)
                                                                                                                                         tatCCCCCTGGGAGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCCCACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTccc              >  2:1353132/1‑150 (MQ=255)
                                                                                                                                             tcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGCACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGga          >  1:520943/1‑150 (MQ=255)
                                                                                                                                             tcCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGga          >  1:996163/1‑150 (MQ=255)
                                                                                                                                                cAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGagag       >  1:1465958/1‑150 (MQ=255)
                                                                                                                                                     ggATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCa  <  1:1403670/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGGGGCGTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCCCAGGGGGATATGACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCA  >  NC_000913/4307149‑4307447

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