Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,889,915 A→C 5.2% intergenic (‑169/+36) 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,9150AC 5.2% 248.5 / 4.0 89intergenic (‑169/+36)fadD/yeaYacyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)/Slp family lipoprotein, RpoE‑regulated
Reads supporting (aligned to +/- strand):  ref base A (46/37);  new base C (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 = 7.23e-01

GTGTTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGG  >  NC_000913/1889766‑1890058
                                                                                                                                                     |                                                                                                                                                
gtgtTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGa                                                                                                                                                  >  2:549706/1‑150 (MQ=255)
  gtTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGccc                                                                                                                                                >  1:839589/1‑150 (MQ=255)
     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: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: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                                      <  2:209362/150‑1 (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:964469/150‑1 (MQ=255)
                                                                                                                tatttTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGgccagc                                 >  1:722077/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                
GTGTTCGTCGTCATCGCAACCCCAAATGATATACATGCCGTTAACATAATATTAACTCATCATACCAGCTTGATAATTACCCAACGAAAAGGTTGCGAAGCGCGTCACTATTTATTTTTATCTTTACCGTAAGAATGCAGAAACAGCGGACCAGCCGCTGTTTCTTTTTCTTTAAAAACAAGCAATTATTCAGTTACAACTGTTTGTACTCTCGCGGGGCCGGGATTATACCAGCCCCATCCGCCGTAGCCATAGGGCCAGCCACGACCGCCATAAAACCATGGATCAATCGG  >  NC_000913/1889766‑1890058

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: