Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,937,547 T→A 5.4% I27N (ATT→AAT)  fucK → L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,937,5470TA 5.4% 239.5 / 3.1 74I27N (ATT→AAT) fucKL‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base T (41/29);  new base A (2/2);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCG  >  NC_000913/2937400‑2937691
                                                                                                                                                    |                                                                                                                                                
cGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAAttg                                                                                                                                                <  1:1379612/150‑1 (MQ=255)
 gggCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAAttgt                                                                                                                                               <  2:851051/150‑1 (MQ=255)
   gCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAAttgttg                                                                                                                                             <  2:606501/150‑1 (MQ=255)
    ccTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGc                                                                                                                                            >  2:1309011/1‑150 (MQ=255)
        aaCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCTACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCcgc                                                                                                                                        <  2:702424/150‑1 (MQ=255)
         aCAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAACAATTTTGccctgc                                                                                                                                       <  2:310511/150‑4 (MQ=255)
           aagcaagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGcc                                                                                                                                     <  2:1085644/150‑1 (MQ=255)
            agcaagAGTGGTTAGCCGGATAAGCTATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGTAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCGCGGGCAACAATTTTGccctgccg                                                                                                                                     <  1:162343/150‑6 (MQ=255)
            agcaagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGTAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCt                                                                                                                                    >  2:1056911/1‑150 (MQ=255)
                agagTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAAc                                                                                                                                >  1:1208652/1‑150 (MQ=255)
                agagTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAAc                                                                                                                                >  2:606276/1‑150 (MQ=255)
                   gTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGCTGCCCGCGCCTCAACGcc                                                                                                                             <  2:623837/150‑1 (MQ=255)
                      gTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTaa                                                                                                                          >  1:606839/1‑150 (MQ=255)
                      gTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTaa                                                                                                                          >  1:1086308/1‑150 (MQ=255)
                       ttAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAAt                                                                                                                         >  1:380514/1‑150 (MQ=255)
                          gCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGcc                                                                                                                      >  1:415593/1‑150 (MQ=255)
                            cGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTAAACGCCTAATGCCAg                                                                                                                    <  2:918484/150‑1 (MQ=255)
                             ggATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGc                                                                                                                   <  1:1381692/150‑1 (MQ=255)
                              gATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCg                                                                                                                  <  1:964480/150‑1 (MQ=255)
                                  aGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGatat                                                                                                              >  1:750060/1‑150 (MQ=255)
                                    cAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATcg                                                                                                            >  2:645339/1‑150 (MQ=255)
                                                  tataTTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAaacaac                                                                                              <  2:556810/150‑1 (MQ=255)
                                                        gCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTgg                                                                                        >  1:1144902/1‑150 (MQ=255)
                                                         cAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGc                                                                                       <  1:1309011/150‑1 (MQ=255)
                                                           ggAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCAc                                                                                     >  2:265953/1‑150 (MQ=255)
                                                               cGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGt                                                                                 <  2:1367830/150‑1 (MQ=255)
                                                                   tATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGtt                                                                             <  1:967922/150‑2 (MQ=255)
                                                                   tATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATAGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGta                                                                             >  2:1234244/1‑149 (MQ=255)
                                                                     tGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCtt                                                                           >  1:309462/1‑150 (MQ=255)
                                                                      gAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCttt                                                                          >  2:454613/1‑150 (MQ=255)
                                                                      gAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCttt                                                                          <  1:952780/150‑1 (MQ=255)
                                                                      gAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCttt                                                                          <  2:309462/150‑1 (MQ=255)
                                                                       aaaCAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTa                                                                         >  1:1152944/1‑150 (MQ=255)
                                                                       aaaCAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTa                                                                         >  1:11969/1‑150 (MQ=255)
                                                                       aaaCAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTa                                                                         <  2:380514/150‑1 (MQ=255)
                                                                        aaCAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAg                                                                        <  2:1239313/150‑1 (MQ=255)
                                                                          cAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGAc                                                                      >  1:472407/1‑150 (MQ=255)
                                                                          cAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGAc                                                                      <  2:415593/150‑1 (MQ=255)
                                                                           aagaagTTAACCTGGTACTCGACTGTCGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGAAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACg                                                                     >  2:198464/1‑150 (MQ=255)
                                                                            agaagTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGc                                                                    >  2:229619/1‑150 (MQ=255)
                                                                                    tCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGc                                                            <  2:1208652/150‑1 (MQ=255)
                                                                                    tCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGAACCAGTGGTCTTTAGACGCCATTTTga                                                            <  1:557712/150‑2 (MQ=255)
                                                                                       tGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACCCCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAAc                                                         >  1:224425/1‑150 (MQ=255)
                                                                                        ggTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACg                                                        >  1:538956/1‑150 (MQ=255)
                                                                                        ggTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACg                                                        >  1:62710/1‑150 (MQ=255)
                                                                                           aCTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCtt                                                     >  1:443115/1‑150 (MQ=255)
                                                                                             tCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTg                                                   <  1:467378/150‑1 (MQ=255)
                                                                                              cGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGc                                                  >  2:1236213/1‑150 (MQ=255)
                                                                                                aCTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTg                                                <  1:1056911/150‑1 (MQ=255)
                                                                                                aCTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTg                                                <  1:1390798/150‑1 (MQ=255)
                                                                                                         cgACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTc                                       <  2:472407/150‑1 (MQ=255)
                                                                                                         cgACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTc                                       >  2:1348276/1‑150 (MQ=255)
                                                                                                            ccAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGc                                    >  1:715840/1‑150 (MQ=255)
                                                                                                               aTGTCAGGGCCATCGCGGTTAAACGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGAAGGAAAACAACACATGGCACCAGTGGTATTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCaaa                                 >  2:963748/1‑150 (MQ=255)
                                                                                                                tGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAAt                                >  2:1217247/1‑150 (MQ=255)
                                                                                                                    aGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAAt                            >  1:1159121/1‑150 (MQ=255)
                                                                                                                     gggCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATa                           <  1:645339/150‑1 (MQ=255)
                                                                                                                     gggCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATa                           <  2:632405/150‑1 (MQ=255)
                                                                                                                      ggCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAg                          <  2:1144902/150‑1 (MQ=255)
                                                                                                                         cATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGa                       >  1:677910/1‑150 (MQ=255)
                                                                                                                         cATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGa                       >  2:863196/1‑150 (MQ=255)
                                                                                                                         cATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGa                       <  1:265953/150‑1 (MQ=255)
                                                                                                                            cgcgGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAact                    >  1:123870/1‑150 (MQ=255)
                                                                                                                             gcgGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAactg                   >  1:210312/1‑150 (MQ=255)
                                                                                                                              cgGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAactga                  >  1:1361935/1‑150 (MQ=255)
                                                                                                                                gTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAactgact                >  2:828540/1‑150 (MQ=255)
                                                                                                                                gTTAATCGGC‑AGGGCAAAATTGTTGCACGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAactgact                >  2:1263550/1‑150 (MQ=255)
                                                                                                                                  tAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAActgactga              >  1:315376/1‑150 (MQ=255)
                                                                                                                                        gC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCca        <  2:677910/150‑1 (MQ=255)
                                                                                                                                         cgcgggCAACAATTTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGccc        >  1:1180264/4‑149 (MQ=255)
                                                                                                                                         cgcgggCAACAATTTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCca        >  1:1100940/4‑150 (MQ=255)
                                                                                                                                         c‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCcac       <  1:606276/150‑1 (MQ=255)
                                                                                                                                            gggCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACAt     >  1:679141/1‑150 (MQ=255)
                                                                                                                                             ggCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATc    <  1:1217247/150‑1 (MQ=255)
                                                                                                                                               cAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGcc         <  2:1052490/143‑1 (MQ=255)
                                                                                                                                               cAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGcc         >  1:1052490/1‑143 (MQ=255)
                                                                                                                                               cAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCcg  >  1:227565/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGC‑AGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCG  >  NC_000913/2937400‑2937691

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