Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,889,917 C→G 5.6% intergenic (‑171/+34) 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,9170CG 5.6% 292.5 / 6.9 89intergenic (‑171/+34)fadD/yeaYacyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)/Slp family lipoprotein, RpoE‑regulated
Reads supporting (aligned to +/- strand):  ref base C (46/38);  new base G (3/2);  total (49/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 = 4.33e-01

GTTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGC  >  NC_000913/1889768‑1890065
                                                                                                                                                     |                                                                                                                                                     
gtTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGccc                                                                                                                                                       >  1:839589/1‑150 (MQ=255)
   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)
                       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                                                                                                                               <  1:1075668/115‑1 (MQ=255)
                                                          catACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTaa                                                                                                                               >  2:1075668/1‑115 (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:1446360/1‑150 (MQ=255)
                                                                        aaTTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCggg                                                                              <  1:1218243/150‑1 (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:844275/150‑1 (MQ=255)
                                                                                aaCGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATAcc                                                                      <  1:1352523/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:209362/150‑1 (MQ=255)
                                                                                                         cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                             >  2:1633474/1‑150 (MQ=255)
                                                                                                         cTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAggg                                             <  1:41963/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                                    <  2:616388/132‑1 (MQ=255)
                                                                                                                                    aTGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACg                                    >  1:616388/1‑132 (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         >  2:524726/1‑150 (MQ=255)
                                                                                                                                             cAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCgg         >  1:989119/1‑150 (MQ=255)
                                                                                                                                                    ccAGCCGCTGTTTCTTTTTCTTTAAGAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAggc  >  2:929524/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGGCTGAGGC  >  NC_000913/1889768‑1890065

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: