Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,681,500 G→A P416S (CCT→TCT)  yphG ← DUF5107 domain‑containing protein YphG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,681,5000GA98.2% 181.3 / ‑4.4 54P416S (CCT→TCT) yphGDUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/31);  total (22/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 = 9.04e-01

CCAGCGCCACGCCGCGCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGGATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATAAGGCAGGAAATACTGCTCAAAAC  >  NC_000913/2681275‑2681739
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
ccAGCGCCACGCCGCGCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTg                                                                                                                                                                                                                          <  1:120474/250‑1 (MQ=255)
  aGCGCCACGCCGCGCAGGTAGTACTCGAACGGTGAACGGCAGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAAGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGGTCTCAATGGTCAGCCTTTCAGAATTGATACAGTGCAaccccccatggat                                                                                                                                                                                                                       >  2:181721/1‑240 (MQ=255)
     gcCACGCCGCTCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGc                                                                                                                                                                                                                    <  1:198253/251‑1 (MQ=255)
     gcCACGCCGCGCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGc                                                                                                                                                                                                                    >  2:193935/1‑251 (MQ=255)
       cACGCCGCGCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGCCGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCACTGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCgg                                                                                                                                                                                                                  >  1:85229/1‑251 (MQ=255)
        acgccgatcAGGTAGTAATCGAACGGTGAACGGCCGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGAGGCGATTTGGCAACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGCAAAATATTGCCGTCGGCATCGGAGAGCTCAACGGTCAGCCTTTCAGAATTGATACCGTGCAACCCGCCCTGGATGGCggg                                                                                                                                                                                                                 >  2:221159/9‑250 (MQ=255)
                    tagtaATCGAAAGGTGAACGGCTGGCGGGCTGATCTTGCTCAAGATGCTGACCGCTAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCAAGTGGCGCTTTGGCGACGACCGGTACCGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGCGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCACAACGCCCTGGATGGCGGTCGCAGGCATc                                                                                                                                                                                                       >  2:44557/1‑249 (MQ=255)
                           cGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACg                                                                                                                                                                                              <  1:48118/251‑1 (MQ=255)
                            gAACGGTGAACGGCAGGCGTGCTGATCTTGCTCCAGATACTGACCGATAAACCAGGCTTCATATGTACTGGTAATGTCTTGGGCTGCCCGTAGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTCCAATATTGCCGTCGGCCTCGGAGAGCTCAATGATCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCAtacgtgccccgg                                                                                                                                                                                             >  2:315134/1‑240 (MQ=255)
                                ggTGAACGGCTGGGGTGATGATATTGCTCCAGATGCTGACCGATAAGCCAGGCTGCAGCTGTACTGGTAATGTCGTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAGTTTATCCCGTGCAACACGGTCTGGAGGGCGGTCGCAGGCATCAGTGGCACGGcatc                                                                                                                                                                                         <  2:159562/251‑1 (MQ=255)
                                  tGAACGGCTAGCGTGATGAAATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATATGTACTGGTAATGTCTTGTGCTGCCAGAGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCACTGGATGGCGGACGCAGGCATCAAGGCCACGGcatcat                                                                                                                                                                                      >  2:306993/1‑251 (MQ=255)
                                    aaCGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCAtcatca                                                                                                                                                                                    <  1:119673/251‑1 (MQ=255)
                                       ggCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGt                                                                                                                                                                                 >  2:121901/1‑251 (MQ=255)
                                        gCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATTGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTa                                                                                                                                                                                <  1:63666/251‑1 (MQ=255)
                                        gCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAgg                                                                                                                                                                                                           <  2:19139/225‑1 (MQ=255)
                                        gCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAgg                                                                                                                                                                                                           >  1:19139/1‑225 (MQ=255)
                                                            ccAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGAt                                                                                                                                                             <  1:37434/250‑1 (MQ=255)
                                                                  gCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCgg                                                                                                                                                      <  2:285934/251‑1 (MQ=255)
                                                                     gACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCggctc                                                                                                                                                   >  2:171779/1‑248 (MQ=255)
                                                                        cGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGcc                                                                                                                                                >  1:334882/1‑251 (MQ=255)
                                                                         gATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCa                                                                                                                                               <  2:282997/251‑1 (MQ=255)
                                                                                      tCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACACTATTGCCGTCGGCATCGGAGCGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGTATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCACAAGGCGGTATCCGtt                                                                                                                                  >  2:285845/1‑251 (MQ=255)
                                                                                            gTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGGTCAAc                                                                                                                              <  1:106181/249‑1 (MQ=255)
                                                                                             gaCTGGTAATGTCTTGTGCTGCCATTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGga                                                                                                                           <  2:166860/250‑1 (MQ=255)
                                                                                                tGGTAATATCTTGAGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTAAAATATTGCCGTCGGCATCGGAGCGCTCAATGGTAAGCCTTTCAGAATTGATACCGTGCAACACGCCTTGGATGGCGTGCGCAGGCATCAGTGCCACGGCTCAATCAAGTAACGCGTTGCATTTACCGATTTCGCGGAgcgcacggcgggatcagtacaaagcaga                                                                                                                         >  2:12073/1‑226 (MQ=255)
                                                                                                           tgtgCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACa                                                                                                              <  2:203109/250‑1 (MQ=255)
                                                                                                              gctgcCAGTGGCGCTTTGGCAACGTCCGGTAACGGCACCGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGACGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGccc                                                                                                          >  1:226400/1‑251 (MQ=255)
                                                                                                                  ccAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGCTAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCAc                                                                                                      >  2:149933/1‑251 (MQ=255)
                                                                                                                  ccAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCACTGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCAc                                                                                                      >  2:49554/1‑251 (MQ=255)
                                                                                                                    aGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCCTACAGCCCCCACTc                                                                                                    >  2:155908/1‑251 (MQ=255)
                                                                                                                     gTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTc                                                                                                    >  1:55446/1‑250 (MQ=255)
                                                                                                                      tGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgagag                                                                                                                        >  2:246597/1‑229 (MQ=255)
                                                                                                                      tGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGgagag                                                                                                                        <  1:246597/229‑1 (MQ=255)
                                                                                                                                    cGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTtacta                                                                                    <  1:143570/251‑1 (MQ=255)
                                                                                                                                         ggTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTGTCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTg                                                                               <  1:111058/251‑1 (MQ=255)
                                                                                                                                            aaCGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGa                                                                             <  1:275406/250‑1 (MQ=255)
                                                                                                                                                         gCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGc                                                                <  1:195932/250‑1 (MQ=255)
                                                                                                                                                          cGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCAt                                                              >  1:55039/1‑251 (MQ=255)
                                                                                                                                                            gCTGATGTTCCTGATAACTCAGTACAATATTTCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATTCCGTTCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCAGTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATGCAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATcg                                                            <  1:200530/251‑1 (MQ=255)
                                                                                                                                                                  gTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAg                                                      <  1:155908/251‑1 (MQ=255)
                                                                                                                                                                           aaCTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAAGTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTGCGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGc                                             <  2:136187/251‑1 (MQ=255)
                                                                                                                                                                            aCTCAGTACAATATTGCCGTCGGCAACGGAGAGCTAAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCGCGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATTGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCa                                            <  2:261155/251‑1 (MQ=255)
                                                                                                                                                                              ccAGTACAATAGTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCTACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCAcc                                          <  2:109537/250‑1 (MQ=255)
                                                                                                                                                                                aGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCAt                                        <  1:193935/251‑1 (MQ=255)
                                                                                                                                                                                 gTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATg                                       >  1:37292/1‑251 (MQ=255)
                                                                                                                                                                                  gaCAATATTGCCTTCGGCATCGGAGAGCTCAATGTTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCATTGCCACGGCATCATCAAGTAACGCGTTGCATTTACGTATGTCGCAGATCGCCATGCTGTGTCCGTACCTCTGAGAGTTGGCATACAGCCCCCACTCAATCCCCCTCTTACTACTCTGGAGGTTGATCACCGCATCTCTGGAGGCATTTTGcgctatgc                                      <  1:56936/250‑8 (MQ=255)
                                                                                                                                                                                   aCAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACTCTGGAGGTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGc                                      <  1:166753/250‑1 (MQ=255)
                                                                                                                                                                                           gCCGTCGGCATCGGAGAGCGCAATGGTCAGCCTGTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGaa                              <  2:190008/250‑1 (MQ=255)
                                                                                                                                                                                                cGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGGATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATa                         <  1:241149/250‑1 (MQ=255)
                                                                                                                                                                                                cGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTGCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTATTATCACCGCATCGCGGGAGGCATTTTGCACCAGGCCCAAAGAATGATa                         <  1:306993/250‑1 (MQ=255)
                                                                                                                                                                                                   cATCGGAGAGCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGGGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTGCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCTCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATAAGGc                     <  1:343241/251‑1 (MQ=255)
                                                                                                                                                                                                           agCTCAATGGTCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCCGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCACTCCCCCGCTTACTCCGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTTCACCATCCCCATAGAATGATAAGAAAGGAAATa             <  1:27491/251‑1 (MQ=255)
                                                                                                                                                                                                                     tCAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATAAGGCAGGCAATACTGCTCaaa    >  2:221591/1‑250 (MQ=255)
                                                                                                                                                                                                                      cAGCCTTTCAGAATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATAAGGCAGGAAATACTGCTCAAAAc  <  1:149933/251‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
CCAGCGCCACGCCGCGCAGGTAGTAATCGAACGGTGAACGGCTGGCGTGATGATATTGCTCCAGATGCTGACCGATAAACCAGGCTTCATCTGTACTGGTAATGTCTTGTGCTGCCAGTGGCGCTTTGGCGACGTCCGGTAACGGCAACGCTTGCGGCTGATGTTCCTGATAACTCAGTACAATATTGCCGTCGGCATCGGAGAGCTCAATGGTCAGCCTTTCAGGATTGATACCGTGCAACACGCCCTGGATGGCGGTCGCAGGCATCAGTGCCACGGCATCATCAAGTAACGCGTTGCATTTACCGATTTCGCGGATCGCCAGGCGGTATCCGTTCAACGGAGAGATGGCATACAGCCCCCACTCAATCCCCCGCTTACTACGCTGGAGTTTTATCACCGCATCGCGGGAGGCATTTTGCACCATGCCCAAAGAATGATAAGGCAGGAAATACTGCTCAAAAC  >  NC_000913/2681275‑2681739

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