Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,555 C→A 6.6% D107Y (GAT→TAT)  gatD ← galactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,5550CA 6.6% 312.8 / 5.0 89D107Y (GAT→TAT) gatDgalactitol‑1‑phosphate dehydrogenase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (33/49);  new base A (5/2);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCACCGCATCAA  >  NC_000913/2172426‑2172680
                                                                                                                                 |                                                                                                                             
gACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTgcg                                                                                                                     >  1:128150/1‑140 (MQ=255)
    aaCGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACt                                                                                                                 <  2:789887/140‑1 (MQ=255)
    aaCGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATACGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACt                                                                                                                 <  2:929986/140‑1 (MQ=255)
    aaCGCTAATCGGCTCAATAAAAGCCCCATCCACAATTGCCAAATCCGTTGGTAGAGAAAAGACATTTTTTCGCTTAACGACAATATATTCTGCAAATCCACCATCTCGCCGCGAGCCAATAAAATCATATTTTGCGCACt                                                                                                                 <  2:1157058/140‑1 (MQ=255)
           aTCGGCTCAATAATAGCCCCATCCGCTATAGGCAGATCCGTGGGGAGAGCAAAGACATTTTTTCCCTTAACGACAATATATTCAGCAAAACCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGAACTGGGAATa                                                                                                          >  2:681693/1‑140 (MQ=255)
           aTCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATa                                                                                                          >  1:73166/1‑140 (MQ=255)
            tCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATaa                                                                                                         <  1:1073720/140‑1 (MQ=255)
              gtctataTAAAAGCAACATCCTAAATAGGCATATCTTTGGGTAGAGCAAACACATTTTTTCGATTAACGACTATATATCCAGCAAATCCTCCCTCCCGCCGCTAGCCAATTAAATCATATTTTGCGCACTGTGAATaaaa                                                                                                       <  1:949680/134‑1 (MQ=255)
                cTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAAcc                                                                                                     >  2:1365281/1‑140 (MQ=255)
                       aaaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAg                                                                                              <  2:971729/140‑1 (MQ=255)
                       aaaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAg                                                                                              <  2:493820/140‑1 (MQ=255)
                       aaaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAg                                                                                              <  2:432785/140‑1 (MQ=255)
                       aaaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAg                                                                                              <  2:310500/140‑1 (MQ=255)
                       aaaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAg                                                                                              <  2:310219/140‑1 (MQ=255)
                        aaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGa                                                                                             <  2:1023679/140‑1 (MQ=255)
                        aaGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGa                                                                                             <  2:1478372/140‑1 (MQ=255)
                                tCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTgg                                                                                     <  1:850399/140‑1 (MQ=255)
                                tCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTgg                                                                                     <  1:1161407/140‑1 (MQ=255)
                                tCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTgg                                                                                     >  2:183814/1‑140 (MQ=255)
                                tCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTgg                                                                                     <  1:850397/140‑1 (MQ=255)
                                tCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTgg                                                                                     <  2:1348021/140‑1 (MQ=255)
                                 ccTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGa                                                                                    <  1:1600006/140‑1 (MQ=255)
                                  ctctATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAANTCCACAATCACGCCGCGAGCCAATCAAATCATATTTTTCGCCCTGGGATTAAAACCCTTTCAGACACTCTGGaa                                                                                   <  2:827881/140‑2 (MQ=255)
                                     aaTAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                <  2:955325/140‑1 (MQ=255)
                                     aaTAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                >  1:798996/1‑140 (MQ=255)
                                     aaTAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                <  2:1401337/140‑1 (MQ=255)
                                        aGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTaa                                                                             <  1:1081057/140‑1 (MQ=255)
                                          gCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTaaaa                                                                           <  1:1552685/140‑1 (MQ=255)
                                          gCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTaaaa                                                                           <  1:812448/140‑1 (MQ=255)
                                          gCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTaaaa                                                                           <  1:812449/140‑1 (MQ=255)
                                          gCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTaaaa                                                                           <  1:847841/140‑1 (MQ=255)
                                           cATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAAc                                                                          <  1:1553047/140‑1 (MQ=255)
                                           cATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAAc                                                                          >  2:1488756/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAATTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACa                                                                         >  1:86504/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCCATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGAAAAGTAAAACa                                                                         >  2:1320023/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACa                                                                         >  1:1440688/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACa                                                                         >  1:1439001/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACa                                                                         >  2:1320339/1‑140 (MQ=255)
                                            atatCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATAAAACCCGTTCAGACACTCTGGACAATTAAAACa                                                                         >  2:948916/1‑140 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAggg                                                                      >  1:678133/1‑140 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAggg                                                                      >  1:961990/1‑140 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATaaaa                                                                                                       <  1:309331/107‑1 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATaaaa                                                                                                       >  2:309264/1‑107 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATaaaa                                                                                                       >  2:309331/1‑107 (MQ=255)
                                               tCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATaaaa                                                                                                       <  1:309264/107‑1 (MQ=255)
                                                     ggTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGtaataa                                                                >  1:284844/1‑140 (MQ=255)
                                                          agCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGca                                                           >  1:1104951/1‑140 (MQ=255)
                                                               aGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAg                                                      >  1:108252/1‑140 (MQ=255)
                                                                gACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAgg                                                     <  2:709396/140‑1 (MQ=255)
                                                                  cATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCa                                                   <  2:499819/140‑1 (MQ=255)
                                                                  cATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCa                                                   <  2:220118/140‑1 (MQ=255)
                                                                    ttttttCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAAc                                                 >  2:1111873/1‑140 (MQ=255)
                                                                      ttttCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCg                                               <  2:642739/140‑1 (MQ=255)
                                                                        ttCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCa                                             <  1:420036/140‑1 (MQ=255)
                                                                        ttCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCa                                             <  1:420054/140‑1 (MQ=255)
                                                                        ttCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTACAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCa                                             <  1:905135/140‑1 (MQ=255)
                                                                           gCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCg                                          >  1:1253942/1‑140 (MQ=255)
                                                                           gCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCg                                          >  1:544921/1‑140 (MQ=255)
                                                                           gCTTAACGACAATAAAATCAGCAAATCCACCATCACGCCGCGAGCCAAAAAAATCATATTATGCGCACTGGGAATAAAACCCATTCAGACACTCAGGACAAGTAAAACAGGTTAATAACGGCACACAGGCAACCGCATCg                                          >  2:1317659/1‑140 (MQ=255)
                                                                            cTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGc                                         <  1:589527/140‑1 (MQ=255)
                                                                              tAACGACAATATATTCAGCAAATCCACCATCTCGCCGCGAGCCAATAAAATCATATTTTGCGCTCTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCa                                       <  2:1559180/140‑1 (MQ=255)
                                                                              tAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCa                                       <  2:665994/140‑1 (MQ=255)
                                                                                      atatatTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTaa                               <  1:1519990/140‑1 (MQ=255)
                                                                                      atatatTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACAGCATCGCCAGGGTGTaa                               <  1:877942/140‑1 (MQ=255)
                                                                                            tCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAatcatc                         >  1:136407/1‑140 (MQ=255)
                                                                                            tCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAatcatc                         >  1:92145/1‑140 (MQ=255)
                                                                                              aGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCaa                       <  1:522901/140‑1 (MQ=255)
                                                                                              aGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCaa                       <  1:523076/140‑1 (MQ=255)
                                                                                              aGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCaa                       <  1:126496/140‑1 (MQ=255)
                                                                                               gCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAac                      >  2:1359176/1‑140 (MQ=255)
                                                                                                cAAATCCACCATCACGCCGCGAGCCAATTAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAaca                     <  2:1370552/140‑1 (MQ=255)
                                                                                                cAAATCCACCATCACGCCGCGAGCCAATAACATCATATTTTGCGCACTGGGAATAAAACCCTTTCTGACACTCGGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAaca                     <  2:428944/140‑1 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                >  1:1684641/1‑78 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                >  1:1684502/1‑78 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                <  2:1684641/78‑1 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                <  2:1684502/78‑1 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                <  2:1682887/78‑1 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAg                                                                                >  1:1682887/1‑78 (MQ=255)
                                                                                                   aTCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACACTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCg                  >  2:776017/1‑140 (MQ=255)
                                                                                                          ctCACGCCGCGAGCCAATAATATCAAATTTAGCGCACTGGGAATAAAACCCTGTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCAACGCCAGGGTGTAAATCAGCAACACCGGATCCCa           >  2:1297760/2‑140 (MQ=255)
                                                                                                          aTCACGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACTGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCa           >  2:1305004/1‑140 (MQ=255)
                                                                                                          aTCAAGCCGCGAGCCAATAAAATAATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCa           >  2:1304424/1‑140 (MQ=255)
                                                                                                           tCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCAc          <  1:1343159/140‑1 (MQ=255)
                                                                                                           tCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCAc          <  2:678133/140‑1 (MQ=255)
                                                                                                           tCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCAc          <  1:537454/140‑1 (MQ=255)
                                                                                                             aCGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCACCg        >  2:914190/1‑140 (MQ=255)
                                                                                                              cgccgcGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCACCGc       >  1:351359/1‑140 (MQ=255)
                                                                                                                 cgcgAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCACCGCATc    >  2:592812/1‑140 (MQ=255)
                                                                                                                   cgAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCCTCGCCAGGGTGTAAATCATCAACACCGGATCCCTCCGCATCag  >  2:282763/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                             
GACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGACATTTTTTCGCTTAACGACAATATATTCAGCAAATCCACCATCACGCCGCGAGCCAATAAAATCATATTTTGCGCACTGGGAATAAAACCCTTTCAGACACTCTGGACAAGTAAAACAGGGTAATAACGGCACACAGGCAACCGCATCGCCAGGGTGTAAATCATCAACACCGGATCCCACCGCATCAA  >  NC_000913/2172426‑2172680

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: