Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,172,796 T→C 5.9% G89G (GGT→GGC mdaB → NADPH quinone reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,172,7960TC 5.9% 346.1 / 7.3 102G89G (GGT→GGCmdaBNADPH quinone reductase
Reads supporting (aligned to +/- strand):  ref base T (47/49);  new base C (3/3);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 108 total aligned reads displayed.

GCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGATCAGTTCT  >  NC_000913/3172648‑3172944
                                                                                                                                                    |                                                                                                                                                    
cccaaggtCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGt                                                                                                                                                      <  4:1168669/142‑1 (MQ=255)
               aCAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTAt                                                                                                                                       <  2:799316/149‑1 (MQ=255)
                 aGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGc                                                                                                                                     >  1:737478/1‑149 (MQ=255)
                     aCTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCt                                                                                                                                           >  1:859352/1‑139 (MQ=255)
                              tCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATgg                                                                                                                            <  4:937652/145‑1 (MQ=255)
                              tCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATgg                                                                                                                            >  3:937652/1‑145 (MQ=255)
                                aaaGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTa                                                                                                                       <  1:1009378/148‑1 (MQ=255)
                                aaaGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTa                                                                                                                       >  2:1009378/1‑148 (MQ=255)
                                 aaGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTAcc                                                                                                                     <  2:136970/149‑1 (MQ=255)
                                       aaGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGTCGCTGTATGCCAGCGATGGTCGTACCCGCaaa                                                                                                               <  1:1465510/149‑1 (MQ=255)
                                         gTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGa                                                                                                             >  3:808348/1‑149 (MQ=255)
                                            cAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTc                                                                                                                          <  1:1133868/133‑1 (MQ=255)
                                            cAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTc                                                                                                                          >  2:1133868/1‑133 (MQ=255)
                                             aaaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGc                                                                                                                  <  2:187530/140‑1 (MQ=255)
                                             aaaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGc                                                                                                                  >  1:187530/1‑140 (MQ=255)
                                             aaaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCg                                                                                                         >  4:468081/1‑149 (MQ=255)
                                             aaaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCg                                                                                                         >  3:478042/1‑149 (MQ=255)
                                              aaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTc                                                                                                                          <  4:1205698/131‑1 (MQ=255)
                                              aaaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTc                                                                                                                          >  3:1205698/1‑131 (MQ=255)
                                               aaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGc                                                                                                                                     <  3:775665/119‑1 (MQ=255)
                                               aaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGc                                                                                                                                     >  4:775665/1‑119 (MQ=255)
                                               aaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTc                                                                                                       <  4:808348/149‑1 (MQ=255)
                                               aaCTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTc                                                                                                       >  3:1183659/1‑149 (MQ=255)
                                                   ttCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGaaa                                                                                                   >  4:276525/1‑149 (MQ=255)
                                                   ttCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGaaa                                                                                                   <  3:276525/149‑1 (MQ=255)
                                                         gggCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATAc                                                                                             >  4:493868/1‑149 (MQ=255)
                                                         gggCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATAc                                                                                             <  3:493868/149‑1 (MQ=255)
                                                         gggCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATAc                                                                                             <  4:816263/149‑1 (MQ=255)
                                                          ggCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACg                                                                                            <  2:737478/149‑1 (MQ=255)
                                                           gCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATg                                                                                                                             <  1:7715/115‑1 (MQ=255)
                                                           gCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATg                                                                                                                             >  2:7715/1‑115 (MQ=255)
                                                             tGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATg                                                                                                                             <  3:1095002/113‑1 (MQ=255)
                                                             tGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATg                                                                                                                             >  4:1095002/1‑113 (MQ=255)
                                                                    gtgATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCcggcgg                                                                                  >  4:117896/1‑149 (MQ=255)
                                                                       aTCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGt                                                                                                                                                      <  1:661502/78‑1 (MQ=255)
                                                                       aTCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGt                                                                                                                                                      >  2:661502/1‑78 (MQ=255)
                                                                       aTCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCt                                                                               >  1:903162/1‑149 (MQ=255)
                                                                            gCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTAc                                                                          >  2:1395154/1‑149 (MQ=255)
                                                                              aGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAg                                                                        <  4:1183659/149‑1 (MQ=255)
                                                                                 tGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACgg                                                                                                                                                 >  1:120773/1‑73 (MQ=255)
                                                                                 tGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACgg                                                                                                                                                 <  2:120773/73‑1 (MQ=255)
                                                                                     aGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTc                                                                                                       <  2:943235/111‑1 (MQ=255)
                                                                                          ggtggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGaaaaaa                                                                                                <  3:333298/113‑1 (MQ=255)
                                                                                          ggtggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGaaaaaa                                                                                                >  4:333298/1‑113 (MQ=255)
                                                                                           gtggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATa                                                                                              >  1:1388842/1‑114 (MQ=255)
                                                                                           gtggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATa                                                                                              <  2:1388842/114‑1 (MQ=255)
                                                                                            tggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAatatat                                                          <  3:117896/149‑1 (MQ=255)
                                                                                             ggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACgg                                                                                           >  4:765528/1‑115 (MQ=255)
                                                                                             ggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATg                                                         >  2:755404/1‑149 (MQ=255)
                                                                                             ggATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGACGAAAAAATACgg                                                                                           <  3:765528/115‑1 (MQ=255)
                                                                                              gATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGc                                                        <  2:903162/149‑1 (MQ=255)
                                                                                                tGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCtt                                                      >  4:310195/1‑149 (MQ=255)
                                                                                                 gggCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCttt                                                     <  1:755404/149‑1 (MQ=255)
                                                                                                  ggCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTtc                                                    <  3:468081/149‑1 (MQ=255)
                                                                                                        ccGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGAc                                              >  3:620029/1‑149 (MQ=255)
                                                                                                        ccGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGAc                                              <  4:620029/149‑1 (MQ=255)
                                                                                                            ggACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAata                                                             >  3:1331243/1‑130 (MQ=255)
                                                                                                            ggACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAata                                                             <  4:1331243/130‑1 (MQ=255)
                                                                                                                  tGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTg                                                                              <  4:744844/107‑1 (MQ=255)
                                                                                                                  tGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTg                                                                              >  3:744844/1‑107 (MQ=255)
                                                                                                                      aaaaTACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCaa                                <  3:310195/149‑1 (MQ=255)
                                                                                                                           aCATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATa                                                                                              >  1:1404566/1‑82 (MQ=255)
                                                                                                                           aCATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATa                                                                                              <  2:1404566/82‑1 (MQ=255)
                                                                                                                             aTTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATcc                                                                                                          >  3:1187422/1‑68 (MQ=255)
                                                                                                                             aTTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATcc                                                                                                          <  4:1187422/68‑1 (MQ=255)
                                                                                                                             aTTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAAc                                       >  2:313724/1‑135 (MQ=255)
                                                                                                                             aTTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAAc                                       <  1:313724/135‑1 (MQ=255)
                                                                                                                                gatgatGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAatatat                                                          <  1:808501/113‑1 (MQ=255)
                                                                                                                                gatgatGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAatatat                                                          >  2:808501/1‑113 (MQ=255)
                                                                                                                                gatgatGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCtt                      <  1:1395154/149‑1 (MQ=255)
                                                                                                                                gatgatGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCtt                      >  2:591938/1‑149 (MQ=255)
                                                                                                                                    atgTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAgg                                                                       <  4:406207/96‑1 (MQ=255)
                                                                                                                                    atgTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAatata                                                           >  2:383156/1‑108 (MQ=255)
                                                                                                                                    atgTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAatata                                                           <  1:383156/108‑1 (MQ=255)
                                                                                                                                      gTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATcc                                                                                                          >  1:1051940/1‑59 (MQ=255)
                                                                                                                                      gTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATcc                                                                                                          <  2:1051940/59‑1 (MQ=255)
                                                                                                                                       taatCACCGAAGGTCACGGGACGCTGTATGCCAGAGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACa                                                                         >  3:974100/4‑91 (MQ=255)
                                                                                                                                       tATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCcg                                                                                      >  2:55747/1‑78 (MQ=255)
                                                                                                                                       tATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCcg                                                                                      <  1:55747/78‑1 (MQ=255)
                                                                                                                                       tATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACa                                                                         <  4:974100/91‑1 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              <  3:53970/101‑1 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              >  2:792001/1‑101 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              <  1:1412368/101‑1 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              >  2:1412368/1‑101 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              <  1:792001/101‑1 (MQ=255)
                                                                                                                                        aTTCACCGAAGGCCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAat                                                              >  4:53970/1‑101 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGc                                                                     >  3:839186/1‑93 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGc                                                                     <  4:839186/93‑1 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAata                                                             <  2:1403216/101‑1 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAata                                                             >  1:1403216/1‑101 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTc                     <  4:282446/141‑1 (MQ=255)
                                                                                                                                         ttCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTc                     >  3:282446/1‑141 (MQ=255)
                                                                                                                                            aCCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCAc                   >  1:626389/1‑140 (MQ=255)
                                                                                                                                            aCCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCAc                   <  2:626389/140‑1 (MQ=255)
                                                                                                                                            aCCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGa          >  3:1463166/1‑149 (MQ=255)
                                                                                                                                              cGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGaaaaa            <  2:1179939/145‑1 (MQ=255)
                                                                                                                                              cGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGaaaaa            >  1:1179939/1‑145 (MQ=255)
                                                                                                                                                    tCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGATCa       <  2:1235501/144‑1 (MQ=255)
                                                                                                                                                    tCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGATCa       >  1:1235501/1‑144 (MQ=255)
                                                                                                                                                    tCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGATCAGttct  >  4:603812/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATCGTTCGCGCCGACAGCGACTACGATGTCAAAGCGGAAGTACAAAACTTTCTCTGGGCTGATGTGGTGATCTGGCAGATGCCAGGCTGGTGGATGGGCGCGCCGTGGACAGTGAAAAAATACATTGATGATGTATTCACCGAAGGTCACGGGACGCTGTATGCCAGCGATGGTCGTACCCGCAAAGATCCGTCGAAAAAATACGGTTCCGGCGGCCTGGTACAGGGCAAAAAATATATGCTTTCTCTGACCTGGAACGCACCAATGGAAGCCTTCACCGAAAAAGATCAGTTCT  >  NC_000913/3172648‑3172944

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