Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,083 G→T 6.6% P34T (CCC→ACC)  ykgJ ← UPF0153 cysteine cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,0830GT 6.6% 190.9 / 4.0 75P34T (CCC→ACC) ykgJUPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base T (2/3);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

TTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTAT  >  NC_000913/303934‑304232
                                                                                                                                                     |                                                                                                                                                     
ttCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCggg                                                                                                                                                       <  2:1060570/150‑1 (MQ=255)
tgctctGCATGTGGCCGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCggg                                                                                                                                                       <  2:688941/148‑1 (MQ=255)
  ctctGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATAGGGGGTTTTACTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGaa                                                                                                                                                     >  2:918530/1‑150 (MQ=255)
   tctGCTTGTGGACGATCGATTTTTATTTATCGTGCAGCAGGCATTTTTGCCCGGGGGTCCTGCTAGGGCTATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGTATATTTGTTCAGTGAGCCTGGCGGGAAt                                                                                                                                                    <  2:863168/150‑1 (MQ=255)
          gaggACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCACGTGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCg                                                                                                                                             <  1:42331/148‑1 (MQ=255)
           tgGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGGTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGc                                                                                                                                            <  1:1191206/150‑1 (MQ=255)
            ggACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCTGTGAGCCTGGCGGGAATAGTACCGcc                                                                                                                                           >  2:353351/1‑150 (MQ=255)
             gACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCa                                                                                                                                          >  1:1185388/1‑150 (MQ=255)
               cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGTACTATTCCCGCCAgg                                                                                                                                        <  1:964145/150‑2 (MQ=255)
               cgatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAga                                                                                                                                        >  1:33989/1‑149 (MQ=255)
                gatcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGGACCGCTCATGCAGCGGTGAAAAGGGGATATTTGTTCAGTGAG‑CTGGCGGTACTATTCCCGCCAggct                                                                                                                                      <  1:35350/150‑4 (MQ=255)
                 atcgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCCTGCATCGGTGAAAAGGGGATAATTGTTCAGTGAG‑CTGGCGGTACTATTCCCGCCAggctc                                                                                                                                     <  2:152335/150‑5 (MQ=255)
                 atcgatTTTTATAATTCGTTCAGCAGGCATTTTTGCCCGGGGCTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAt                                                                                                                                      <  2:525756/150‑1 (MQ=255)
                   cgatTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCAtcg                                                                                                                                    >  1:1207805/1‑150 (MQ=255)
                            tatatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCg                                                                                                                           >  1:70792/1‑150 (MQ=255)
                              tatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCCTGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGc                                                                                                                         >  2:1166992/1‑150 (MQ=255)
                              tatCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGc                                                                                                                         <  1:898031/150‑1 (MQ=255)
                                tCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGccc                                                                                                                       >  2:342319/1‑150 (MQ=255)
                                 cGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCa                                                                                                                      <  1:1290974/150‑1 (MQ=255)
                                  gTGCAGCAGGCATGTTTTCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTGCGGCCCAg                                                                                                                     <  2:1114262/150‑1 (MQ=255)
                                    gcagcaGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATAGATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTa                                                                                                                   <  2:126717/150‑1 (MQ=255)
                                       gcagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCGGTGAGCCTGGCGGGAATAGTACCGCCCGCATCGTCGGCTTCGGCCCAGTaaaa                                                                                                                >  2:635765/1‑150 (MQ=255)
                                        cagGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTaaaaa                                                                                                               >  1:1139840/1‑150 (MQ=255)
                                          ggCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAaga                                                                                                             <  2:90467/150‑1 (MQ=255)
                                               ttatGCCAGGGGTTCCTGCAAGGGCAATACATCGGCGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAATAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGAcgc                                                                                                        <  2:1064299/147‑1 (MQ=255)
                                                tttGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGAcgcg                                                                                                       >  1:279395/1‑150 (MQ=255)
                                                tttGCCCGGGGGGCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGGGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGAcgcg                                                                                                       <  1:1157965/150‑1 (MQ=255)
                                                 ttGCCCGGGGTGCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCGCATGCAGCGGTGAAAAGGGGATATTTGTGCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCgg                                                                                                      <  2:625861/150‑1 (MQ=255)
                                                  tGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGa                                                                                                     <  1:894309/150‑1 (MQ=255)
                                                    cccGGGGTTCCTGCAAGGGCAATACATCGGGGGGTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTCCCGCCAGCATCGTCGGCTTCGGCCCAGAAAAAAGAGACGCCGaaa                                                                                                   >  2:968364/1‑150 (MQ=255)
                                                         ggTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCGCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCGGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGc                                                                                              <  2:896782/150‑1 (MQ=255)
                                                          gTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGca                                                                                             >  1:238329/1‑150 (MQ=255)
                                                            tCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGcaca                                                                                           <  2:1261833/150‑1 (MQ=255)
                                                              cTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAACAGACGCGGAAAAATGaacaac                                                                                         >  2:227229/1‑150 (MQ=255)
                                                                gCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTTCCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAg                                                                                       <  1:10571/150‑1 (MQ=255)
                                                                  aaGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGc                                                                                     <  2:878356/150‑1 (MQ=255)
                                                                    gggCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCACAGTAAAAAGAGACGCGGATAAATGCACAACAGGcac                                                                                   >  2:536003/1‑150 (MQ=255)
                                                                      gCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGcacca                                                                                 >  1:1110558/1‑150 (MQ=255)
                                                                       cAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGcaccac                                                                                >  1:1274837/1‑150 (MQ=255)
                                                                            cATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTc                                                                           >  2:1275478/1‑150 (MQ=255)
                                                                               cGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTcatg                                                                        >  1:1226427/1‑150 (MQ=255)
                                                                                    gTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCTTgg                                                                   <  1:701414/150‑1 (MQ=255)
                                                                                    gTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATgg                                                                   <  1:222270/150‑1 (MQ=255)
                                                                                    gTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGACAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTcatgcatgt                                                                   >  1:660610/1‑149 (MQ=255)
                                                                                      ttttCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGAt                                                                 >  1:109810/1‑150 (MQ=255)
                                                                                          cTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAg                                                             >  1:1147511/1‑150 (MQ=255)
                                                                                                ggTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCt                                                       >  1:1178728/1‑150 (MQ=255)
                                                                                                 gTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTc                                                      <  2:432663/150‑1 (MQ=255)
                                                                                                  gaCCGCTCATGCATCGGTGAAAAGGGGATATTTGTGCAGGGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAATTGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCa                                                     <  1:888040/149‑1 (MQ=255)
                                                                                                    ccGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATa                                                   >  1:591900/1‑150 (MQ=255)
                                                                                                      gCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAAt                                                 >  2:1157051/1‑150 (MQ=255)
                                                                                                        tCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGG‑‑‑CACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                            <  2:303244/150‑1 (MQ=255)
                                                                                                         cATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCa                                              <  1:298402/150‑1 (MQ=255)
                                                                                                           tGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCAcc                                            >  1:141591/1‑150 (MQ=255)
                                                                                                             cATCGGTGAACAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAg                                          <  2:680511/150‑1 (MQ=255)
                                                                                                             cATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAg                                          <  2:309256/150‑1 (MQ=255)
                                                                                                                cGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAACTGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTAATAATTCACCAGTAc                                       >  2:858072/1‑150 (MQ=255)
                                                                                                                 ggTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACTCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACa                                      <  1:631148/150‑1 (MQ=255)
                                                                                                                  gTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAg                                     >  1:412872/1‑150 (MQ=255)
                                                                                                                     aaaaGGGGATATTTGTTCAGTGAGCCTGGCGTGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTa                                  <  1:1184933/150‑1 (MQ=255)
                                                                                                                       aaGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaa                                <  2:238329/150‑1 (MQ=255)
                                                                                                                        aGGGGATATTTGTTCAGTGCGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTaaaa                               >  1:885759/1‑150 (MQ=255)
                                                                                                                         ggggATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGCCGCGGAAAAATGCACAACAGGCACCACACGGCATGCAAGGATTCAGATTGCTCATAATTCCCCATTACAGCTaaaaa                              >  1:863210/1‑150 (MQ=255)
                                                                                                                           ggATAGTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGATTCGGACCAGTAAAAAGAGCCGCGGAAAAATGCCCAACAGGCACCACAAGTCATGCATGGATTCAGATTGCTCCTAAATCACCACTACAGCTTAAATTc                            >  1:71322/1‑150 (MQ=255)
                                                                                                                            gATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCg                           <  2:367896/150‑1 (MQ=255)
                                                                                                                             atatTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGt                          <  2:33989/150‑1 (MQ=255)
                                                                                                                              tatTTGTTCAGTGAGCCTGGCGGGATTAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTa                         <  1:1166992/150‑1 (MQ=255)
                                                                                                                              tatTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTa                         >  1:1245367/1‑150 (MQ=255)
                                                                                                                                      cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAgagacc                 >  2:502015/1‑148 (MQ=255)
                                                                                                                                      cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAAc                 <  2:1185388/150‑1 (MQ=255)
                                                                                                                                      cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAAc                 >  2:435870/1‑150 (MQ=255)
                                                                                                                                             gcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGTGAAATGCACAACCGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACCGCTAAACATCGTAAAGAGAACCGCGGTa          >  1:1149963/2‑150 (MQ=255)
                                                                                                                                             gcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCAGAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTa          >  2:1077677/2‑150 (MQ=255)
                                                                                                                                              cTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTAc         >  1:931312/1‑150 (MQ=255)
                                                                                                                                              cTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTAc         >  1:413190/1‑150 (MQ=255)
                                                                                                                                                ggCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCCCCACACGTCATGCATGGATTCAGATTGCTCCTAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACAt       >  2:895559/1‑150 (MQ=255)
                                                                                                                                                  cGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTGAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCg     >  1:474203/1‑150 (MQ=255)
                                                                                                                                                     gAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCCCCACACGTCATGCATGGATTCAGATTGCTCATAAATCACCAGTACCGCTAAAAATCGTAAAGAGAACAGGGGTACATCGTAt  >  1:165455/1‑150 (MQ=255)
                                                                                                                                                     gAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCCCCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTAt  >  2:979253/1‑150 (MQ=255)
                                                                                                                                                     gAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTAt  <  2:1207805/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTAT  >  NC_000913/303934‑304232

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