Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 931,428 A→G 5.2% L208P (CTG→CCG)  trxB ← thioredoxin reductase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913931,4280AG 5.2% 244.3 / 2.4 77L208P (CTG→CCG) trxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base G (2/2);  total (40/37)
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.00e+00

CCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAGCGCCTCTTCAAC  >  NC_000913/931286‑931576
                                                                                                                                              |                                                                                                                                                    
cccTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTa                                                                                                                                                 <  4:61366/148‑1 (MQ=255)
        aaaTCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGgtg                                                                                                                                        >  2:350704/1‑149 (MQ=255)
          aTCTCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGtgtg                                                                                                                                      <  3:209670/149‑1 (MQ=255)
              cAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGCACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAg                                                                                                                                   >  5:157533/1‑148 (MQ=255)
               aGTATTCTGGCTGTGACCGATAGCAACAAACAGACCGGAAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGaa                                                                                                                                 <  7:360579/149‑1 (MQ=255)
                   ttCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAatgat                                                                                                                             <  5:45253/149‑1 (MQ=255)
                   ttCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAatga                                                                                                                              <  1:187115/148‑1 (MQ=255)
                   ttCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACGCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGGACGGTTGGTGTGCAGAatgat                                                                                                                             <  3:259111/149‑1 (MQ=255)
                      gggCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGt                                                                                                                           <  2:176449/148‑1 (MQ=255)
                            tgACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGtt                                                                                                                    >  2:196799/1‑149 (MQ=255)
                             gACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTAGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTtc                                                                                                                   >  7:330447/1‑149 (MQ=255)
                                   tAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACt                                                                                                             <  7:306686/149‑1 (MQ=255)
                                    aGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACtt                                                                                                            >  2:390787/1‑149 (MQ=255)
                                       aacaaacaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCCCTTTAt                                                                                                         >  3:272656/1‑149 (MQ=255)
                                       aacaaacaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTAt                                                                                                         >  4:271057/1‑149 (MQ=255)
                                        acaaacaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTGTCCACTTTATc                                                                                                        >  5:19367/1‑149 (MQ=255)
                                        acaaacaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATc                                                                                                        >  1:217587/1‑149 (MQ=255)
                                          aaacaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCa                                                                                                      >  8:52263/1‑149 (MQ=255)
                                            acaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGCCGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCAt                                                                                                     >  5:333390/1‑148 (MQ=255)
                                            acaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGACCGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCAt                                                                                                     >  4:100715/1‑148 (MQ=255)
                                            acaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATc                                                                                                    <  3:251078/149‑1 (MQ=255)
                                            acaGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATc                                                                                                    <  8:221064/149‑1 (MQ=255)
                                              aGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCCCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGAGTACGGTTGGTGTGCAGAATGATGTTGACGTTCTCCACATTATCCATCAg                                                                                                  >  8:149193/1‑149 (MQ=255)
                                               gACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCAATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAgg                                                                                                 >  7:320977/1‑149 (MQ=255)
                                                  cGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGgcgc                                                                                              >  5:215447/1‑149 (MQ=255)
                                                    gCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACg                                                                                                                                               >  4:192796/1‑98 (MQ=255)
                                                    gCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACg                                                                                                                                               <  3:192796/98‑1 (MQ=255)
                                                      aaCGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTaa                                                                                          >  2:346317/1‑149 (MQ=255)
                                                       aCGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGATCGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAAt                                                                                         <  5:353095/149‑1 (MQ=255)
                                                       aCGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAAt                                                                                         <  3:299634/149‑1 (MQ=255)
                                                            gAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACgg                                                                                                                                              >  8:334998/1‑91 (MQ=255)
                                                            gAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACgg                                                                                                                                              <  7:334998/91‑1 (MQ=255)
                                                                gACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGAtttt                                                                                <  8:47586/149‑1 (MQ=255)
                                                                 aCTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGAttttt                                                                               >  5:242386/1‑149 (MQ=255)
                                                                     gATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTc                                                                             <  6:242386/147‑1 (MQ=255)
                                                                      aTGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTACCCATCAGGCGCTTAATGa                                                                                       >  8:274600/1‑136 (MQ=255)
                                                                      aTGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTACCCATCAGGCGCTTAATGa                                                                                       <  7:274600/136‑1 (MQ=255)
                                                                         ttATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTATTGAGGATTTTTTCcgcgcg                                                                       <  6:215447/149‑1 (MQ=255)
                                                                         ttATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCcgcgcg                                                                       >  3:80911/1‑149 (MQ=255)
                                                                           aTCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCgg                                                                      >  5:169238/1‑148 (MQ=255)
                                                                           aTCGCTGGTTTGCGTTTCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCCTCACGCGCTTAATGAGGATTTTTTCCGCGCgg                                                                      <  7:149193/148‑1 (MQ=255)
                                                                            tCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGa                                                                     >  7:246778/1‑148 (MQ=255)
                                                                                  ttttGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCcgc                                                                          <  3:267126/137‑1 (MQ=255)
                                                                                  ttttGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCcgc                                                                          >  4:267126/1‑137 (MQ=255)
                                                                                  ttttGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               >  5:226182/1‑148 (MQ=255)
                                                                                  ttttGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               >  7:58187/1‑148 (MQ=255)
                                                                                    ttGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCa                                                             <  2:217587/148‑1 (MQ=255)
                                                                                              cgcAGACGAACGGCAGTGACACCCTGTTGATCGCCGGTCACTTCTTCCAGCGTTCGGTTGGTGTGCAGAATGAGGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAAt                                                  <  6:169238/149‑1 (MQ=255)
                                                                                              cgcAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAAt                                                  >  7:207706/1‑149 (MQ=255)
                                                                                              cgcAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAAt                                                  >  2:2910/1‑149 (MQ=255)
                                                                                               gcAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATc                                                 <  1:196799/149‑1 (MQ=255)
                                                                                               gcAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATc                                                 <  8:207706/149‑1 (MQ=255)
                                                                                               gcAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATc                                                 <  3:271057/149‑1 (MQ=255)
                                                                                                  gACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGa                                              <  6:226182/149‑1 (MQ=255)
                                                                                                    cGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATg                                            <  1:390787/149‑1 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCGGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               <  3:312171/129‑1 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCGGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               <  3:317475/129‑1 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCGGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               >  4:312171/1‑129 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCGGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGt                                                               >  4:317475/1‑129 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCg                                                                >  8:58845/1‑128 (MQ=255)
                                                                                                     gAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCg                                                                <  7:58845/128‑1 (MQ=255)
                                                                                                         gCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACtt                                       <  1:2910/149‑1 (MQ=255)
                                                                                                                    ccATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGAtttttt                                                                              >  7:156472/1‑99 (MQ=255)
                                                                                                                    ccATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGAtttttt                                                                              <  8:156472/99‑1 (MQ=255)
                                                                                                                    ccATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATg                            >  5:271541/1‑149 (MQ=255)
                                                                                                                     cATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGt                           >  6:43571/1‑149 (MQ=255)
                                                                                                                      aTTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCAt                                                                                                     <  5:107624/74‑1 (MQ=255)
                                                                                                                      aTTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCAt                                                                                                     >  6:107624/1‑74 (MQ=255)
                                                                                                                       tttGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTa                         <  6:157533/149‑1 (MQ=255)
                                                                                                                                   tCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAgc             <  1:346317/149‑1 (MQ=255)
                                                                                                                                   tCACTTCTTCCAGCGTACGGTTGGTGTGCAGAACGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAgc             >  4:120394/1‑149 (MQ=255)
                                                                                                                                       ttcttcCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTa                         >  5:319754/1‑133 (MQ=255)
                                                                                                                                       ttcttcCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGAt                  <  3:249869/140‑1 (MQ=255)
                                                                                                                                       ttcttcCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGAt                  >  4:249869/1‑140 (MQ=255)
                                                                                                                                       ttcttcCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAGCGCct         >  2:58690/1‑149 (MQ=255)
                                                                                                                                             cAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACa                     <  6:60826/131‑1 (MQ=255)
                                                                                                                                             cAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACa                     >  5:60826/1‑131 (MQ=255)
                                                                                                                                              aGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATGTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAGCGCCTCTTCAAc  <  4:272656/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGACCGGCAACGTCGAGTGACTCGATGTTATCGCTGTTTTGCGTATCGCGCAGACGAACGCCAGTGACACCCATTTGATCGCCGGTCACTTCTTCCAGCGTACGGTTGGTGTGCAGAATGATGTTGCCGTTCTCCACTTTATCCATCAGGCGCTTAATGAGGATTTTTTCCGCGCGGAAACCGTCACGGCGGTGAATCAGATGCACTTCCGAAGCGATGTTAGACAGATACAGCGCCTCTTCAAC  >  NC_000913/931286‑931576

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: