Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 304,085 A→C 6.5% I33S (ATT→AGT)  ykgJ ← UPF0153 cysteine cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913304,0850AC 6.5% 214.1 / 5.7 77I33S (ATT→AGT) ykgJUPF0153 cysteine cluster protein
Reads supporting (aligned to +/- strand):  ref base A (39/32);  new base C (2/3);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

CTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATG  >  NC_000913/303936‑304233
                                                                                                                                                     |                                                                                                                                                    
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: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)
                                                                                                                                                     |                                                                                                                                                    
CTCTGCATGTGGACGATCGATTTTTATATATCGTGCAGCAGGCATTTTTGCCCGGGGTTCCTGCAAGGGCAATACATCGGGGGTTTTTCTGATTGGTACCGCTCATGCATCGGTGAAAAGGGGATATTTGTTCAGTGAGCCTGGCGGGAATAGTACCGCCAGCATCGTCGGCTTCGGCCCAGTAAAAAGAGACGCGGAAAAATGCACAACAGGCACCACACGTCATGCATGGATTCAGATTGCTCATAATTCACCAGTACAGCTAAAAATCGTAAAGAGAACAGCGGTACATCGTATG  >  NC_000913/303936‑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: ‑