Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,088 G→T 6.4% T32N (ACT→AAT)  ykgJ ← UPF0153 cysteine cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,0880GT 6.4% 220.4 / 3.9 77T32N (ACT→AAT) ykgJUPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base T (2/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATG  >  NC_000913/303944‑304233
                                                                                                                                                |                                                                                                                                                 
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:435870/1‑150 (MQ=255)
                                                                                                                            cAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAAc                  <  2:1185388/150‑1 (MQ=255)
                                                                                                                                   gcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGTGAAATGCACAACCGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACCGCTAAACATCGTAAAGAGAACCGCGGTa           >  1:1149963/2‑150 (MQ=255)
                                                                                                                                   gcTGGCGGTACTATTCCCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCAGAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTa           >  2:1077677/2‑150 (MQ=255)
                                                                                                                                    cTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTAc          >  1:413190/1‑150 (MQ=255)
                                                                                                                                    cTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTAc          >  1:931312/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)
                                                                                                                                            aaTAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATg  <  2:1112177/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATG  >  NC_000913/303944‑304233

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