Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,414,192 G→A intergenic (+1/‑83) yjfC → / → aidB putative acid‑‑amine ligase YjfC/putative acyl‑CoA dehydrogenase AidB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,414,1920GA90.2% 131.8 / 4.2 51intergenic (+1/‑83)yjfC/aidBputative acid‑‑amine ligase YjfC/putative acyl‑CoA dehydrogenase AidB
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (20/26);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

AAGGCGGCAACGTCACCATTTTTGACGGCAAGAATAACGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAAT  >  NC_000913/4413950‑4414437
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
aaGGCGGCAACGTCACCATTTTTGACGGCAAGAATAACGTTGTTGACCACGCTGATGGGGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGtttt                                                                                                                                                                                                                                               <  2:314634/251‑1 (MQ=255)
                  tttttGACGGCAAGAATAACGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTACGGAATGGGGATCCGTGACGATAACACAATGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAACTGTTTTAGCAATCTCTTTCTGTc                                                                                                                                                                                                                              >  1:454094/1‑250 (MQ=255)
                        aCGGCAAGAATAACGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCACCCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTAATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCCGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATc                                                                                                                                                                                                                       >  2:219878/1‑251 (MQ=255)
                        aCGGCAAGAATAACGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATc                                                                                                                                                                                                                       <  1:390787/251‑1 (MQ=255)
                                   aaCGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTg                                                                                                                                                                                                             >  1:138823/1‑250 (MQ=255)
                                                   cTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGAc                                                                                                                                                                                            >  1:354206/1‑251 (MQ=255)
                                                            aTTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTg                                                                                                                                                                                    >  1:143959/1‑250 (MQ=255)
                                                                      tGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATggaggg                                                                                                                                                                            >  1:191041/1‑248 (MQ=255)
                                                                        aaCCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGacc                                                                                                                                                                       >  2:434978/1‑250 (MQ=255)
                                                                         aCCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTGTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGacac                                                                                                                                                                      <  2:5145/251‑1 (MQ=255)
                                                                             atgatCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCAGCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTg                                                                                                                                                                  <  1:346555/251‑1 (MQ=255)
                                                                               gatCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGCCACAGTGc                                                                                                                                                                 >  1:367497/1‑250 (MQ=255)
                                                                                   tACCAGGCGTTTCAACCTCTGCCGCAGGTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATAGTGACTGCCATTGATGGAGGGAGACAAAGTGCACTgg                                                                                                                                                            >  1:496084/1‑251 (MQ=255)
                                                                                    aCCAGGCTGTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACt                                                                                                                                                              <  2:151627/248‑1 (MQ=255)
                                                                                    aCCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCAc                                                                                                                                                               <  1:248861/247‑1 (MQ=255)
                                                                                     ccAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACTCTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCACTCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGc                                                                                                                                                           >  1:123061/1‑250 (MQ=255)
                                                                                               cAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGGAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACAcc                                                                                                                                                <  1:305753/251‑1 (MQ=255)
                                                                                                       gCCGCGGTTTTGCGATAGCTACACACTCATCGGTAGCTGGAGTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACCCCGTTTTTaa                                                                                                                                        <  1:288995/251‑1 (MQ=255)
                                                                                                            ggTTTGGCGATATCTACACACTCATCTGTAGCTGGATTGGCGATGATGAAGAGTGCGGAATGGGGAGCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCCCTGGCAAACTCACACCGTTTTTAATCAAc                                                                                                                                   <  2:454094/251‑1 (MQ=255)
                                                                                                               ttGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCt                                                                                                                                 <  1:434978/250‑1 (MQ=255)
                                                                                                                tGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCGGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCta                                                                                                                                <  1:26255/250‑1 (MQ=255)
                                                                                                                   cGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTAATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCCGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAcc                                                                                                                            <  1:219878/251‑1 (MQ=255)
                                                                                                                    gATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATAcc                                                                                                                            >  1:267689/1‑250 (MQ=255)
                                                                                                                     aTAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCAt                                                                                                                          >  1:271031/1‑251 (MQ=255)
                                                                                                                       aGCTACACACTCATCGGTAGCAGTATGGGCGATGATTAAGCTTGCGGAATGTGGATCCGTTAAGATAACACACTGATCCACAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATTTCAGTGACCATACTAACGTTAACTGCTATTGATGGAGGGAGACACAGTGCACTGGAAAACTCACACCGTTTTTAATCAACCTATTCCAtt                                                                                                                         <  2:306422/250‑1 (MQ=255)
                                                                                                                         cTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTaa                                                                                                                       >  1:106927/1‑250 (MQ=255)
                                                                                                                                   aTCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAAc                                                                                                            >  1:109110/1‑251 (MQ=255)
                                                                                                                                        tAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTtat                                                                                                        >  2:136647/1‑250 (MQ=255)
                                                                                                                                             ggatttccGATGATTAGGCGTGCGGCATGGGGATCCGTGAAGATAACACACTGAGCACTAAAGACACCTCACGTTTCGTTCCGCATTACACTGCTGGATAAAAATGATGGAGCAATCTCTTTCTTGCATGAATCCGTGGCAGTGACCATACGAATGGTGACTGCTATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTGATCCCTGt                                                                                                  <  2:315022/244‑1 (MQ=255)
                                                                                                                                              gATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCACCTTATACCTGTc                                                                                                 >  2:214683/1‑251 (MQ=255)
                                                                                                                                                 tGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTACTGGATAAGAATGTTTTAGCAATCTCTTTCTg                                                                                                                                                                                                                                >  2:107068/1‑121 (MQ=255)
                                                                                                                                                 tGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTACTGGATAAGAATGTTTTAGCAATCTCTTTCTg                                                                                                                                                                                                                                <  1:107068/121‑1 (MQ=255)
                                                                                                                                                     gatgatGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGgc                                                                                          <  2:398377/251‑1 (MQ=255)
                                                                                                                                                            aaGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAc                                                                                                                                                                                                           >  2:165696/1‑131 (MQ=255)
                                                                                                                                                            aaGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGAc                                                                                                                                                                                                           <  1:165696/131‑1 (MQ=255)
                                                                                                                                                                    ggAATGGGGATCCGTGAAGATACCACACTGATCACCAAAGACACCTAACGTTTCGTTCAGCATTACAATGCTGGATCAAAATGATTTAGCAATCTCTTTCTGTCATGAATCCATGGCCGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGAAACAGTGCACTGGCACACTCACACAGTTTTTAATAAACATATCCCTTTAAATAACAGCAACTTATACATATATGATGGCGCGATCTGCGCAGCg                                                                           >  2:302998/1‑251 (MQ=255)
                                                                                                                                                                    ggAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                                                                           >  2:386149/1‑251 (MQ=255)
                                                                                                                                                                    ggAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCAATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATCGCGCGCTCTGCGAAGCg                                                                           <  2:143959/251‑1 (MQ=255)
                                                                                                                                                                    ggAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTCTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCg                                                                           >  1:76247/1‑251 (MQ=255)
                                                                                                                                                                                    aaGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGTCAGTGACCATACTAATGGTGACTGACATTGATGGAGGGAGACACAGTGCACTGTCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCGCTGCGAAGCGGTAACGCGTGAAGGTg                                                           <  1:415712/251‑1 (MQ=255)
                                                                                                                                                                                          acacacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTgctggct                                                     <  1:103775/251‑1 (MQ=255)
                                                                                                                                                                                             cacTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTggg                                                  <  2:124799/251‑1 (MQ=255)
                                                                                                                                                                                                                  cACGTTTCGTTCCGCATTACATTGCGGGATAAAATTGTTTTAGCAATCTCGTTTTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGATACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATTCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTCCTGGCTGGGATAGCGATTTTCTAGCCAGt                              <  2:344343/250‑1 (MQ=255)
                                                                                                                                                                                                                   ccGTTTCGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTAt                            <  1:220003/250‑1 (MQ=255)
                                                                                                                                                                                                                         cGTTCCGCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTca                      <  1:386149/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCATTACATTGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGACGGTGCTGGCTGGGCTAGCGATTTTCTAACCGGTATTGGTCAGCAGtt                >  2:100751/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cATTGCTGGATAAAAATGTTTTAGCAATCTCTGTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAAc          <  2:271031/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCa                                                                                                                                     <  2:307324/125‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTGGATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGGTACCATAATAATGGTGACTGACATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCt      <  2:116336/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gATAAAAATGTTTTAGCAATCTCTTTCTGTCATGAATACATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGACGCGGGAACGCGTGAAGGTGCTGGATGGGATAGCGATTTTCTAGCCAGAATTGGTCAGCAGTTAGGAACGGCTGAAt  >  1:40262/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gATAAAAATGTTTTAGCAATCACTTTCCGTCATGAATCCATGGCTGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATCCCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAAt  <  2:335386/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
AAGGCGGCAACGTCACCATTTTTGACGGCAAGAATAACGTTGTTGACCACGCTGATGGTGATTACGCCGATGAACCGATGATCTACCAGGCGTTTCAACCTCTGCCGCGGTTTGGCGATAGCTACACACTCATCGGTAGCTGGATTGTCGATGATGAAGCGTGCGGAATGGGGATCCGTGAAGATAACACACTGATCACCAAAGACACCTCACGTTTCGTTCCGCATTACATTGCTGGATAAGAATGTTTTAGCAATCTCTTTCTGTCATGAATCCATGGCAGTGACCATACTAATGGTGACTGCCATTGATGGAGGGAGACACAGTGCACTGGCAAACTCACACCGTTTTTAATCAACCTATACCATTAAATAACAGCAACTTATACCTGTCTGATGGCGCGCTCTGCGAAGCGGTAACGCGTGAAGGTGCTGGCTGGGATAGCGATTTTCTAGCCAGTATTGGTCAGCAGTTAGGAACGGCTGAAT  >  NC_000913/4413950‑4414437

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