Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,889,918 A→G 5.5% intergenic (‑172/+33) 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,9180AG 5.5% 264.9 / 5.0 90intergenic (‑172/+33)fadD/yeaYacyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)/Slp family lipoprotein, RpoE‑regulated
Reads supporting (aligned to +/- strand):  ref base A (47/38);  new base G (3/2);  total (50/40)
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.33e-01

CGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGG  >  NC_000913/1889771‑1890067
                                                                                                                                                   |                                                                                                                                                      
cttcgtcATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                      <  1:45268/148‑1 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                      >  2:90633/1‑150 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                      >  2:964469/1‑150 (MQ=255)
cgtcgtCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGc                                                                                                                                                      >  2:68169/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)
                    aTGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑aa                                                                                                                                 <  1:1151128/150‑1 (MQ=255)
                        tataCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAAc                                                                                                                             <  1:1219221/150‑1 (MQ=255)
                        tataCATGCCG‑TAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACa                                                                                                                            <  2:1420946/150‑1 (MQ=255)
                         ataCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACa                                                                                                                            <  2:1449082/150‑1 (MQ=255)
                           aCATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAg                                                                                                                          >  2:666627/1‑150 (MQ=255)
                            cATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGc                                                                                                                         <  2:642262/150‑1 (MQ=255)
                                 cGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTATACCGTAAGAATGAAGAAACAGTGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAatta                                                                                                                    >  2:577836/1‑150 (MQ=255)
                                  gTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAattat                                                                                                                   <  1:1647318/150‑1 (MQ=255)
                                       cATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGt                                                                                                              >  2:948162/1‑150 (MQ=255)
                                              ttAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACt                                                                                                       <  1:1566822/150‑1 (MQ=255)
                                              ttAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACt                                                                                                       <  2:118241/150‑1 (MQ=255)
                                                 aCTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGtt                                                                                                    <  2:222250/150‑1 (MQ=255)
                                                    catcatACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGt                                                                                                 >  2:501935/1‑150 (MQ=255)
                                                      tcatACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTAc                                                                                               >  2:1218243/1‑150 (MQ=255)
                                                       catACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑aa                                                                                                                                 >  2:1075668/1‑115 (MQ=255)
                                                       catACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑aa                                                                                                                                 <  1:1075668/115‑1 (MQ=255)
                                                             aGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGACGCTGTTTCTTTTTATTT‑AAAAACAAGCAATTATTAAGTTACAACTGTTTGTAATCTcgcc                                                                                        >  2:1545979/1‑149 (MQ=255)
                                                             aGCATGTTATTTACCCATAGAAAAGGTTGAGAAGCGCGTCACTATTTATTATTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTcgcg                                                                                        <  1:1028720/150‑1 (MQ=255)
                                                              gCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCgg                                                                                       >  2:227198/1‑150 (MQ=255)
                                                               cTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCggg                                                                                      >  1:455462/1‑150 (MQ=255)
                                                                 tGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACAGTTTGTACTCTCGCGGGGc                                                                                    >  2:1437548/1‑150 (MQ=255)
                                                                   aTAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCg                                                                                  <  1:68169/150‑1 (MQ=255)
                                                                     acttcaccAACGAAAATGCTGCGAAGCGCGTCACTATTTATTTTTCTCATTACAGTAAGAATGCAGAAACCGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAATCAAGCAAATATTCAGTTACAACTGTTTGTACTCTCTCGGGGCCggg                                                                                <  2:1478962/144‑1 (MQ=255)
                                                                     aaTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCggg                                                                                >  1:1446360/1‑150 (MQ=255)
                                                                     aaTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCggg                                                                                <  1:1218243/150‑1 (MQ=255)
                                                                      aTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGa                                                                               <  1:922257/150‑1 (MQ=255)
                                                                       ttACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGAt                                                                              <  1:1169625/150‑1 (MQ=255)
                                                                           ccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATtata                                                                          >  2:936475/1‑150 (MQ=255)
                                                                           ccAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTAACTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGATTCTTTTTCTTT‑AAAAACAACCAATTATTCAGAAACAACAGTTTGTACTCTCGCGGGGCCGGGATtata                                                                          >  1:616279/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACTACTGTTTGTAGTCTCGCGGGGCCGGGATTATAc                                                                         >  2:278263/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAc                                                                         >  1:1542438/1‑150 (MQ=255)
                                                                            cAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTATTATTATTT‑AAAAACAATCAATTAATCAGTTACAACTGTAAGTACTCTCCCGGGGCAGGGATTATAc                                                                         >  1:1352430/1‑150 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                        <  1:844275/150‑1 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                        <  2:339424/150‑1 (MQ=255)
                                                                             aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                        <  1:1352523/150‑1 (MQ=255)
                                                                                 aaaaGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGcc                                                                    <  1:514865/150‑1 (MQ=255)
                                                                                    aGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCa                                                                 >  2:73894/1‑150 (MQ=255)
                                                                                     ggTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCAt                                                                >  1:1637573/1‑150 (MQ=255)
                                                                                      gTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATc                                                               <  2:455462/150‑1 (MQ=255)
                                                                                      gTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTCTACTCTCGCGGGGCCGGGATTATACCAGCACCATc                                                               >  2:319860/1‑150 (MQ=255)
                                                                                           gAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATccgccg                                                          >  2:979089/1‑150 (MQ=255)
                                                                                                gcgTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTAAACAAGCCCCATCCGCCGTAGcc                                                     >  1:895829/1‑150 (MQ=255)
                                                                                                 cgTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTCGCCa                                                    >  2:186199/1‑150 (MQ=255)
                                                                                                  gTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCTGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCAt                                                   >  2:643797/1‑150 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                               >  2:1633474/1‑150 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                               <  2:209362/150‑1 (MQ=255)
                                                                                                      cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                               <  1:41963/150‑1 (MQ=255)
                                                                                                       tatttattTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgc                                              <  1:501935/150‑1 (MQ=255)
                                                                                                          ttatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccag                                           >  2:155834/1‑150 (MQ=255)
                                                                                                          ttatttTTATCTTTACCCTAAGTATGGATTAACATCGGACCAGCCGGTGGTTGTTTTTCTTTAAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgcca                                            <  1:969164/150‑1 (MQ=255)
                                                                                                           tatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagc                                          <  1:964469/150‑1 (MQ=255)
                                                                                                           tatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagc                                          >  1:722077/1‑150 (MQ=255)
                                                                                                            atttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGATTGTACTCTCGCGGGGCCGGGATTATACCAGCCACATCCGCCGTAGCCATAGTgccagcc                                         >  2:1098195/1‑150 (MQ=255)
                                                                                                                   tCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCg                                  <  1:666627/150‑1 (MQ=255)
                                                                                                                    cTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGAAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGc                                 >  2:226017/1‑150 (MQ=255)
                                                                                                                       tACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGACATAGGGCCAGCCACGACCGCCAt                              >  2:876347/1‑150 (MQ=255)
                                                                                                                        aCCGTAAGAATGCAGAAACAGCGGTCCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCTGTTACAACTGTTTGGACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATa                             <  2:895871/150‑1 (MQ=255)
                                                                                                                        aCCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATa                             <  1:90633/150‑1 (MQ=255)
                                                                                                                         ccGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATaa                            >  1:703648/1‑150 (MQ=255)
                                                                                                                         ccGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATaa                            >  2:1718513/1‑150 (MQ=255)
                                                                                                                              aGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTACTTTACAACTGTTTGTAATCTCGCGGGGCCGGGATTATAACAGCCCCATCCGCCGTAGCCATATGGCCAGCCACGCCCGACATAAAACCa                       >  1:805996/1‑150 (MQ=255)
                                                                                                                               agaaaaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCAt                      >  1:596584/7‑150 (MQ=255)
                                                                                                                                aaTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATg                     >  2:461494/1‑150 (MQ=255)
                                                                                                                                 aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACg                                      >  1:616388/1‑132 (MQ=255)
                                                                                                                                 aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACg                                      <  2:616388/132‑1 (MQ=255)
                                                                                                                                  aaaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGa                   >  1:1412606/4‑150 (MQ=255)
                                                                                                                                   gCAGAATCAGCGGACCAGCAGCGGTTTCTTTTTCTTT‑AAAAACAAGAAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGAt                  <  1:152203/150‑1 (MQ=255)
                                                                                                                                   gCAGAAACAGCGGACCAGCCGCTGTTTC‑TTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                 <  1:484186/150‑1 (MQ=255)
                                                                                                                                   aaaGAAACAGCGGCTGGTCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGAt                  >  1:893983/3‑150 (MQ=255)
                                                                                                                                    cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCTTCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                 >  1:90157/1‑150 (MQ=255)
                                                                                                                                    cAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATc                 <  1:479034/150‑1 (MQ=255)
                                                                                                                                         aCAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCg            <  2:870745/150‑1 (MQ=255)
                                                                                                                                         aCAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGATACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCg            >  2:432115/1‑150 (MQ=255)
                                                                                                                                          cAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg           >  2:524726/1‑150 (MQ=255)
                                                                                                                                          cAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg           >  1:989119/1‑150 (MQ=255)
                                                                                                                                                 ccAGCCGCTGTTTCTTTTTCTTT‑AAGAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggc    >  2:929524/1‑150 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg  >  1:908991/1‑150 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg  <  2:942956/150‑1 (MQ=255)
                                                                                                                                                   aGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggcgg  >  1:1549179/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTT‑AAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGCGG  >  NC_000913/1889771‑1890067

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: ‑