Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,012,938 T→A 5.9% R434R (CGA→CGT yqeB ← XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,9380TA 5.9% 476.1 / 15.4 138R434R (CGA→CGTyqeBXdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base T (60/69);  new base A (4/4);  total (64/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 140 total aligned reads displayed.

GTCGTTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGCCCAAGCG  >  NC_000913/3012791‑3013083
                                                                                                                                                   |                                                                                                                                                 
gtcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGt                                                                                                                                                   >  2:256237/1‑148 (MQ=255)
gtcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGt                                                                                                                                                   >  2:66670/1‑148 (MQ=255)
gtcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTc                                                                                                                                                  <  4:684616/149‑1 (MQ=255)
 tcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCTGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTc                                                                                                                                                  <  1:7359/148‑1 (MQ=255)
 tcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCg                                                                                                                                                 >  4:139834/1‑149 (MQ=255)
 tcgtTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCg                                                                                                                                                 >  3:845239/1‑149 (MQ=255)
     tCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGACGggt                                                                                                                                              >  2:394842/1‑148 (MQ=255)
        acaaGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGggtggt                                                                                                                                           >  2:424121/1‑148 (MQ=255)
         caaGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGACGggtgatc                                                                                                                                          <  2:494786/148‑4 (MQ=255)
          aaGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCGTGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATg                                                                                                                                        >  4:210280/1‑149 (MQ=255)
          aaGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCGGTCGGGTGGTAt                                                                                                                                         >  1:754922/1‑148 (MQ=255)
           aGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGACGggtgatccg                                                                                                                                        <  1:394842/148‑6 (MQ=255)
            gCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGcc                                                                                                                                      <  4:411329/149‑1 (MQ=255)
            gCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCATGCAGCAGGAGCACTGATCACCCGTCGGTTGGTATGc                                                                                                                                       <  2:784169/148‑1 (MQ=255)
             ccACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGccc                                                                                                                                     <  3:826383/149‑1 (MQ=255)
             ccACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGccc                                                                                                                                     <  4:587242/149‑1 (MQ=255)
                      aTCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATAtt                                                                                                                            <  4:416364/149‑1 (MQ=255)
                        ccccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTg                                                                                                                           <  2:744061/148‑1 (MQ=255)
                        ccccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGc                                                                                                                          >  4:589576/1‑149 (MQ=255)
                        ccccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGc                                                                                                                          >  4:43591/1‑149 (MQ=255)
                          ccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGcc                                                                                                                         >  1:474023/1‑148 (MQ=255)
                          ccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCa                                                                                                                        <  3:646523/149‑1 (MQ=255)
                          ccGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGATGGTATGCCCCATAATATTGCCa                                                                                                                        >  4:564318/1‑149 (MQ=255)
                           cGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAg                                                                                                                       <  4:468819/149‑1 (MQ=255)
                             tCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGa                                                                                                                     <  3:28903/149‑1 (MQ=255)
                                   gTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACAcc                                                                                                                <  2:524291/148‑1 (MQ=255)
                                   gTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCg                                                                                                               <  4:524664/149‑1 (MQ=255)
                                    tGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCgg                                                                                                              <  3:1265/149‑1 (MQ=255)
                                     gCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACAcc                                                                                                                >  2:846963/1‑146 (MQ=255)
                                     gCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACAcc                                                                                                                <  1:846963/146‑1 (MQ=255)
                                      cTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGACGggtgatc                                                                                                                                          <  4:412122/119‑4 (MQ=255)
                                      cTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGACGggtgatc                                                                                                                                          >  3:412122/1‑116 (MQ=255)
                                       tttAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTAt                                                                                                           >  4:774723/1‑149 (MQ=255)
                                       tttAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTAt                                                                                                           >  3:846344/1‑149 (MQ=255)
                                        ttAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTAt                                                                                                           >  1:697637/1‑148 (MQ=255)
                                        ttAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTAt                                                                                                           >  1:105645/1‑148 (MQ=255)
                                         tAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAaca                                                                                                                  <  2:347967/140‑1 (MQ=255)
                                         tAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAaca                                                                                                                  >  1:347967/1‑140 (MQ=255)
                                           aTTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATtct                                                                                                        <  2:145091/148‑1 (MQ=255)
                                           aTTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATtct                                                                                                        <  1:871587/148‑1 (MQ=255)
                                            tttCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGggtggt                                                                                                                                           >  3:458298/1‑112 (MQ=255)
                                             ttCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGggtgg                                                                                                                                            <  2:55955/110‑1 (MQ=255)
                                             ttCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGggtgg                                                                                                                                            >  1:55955/1‑110 (MQ=255)
                                             ttCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCt                                                                                                     <  3:43591/149‑1 (MQ=255)
                                             ttCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGATGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCt                                                                                                     <  3:564318/149‑1 (MQ=255)
                                              tCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCt                                                                                                     <  1:256237/148‑1 (MQ=255)
                                              tCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCt                                                                                                     >  1:807033/1‑148 (MQ=255)
                                               cATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTgc                                                                                                   >  4:43938/1‑149 (MQ=255)
                                                  gCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAg                                                                                                                       >  3:891765/1‑126 (MQ=255)
                                                  gCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTgcgc                                                                                                 >  1:166211/1‑148 (MQ=255)
                                                    tCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGACCACCCGACggg                                                                                                                                               >  4:764768/1‑100 (MQ=255)
                                                    tCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGACCACCCGACggg                                                                                                                                               <  3:764768/100‑1 (MQ=255)
                                                      aCCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCa                                                                                                                        >  4:697355/1‑121 (MQ=255)
                                                      aCCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCa                                                                                                                        <  3:697355/121‑1 (MQ=255)
                                                      aCCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAAc                                                                                            >  4:309475/1‑149 (MQ=255)
                                                         aaTCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCa                                                                                          >  2:392038/1‑148 (MQ=255)
                                                          aTCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCaga                                                                                        <  3:309475/149‑1 (MQ=255)
                                                          aTCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCaga                                                                                        <  4:872745/149‑1 (MQ=255)
                                                              aGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTaa                                                                                    >  3:319519/1‑149 (MQ=255)
                                                              aGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTaa                                                                                    >  3:719434/1‑149 (MQ=255)
                                                               ggCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTaa                                                                                    >  1:815618/1‑148 (MQ=255)
                                                                gCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAAt                                                                                  <  3:589576/149‑1 (MQ=255)
                                                                gCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAAt                                                                                  <  3:493778/149‑1 (MQ=255)
                                                                    tCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCAc                                                                               <  1:424121/148‑1 (MQ=255)
                                                                       caTCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGa                                                                           <  3:210280/149‑1 (MQ=255)
                                                                        aTCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGAc                                                                          >  3:520547/1‑149 (MQ=255)
                                                                                   tCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGc                                                               <  4:845239/149‑1 (MQ=255)
                                                                                      ccAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGc                                                                                                                                       >  1:591957/1‑74 (MQ=255)
                                                                                      ccAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGccc                                                             >  2:340933/1‑148 (MQ=255)
                                                                                         gATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCcgcg                                                          <  2:729223/148‑1 (MQ=255)
                                                                                         gATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGa                                                         <  3:550668/149‑1 (MQ=255)
                                                                                          aTCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGa                                                         >  2:473740/1‑148 (MQ=255)
                                                                                          aTCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGAt                                                        <  3:366789/149‑1 (MQ=255)
                                                                                            cGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGAtt                                                       <  1:891957/148‑1 (MQ=255)
                                                                                               cTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGt                                                    >  1:773529/1‑148 (MQ=255)
                                                                                                   tttCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCaa                                               >  4:476943/1‑149 (MQ=255)
                                                                                                   tttCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCaa                                               <  4:283313/149‑1 (MQ=255)
                                                                                                   tttCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCaa                                               <  3:774723/149‑1 (MQ=255)
                                                                                                   tttCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCa                                                <  2:11918/148‑1 (MQ=255)
                                                                                                        cGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGCTCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAAGTAc                                           >  1:506411/1‑148 (MQ=255)
                                                                                                         gTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCg                                         >  3:146574/1‑149 (MQ=255)
                                                                                                              aTCGCATAATGCCTGCAGCAGGAGCACGGATCGCCCGTCGGTTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGa                                    <  4:542380/149‑1 (MQ=255)
                                                                                                              aTCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGa                                    >  3:217455/1‑149 (MQ=255)
                                                                                                                cGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACa                                  >  4:135518/1‑149 (MQ=255)
                                                                                                                      aTGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGt                                                    >  4:799167/1‑125 (MQ=255)
                                                                                                                      aTGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCtt                            >  4:496127/1‑149 (MQ=255)
                                                                                                                      aTGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCt                             >  2:260922/1‑148 (MQ=255)
                                                                                                                      aTGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTCGt                                                    <  3:799167/125‑1 (MQ=255)
                                                                                                                       tGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTCTCAATTACCGCATGACAATCCtt                            >  1:16635/1‑148 (MQ=255)
                                                                                                                          cTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTccc                         <  1:409665/148‑1 (MQ=255)
                                                                                                                          cTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCca                        <  4:217455/149‑1 (MQ=255)
                                                                                                                           tGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCca                        >  2:549090/1‑148 (MQ=255)
                                                                                                                                cagGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGt                   <  1:473740/148‑1 (MQ=255)
                                                                                                                                  gaaGCACGGATCACCCGGCGAGTGGTATGCCCAAACGTATTACCAAGAACACCCGTATTCTCCTGCTCACAACCAGAGTAAAGCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCATCCCCGCAGTaa                 <  1:912998/146‑1 (MQ=255)
                                                                                                                                   gAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCcagc                      <  2:217267/142‑1 (MQ=255)
                                                                                                                                     gCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGc              <  2:649594/148‑1 (MQ=255)
                                                                                                                                      cACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTaa                 <  2:748488/144‑1 (MQ=255)
                                                                                                                                         ataccACCCGACGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGc         >  4:221726/5‑149 (MQ=255)
                                                                                                                                              aCCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCTGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGc              <  3:616358/139‑1 (MQ=255)
                                                                                                                                                ccGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGCCCAAgcg  >  3:217877/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTCGTTCAACAAGCCACGCACCATCCCCGTCAACGGTGCTTTAATTTCATGCTCACCAATCCAGGCAATCACATCGCCCTCTTTCACCAGATCGCCTAATTTCACGTTGGATCGCATAATGCCTGCAGCAGGAGCACGGATCACCCGTCGGGTGGTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGCCCAAGCG  >  NC_000913/3012791‑3013083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: