Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,675,597 C→A 6.9% A55A (GCG→GCT yphC ← putative zinc‑binding dehydrogenase YphC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,675,5970CA 6.9% 239.2 / 5.8 72A55A (GCG→GCTyphCputative zinc‑binding dehydrogenase YphC
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base A (3/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TTCGGGCAAAAACCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTGGTAAATAAGCT  >  NC_000913/2675448‑2675744
                                                                                                                                                     |                                                                                                                                                   
ttCGGGCAAAAACCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTgccgcc                                                                                                                                                     <  1:32908/150‑1 (MQ=255)
  cGGGCAAAAACCACAGCCAGAAATGTGATACACCAGCATGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCcgccgc                                                                                                                                                   >  2:1259492/1‑150 (MQ=255)
  cGGGCAAAAACCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCcgccgc                                                                                                                                                   >  1:352526/1‑150 (MQ=255)
      cAAAAACCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCtgtg                                                                                                                                               >  1:244072/1‑150 (MQ=255)
        aaaaCCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGc                                                                                                                                             <  2:40342/150‑1 (MQ=255)
         aaaCCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCa                                                                                                                                            >  2:330476/1‑150 (MQ=255)
            ccACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACgg                                                                                                                                         >  2:1243342/1‑150 (MQ=255)
            ccACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCACAGCGGCGGCAccc                                                                                                                                         <  1:1361779/150‑3 (MQ=255)
            ccACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCACAGCGGCGGCAccc                                                                                                                                         <  2:124499/150‑3 (MQ=255)
               caGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtgt                                                                                                                                      <  2:1082347/150‑1 (MQ=255)
                 gCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtg                                                                                                                                    <  2:352526/150‑1 (MQ=255)
                  ccAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtga                                                                                                                                   >  1:883408/1‑150 (MQ=255)
                   cAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGCAGGCAGCCTTGCCCCATCGCCACCATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGTTTTATCGTGTTCCGCCGCTGTGGCACGGTGTtgat                                                                                                                                  <  2:897588/150‑1 (MQ=255)
                    aGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatg                                                                                                                                 <  1:805013/150‑1 (MQ=255)
                        aTGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatag                                                                                                                             <  1:899622/150‑1 (MQ=255)
                        aTGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatag                                                                                                                             <  2:457434/150‑1 (MQ=255)
                        aTGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatag                                                                                                                             >  2:57382/1‑150 (MQ=255)
                        aTGTGAAACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACATTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatag                                                                                                                             <  2:1012022/150‑1 (MQ=255)
                         tgtgATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgataga                                                                                                                            <  1:783946/150‑1 (MQ=255)
                             aTACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGAtata                                                                                                                        >  2:199491/1‑150 (MQ=255)
                              tACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAg                                                                                                                       <  1:270012/150‑1 (MQ=255)
                              tACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAg                                                                                                                       <  2:958064/150‑1 (MQ=255)
                                cacCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTg                                                                                                                     >  2:1365816/1‑150 (MQ=255)
                                cacCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTg                                                                                                                     >  2:251607/1‑150 (MQ=255)
                                  ccAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGGGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGa                                                                                                                   <  2:1018849/150‑1 (MQ=255)
                                   cAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGAc                                                                                                                  >  2:294747/1‑150 (MQ=255)
                                       aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCg                                                                                                              <  1:1536531/150‑1 (MQ=255)
                                           ggTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTc                                                                                                          <  1:1131038/150‑1 (MQ=255)
                                                cccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCaa                                                                                                     >  2:692194/1‑150 (MQ=255)
                                                  ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAAt                                                                                                   <  1:330476/150‑1 (MQ=255)
                                                   tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATc                                                                                                  <  1:708635/150‑1 (MQ=255)
                                                    ctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTAATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATcc                                                                                                 >  1:893094/1‑150 (MQ=255)
                                                     tttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCATGGTATACCGGTTTAACGGTTCCCCCCGGTTTGGCACGGTGTTGAGGATAGATATAGGGGACTGCGCTGCCCCAAAGccc                                                                                                >  2:511812/1‑150 (MQ=255)
                                                       tAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGCCCGTTGATAAAGCCCTGGGATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgg                                                                                              >  1:188714/1‑150 (MQ=255)
                                                        aaaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgga                                                                                             >  2:854603/1‑150 (MQ=255)
                                                         aaaTGGCGGCAGCCTTGCCCCATCGCCACAATATGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATAGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggag                                                                                            >  2:1466007/1‑150 (MQ=255)
                                                                ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCa                                                                                     <  1:1037141/150‑1 (MQ=255)
                                                                     ccTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTg                                                                                >  1:1313878/1‑150 (MQ=255)
                                                                        tGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATc                                                                             <  1:854603/150‑1 (MQ=255)
                                                                          ccccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAg                                                                           >  1:1397132/1‑150 (MQ=255)
                                                                               tCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCt                                                                      <  1:632236/150‑1 (MQ=255)
                                                                                cGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTg                                                                     <  1:1206820/150‑1 (MQ=255)
                                                                                cGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTg                                                                     <  1:919700/150‑1 (MQ=255)
                                                                                  ccACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGt                                                                   <  1:251607/150‑1 (MQ=255)
                                                                                    acaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTa                                                                 <  2:244072/150‑1 (MQ=255)
                                                                                       aTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATc                                                              <  2:1068976/150‑1 (MQ=255)
                                                                                        tCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATcc                                                             <  1:56939/150‑1 (MQ=255)
                                                                                           gCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCg                                                          >  2:396218/1‑150 (MQ=255)
                                                                                           gCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGCGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCg                                                          >  2:701401/1‑150 (MQ=255)
                                                                                            cccGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTg                                                                     <  2:406751/138‑1 (MQ=255)
                                                                                            cccGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTg                                                                     >  1:406751/1‑138 (MQ=255)
                                                                                                  cGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCgg                                                   <  2:537914/150‑1 (MQ=255)
                                                                                                   ggTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGc                                                  <  1:294747/150‑1 (MQ=255)
                                                                                                    gTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCa                                                 >  1:61410/1‑150 (MQ=255)
                                                                                                      tCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCAcc                                               <  1:309233/150‑1 (MQ=255)
                                                                                                        aTGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGc                                             >  2:513202/1‑150 (MQ=255)
                                                                                                         tGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCa                                            >  1:372366/1‑150 (MQ=255)
                                                                                                            ccGTTGTTAAAGCCCTGGTATAACTGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGAGATCGCTTCCGCAAATCCCTGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCTGCTTCGGCACCGCAACt                                         <  2:16170/150‑1 (MQ=255)
                                                                                                               ttGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGAACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTcc                                      <  1:257535/150‑1 (MQ=255)
                                                                                                                tGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTccc                                     <  1:274291/150‑1 (MQ=255)
                                                                                                                 gATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCg                                    >  1:1306347/1‑150 (MQ=255)
                                                                                                                   tAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCa                                  <  2:188714/150‑1 (MQ=255)
                                                                                                                    aaaGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAg                                 >  2:912788/1‑150 (MQ=255)
                                                                                                                     aaGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGa                                >  1:170318/1‑150 (MQ=255)
                                                                                                                         ccTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGa                            <  1:1259492/150‑1 (MQ=255)
                                                                                                                              tataACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCCGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTc                       <  1:57382/150‑1 (MQ=255)
                                                                                                                               ataACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCg                      >  1:743004/1‑150 (MQ=255)
                                                                                                                               ataACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCg                      >  1:737504/1‑150 (MQ=255)
                                                                                                                                 aaCGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGaa                    >  2:98302/1‑150 (MQ=255)
                                                                                                                                 aaCGGTTTATCGGGTGCAGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGaa                    >  2:412069/1‑150 (MQ=255)
                                                                                                                                       caacaccgTGCCACAGCGGCGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCt              >  2:340065/8‑150 (MQ=255)
                                                                                                                                         acaccgTGCCACAGCGGCGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTgg            >  1:935415/6‑150 (MQ=255)
                                                                                                                                         acaccgTGCCACAGCGGCGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTgg            >  1:479849/6‑150 (MQ=255)
                                                                                                                                            gggTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTGGTaa         >  2:1304065/1‑150 (MQ=255)
                                                                                                                                               tGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTGGTAAATa      >  1:997069/1‑150 (MQ=255)
                                                                                                                                                   gccgcTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTGGTAAATAAGCt  <  2:883408/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCGGGCAAAAACCACAGCCAGAAATGTGATACACCAGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATTTTGATCAGTACCTGGTTAATCCCCGGCGTCGGCACCGCAACTTCCCGCAGATCGACGGTCGAATTTCCTGGTAAATAAGCT  >  NC_000913/2675448‑2675744

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