Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,889,919 G→T 5.4% intergenic (‑173/+32) fadD ← / ← yeaY acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)/Slp family lipoprotein, RpoE‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,889,9190GT 5.4% 292.9 / 5.7 93intergenic (‑173/+32)fadD/yeaYacyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)/Slp family lipoprotein, RpoE‑regulated
Reads supporting (aligned to +/- strand):  ref base G (49/39);  new base T (3/2);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

CGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGC  >  NC_000913/1889771‑1890068
                                                                                                                                                    |                                                                                                                                                      
cttcgtcATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                       <  1:45268/148‑1 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                       >  2:90633/1‑150 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                       >  2:68169/1‑150 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                       >  2:964469/1‑150 (MQ=255)
   ggtcATCGAAACCCCAAATGATATCCATGCCGGTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCTAAGCGCGGCACTATTTATTTTTATCTTTACCGTATGAATGCAGAAACAGCGGACCAGCCGc                                                                                                                                                    <  2:925585/149‑1 (MQ=255)
     tcATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTg                                                                                                                                                  >  1:962268/1‑150 (MQ=255)
         cGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttc                                                                                                                                              >  2:471081/1‑150 (MQ=255)
          gCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGCTGGTCCGCTGtttct                                                                                                                                             <  1:1239385/150‑1 (MQ=255)
          gCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttct                                                                                                                                             >  2:844275/1‑150 (MQ=255)
           cAACCCCAAATGATATACATGCCGTTAACATAATATCAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGCTGGTCCGCTGtttctg                                                                                                                                            <  1:709609/150‑2 (MQ=255)
            aaCCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttcttt                                                                                                                                           <  2:962268/150‑1 (MQ=255)
              ccccAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttcttttc                                                                                                                                         >  2:484186/1‑149 (MQ=255)
                ccAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttctttttct                                                                                                                                       >  1:209362/1‑150 (MQ=255)
                 cAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGATAAGGTTGCGAAGCGAGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGtttctttttctt                                                                                                                                      <  2:339390/150‑1 (MQ=255)
                    aTGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTaa                                                                                                                                  <  1:1151128/150‑1 (MQ=255)
                        tataCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAAc                                                                                                                              <  1:1219221/150‑1 (MQ=255)
                        tataCATGCCGTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACa                                                                                                                             <  2:1420946/150‑1 (MQ=255)
                         ataCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACa                                                                                                                             <  2:1449082/150‑1 (MQ=255)
                           aCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAg                                                                                                                           >  2:666627/1‑150 (MQ=255)
                            cATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGc                                                                                                                          <  2:642262/150‑1 (MQ=255)
                                 cGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTATACCGTAAGAATGAAGAAACAGTGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAatta                                                                                                                     >  2:577836/1‑150 (MQ=255)
                                  gTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAattat                                                                                                                    <  1:1647318/150‑1 (MQ=255)
                                       cATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGt                                                                                                               >  2:948162/1‑150 (MQ=255)
                                              ttAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACt                                                                                                        <  1:1566822/150‑1 (MQ=255)
                                              ttAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACt                                                                                                        <  2:118241/150‑1 (MQ=255)
                                                 aCTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGtt                                                                                                     <  2:222250/150‑1 (MQ=255)
                                                    catcatACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGt                                                                                                  >  2:501935/1‑150 (MQ=255)
                                                      tcatACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTAc                                                                                                >  2:1218243/1‑150 (MQ=255)
                                                       catACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTaa                                                                                                                                  >  2:1075668/1‑115 (MQ=255)
                                                       catACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTaa                                                                                                                                  <  1:1075668/115‑1 (MQ=255)
                                                             aGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGACGCTGTTTCTTTTTATTTAAAAACAAGCAATTATTAAGTTACAACTGTTTGTAATCTcgcc                                                                                         >  2:1545979/1‑149 (MQ=255)
                                                             aGCATGTTATTTACCCATAGAAAAGGTTGAGAAGCGCGTCACTATTTATTATTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTcgcg                                                                                         <  1:1028720/150‑1 (MQ=255)
                                                              gCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCgg                                                                                        >  2:227198/1‑150 (MQ=255)
                                                               cTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCggg                                                                                       >  1:455462/1‑150 (MQ=255)
                                                                 tGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACAGTTTGTACTCTCGCGGGGc                                                                                     >  2:1437548/1‑150 (MQ=255)
                                                                   aTAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCg                                                                                   <  1:68169/150‑1 (MQ=255)
                                                                     acttcaccAACGAAAATGCTGCGAAGCGCGTCACTATTTATTTTTCTCATTACAGTAAGAATGCAGAAACCGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAATCAAGCAAATATTCAGTTACAACTGTTTGTACTCTCTCGGGGCCggg                                                                                 <  2:1478962/144‑1 (MQ=255)
                                                                     aaTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCggg                                                                                 <  1:1218243/150‑1 (MQ=255)
                                                                     aaTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCggg                                                                                 >  1:1446360/1‑150 (MQ=255)
                                                                      aTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGa                                                                                <  1:922257/150‑1 (MQ=255)
                                                                       ttACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGAt                                                                               <  1:1169625/150‑1 (MQ=255)
                                                                           ccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATtata                                                                           >  2:936475/1‑150 (MQ=255)
                                                                           ccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTAACTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGATTCTTTTTCTTTAAAAACAACCAATTATTCAGAAACAACAGTTTGTACTCTCGCGGGGCCGGGATtata                                                                           >  1:616279/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACTACTGTTTGTAGTCTCGCGGGGCCGGGATTATAc                                                                          >  2:278263/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAc                                                                          >  1:1542438/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTATTATTATTTAAAAACAATCAATTAATCAGTTACAACTGTAAGTACTCTCCCGGGGCAGGGATTATAc                                                                          >  1:1352430/1‑150 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                         <  2:339424/150‑1 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                         <  1:1352523/150‑1 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                         <  1:844275/150‑1 (MQ=255)
                                                                                 aaaaGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGcc                                                                     <  1:514865/150‑1 (MQ=255)
                                                                                    aGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCa                                                                  >  2:73894/1‑150 (MQ=255)
                                                                                     ggTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCAt                                                                 >  1:1637573/1‑150 (MQ=255)
                                                                                      gTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATc                                                                <  2:455462/150‑1 (MQ=255)
                                                                                      gTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTCTACTCTCGCGGGGCCGGGATTATACCAGCACCATc                                                                >  2:319860/1‑150 (MQ=255)
                                                                                           gAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATccgccg                                                           >  2:979089/1‑150 (MQ=255)
                                                                                                gcgTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTAAACAAGCCCCATCCGCCGTAGcc                                                      >  1:895829/1‑150 (MQ=255)
                                                                                                 cgTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTCGCCa                                                     >  2:186199/1‑150 (MQ=255)
                                                                                                  gTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCTGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCAt                                                    >  2:643797/1‑150 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                                >  2:1633474/1‑150 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                                <  1:41963/150‑1 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                                <  2:209362/150‑1 (MQ=255)
                                                                                                       tatttattTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgc                                               <  1:501935/150‑1 (MQ=255)
                                                                                                          ttatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccag                                            >  2:155834/1‑150 (MQ=255)
                                                                                                          ttatttTTATCTTTACCCTAAGTATGGATTAACATCGGACCAGCCGGTGGTTGTTTTTCTTTAAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgcca                                             <  1:969164/150‑1 (MQ=255)
                                                                                                           tatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagc                                           >  1:722077/1‑150 (MQ=255)
                                                                                                           tatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagc                                           <  1:964469/150‑1 (MQ=255)
                                                                                                            atttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGATTGTACTCTCGCGGGGCCGGGATTATACCAGCCACATCCGCCGTAGCCATAGTgccagcc                                          >  2:1098195/1‑150 (MQ=255)
                                                                                                                   tCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCg                                   <  1:666627/150‑1 (MQ=255)
                                                                                                                    cTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGAAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGc                                  >  2:226017/1‑150 (MQ=255)
                                                                                                                       tACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGACATAGGGCCAGCCACGACCGCCAt                               >  2:876347/1‑150 (MQ=255)
                                                                                                                        aCCGTAAGAATGCAGAAACAGCGGTCCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCTGTTACAACTGTTTGGACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATa                              <  2:895871/150‑1 (MQ=255)
                                                                                                                        aCCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATa                              <  1:90633/150‑1 (MQ=255)
                                                                                                                         ccGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATaa                             >  2:1718513/1‑150 (MQ=255)
                                                                                                                         ccGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATaa                             >  1:703648/1‑150 (MQ=255)
                                                                                                                              aGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTACTTTACAACTGTTTGTAATCTCGCGGGGCCGGGATTATAACAGCCCCATCCGCCGTAGCCATATGGCCAGCCACGCCCGACATAAAACCa                        >  1:805996/1‑150 (MQ=255)
                                                                                                                               agaaaaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCAt                       >  1:596584/7‑150 (MQ=255)
                                                                                                                                aaTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATg                      >  2:461494/1‑150 (MQ=255)
                                                                                                                                 aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACg                                       >  1:616388/1‑132 (MQ=255)
                                                                                                                                 aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACg                                       <  2:616388/132‑1 (MQ=255)
                                                                                                                                  aaaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGa                    >  1:1412606/4‑150 (MQ=255)
                                                                                                                                   gCAGAATCAGCGGACCAGCAGCGGTTTCTTTTTCTTTAAAAACAAGAAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGAt                   <  1:152203/150‑1 (MQ=255)
                                                                                                                                   gCAGAAACAGCGGACCAGCCGCTGTTTCTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                  <  1:484186/150‑1 (MQ=255)
                                                                                                                                   aaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGAt                   >  1:893983/3‑150 (MQ=255)
                                                                                                                                    cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCTTCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                  >  1:90157/1‑150 (MQ=255)
                                                                                                                                    cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                  <  1:479034/150‑1 (MQ=255)
                                                                                                                                         aCAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCg             <  2:870745/150‑1 (MQ=255)
                                                                                                                                         aCAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGATACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCg             >  2:432115/1‑150 (MQ=255)
                                                                                                                                          cAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg            >  1:989119/1‑150 (MQ=255)
                                                                                                                                          cAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg            >  2:524726/1‑150 (MQ=255)
                                                                                                                                                 ccAGCCGCTGTTTCTTTTTCTTTAAGAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggc     >  2:929524/1‑150 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg   >  1:908991/1‑150 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg   <  2:942956/150‑1 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg   >  1:1549179/1‑150 (MQ=255)
                                                                                                                                                    gCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcggc  <  1:1464906/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGGC  >  NC_000913/1889771‑1890068

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