Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,137,219 C→T G24R (GGG→AGG)  talC ← putative transaldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,2190CT99.1% 1565.0 / ‑3.4 436G24R (GGG→AGG) talCputative transaldolase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (199/232);  total (201/235)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 451 total aligned reads displayed.

CCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTC  >  NC_000913/4137128‑4137308
                                                                                           |                                                                                         
 ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTgg                                                                                         >  2:2297697/1‑93 (MQ=255)
  cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                                        >  1:3676092/1‑93 (MQ=255)
  cTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGc                                                                                        >  1:6462003/1‑93 (MQ=255)
    catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                                      <  1:3072669/93‑1 (MQ=255)
      tcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATg                                                                                    >  1:6747274/1‑93 (MQ=255)
         ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                                 <  1:4768731/93‑1 (MQ=255)
           aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGaa                                                                               <  2:1540019/93‑1 (MQ=255)
            aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                              <  1:5895812/93‑1 (MQ=255)
            aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                              >  1:8914278/1‑93 (MQ=255)
             tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                             <  2:5044076/93‑1 (MQ=255)
             acGCTTTTTGCAGACGCGGCAGCCCTTCCCAGATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                             <  2:4465667/92‑1 (MQ=255)
              ggCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAACGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                            <  1:9160141/92‑1 (MQ=255)
               gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           <  1:8592889/93‑1 (MQ=255)
               gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           >  2:6989697/1‑93 (MQ=255)
               gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           >  1:6188099/1‑93 (MQ=255)
               gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           <  2:7546727/93‑1 (MQ=255)
               gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           <  1:1856085/93‑1 (MQ=255)
                 tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                                         >  2:3783477/1‑93 (MQ=255)
                 tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                                         <  1:3938560/93‑1 (MQ=255)
                 tttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCg                                                                         >  1:5253485/1‑93 (MQ=255)
                  ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGt                                                                        <  1:9743801/93‑1 (MQ=255)
                      gCAGACGCTGCAGCACTTTCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGCATATGCGTGCCa                                                                    >  1:4456402/1‑93 (MQ=255)
                      gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCa                                                                    >  2:3367142/1‑93 (MQ=255)
                       cAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAg                                                                   >  2:468559/1‑93 (MQ=255)
                         gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                 >  2:5382812/1‑93 (MQ=255)
                         gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                 >  1:6592092/1‑93 (MQ=255)
                          aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                                <  2:3657819/93‑1 (MQ=255)
                           cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCATCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                               <  1:771323/93‑1 (MQ=255)
                            gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                              >  2:9290601/1‑93 (MQ=255)
                              ggCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTGGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCg                                                            <  2:7778660/93‑1 (MQ=255)
                               gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCgc                                                           >  1:1167191/1‑92 (MQ=255)
                               gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGa                                                           <  2:5196589/93‑1 (MQ=255)
                                cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                          <  1:1987837/93‑1 (MQ=255)
                                 agcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACt                                                         <  1:4652485/93‑1 (MQ=255)
                                  gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                                        >  2:5853841/1‑93 (MQ=255)
                                   caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                                       <  2:8600217/93‑1 (MQ=255)
                                    aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCt                                                      >  2:7116033/1‑93 (MQ=255)
                                     cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                                     <  2:4310642/93‑1 (MQ=255)
                                      ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                    >  2:4124192/1‑93 (MQ=255)
                                           atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGt                                               >  2:6277592/1‑93 (MQ=255)
                                           atatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGt                                               <  1:5704792/93‑1 (MQ=255)
                                            tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                              <  2:6119515/93‑1 (MQ=255)
                                              tGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAg                                            <  2:1113035/93‑1 (MQ=255)
                                                gACTCGTTGCTGGCAGCGATAATGCTCGGGGTGGTTGGCACCCTGGCAATGGGGAATAGGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCg                                          <  1:8185799/93‑1 (MQ=255)
                                                gACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCg                                          >  1:6935904/1‑93 (MQ=255)
                                                gACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCg                                          >  2:4107473/1‑93 (MQ=255)
                                                 gcTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCgg                                         <  2:3544084/92‑1 (MQ=255)
                                                 aCTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCgg                                         <  1:9605608/93‑1 (MQ=255)
                                                 aCTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCgg                                         <  2:3849640/93‑1 (MQ=255)
                                                  cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        <  2:3007667/93‑1 (MQ=255)
                                                  cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        >  2:2724320/1‑93 (MQ=255)
                                                  cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        >  1:5534566/1‑93 (MQ=255)
                                                   tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                       >  1:7391685/1‑93 (MQ=255)
                                                   tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                       <  1:8827749/93‑1 (MQ=255)
                                                         cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                                 <  2:6356402/93‑1 (MQ=255)
                                                         cTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGa                                 <  1:1675200/93‑1 (MQ=255)
                                                           ggCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATa                               <  1:3309083/93‑1 (MQ=255)
                                                            gCAGGGATAATGCTCGGGTTGGTTGGCCCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                              <  2:8308141/93‑1 (MQ=255)
                                                            gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                              <  2:1742351/93‑1 (MQ=255)
                                                            gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                              <  2:8250583/93‑1 (MQ=255)
                                                            gCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                              <  1:6964011/93‑1 (MQ=255)
                                                             cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                             <  1:3441616/93‑1 (MQ=255)
                                                             cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                             <  1:9225943/93‑1 (MQ=255)
                                                              aGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                            >  1:5805328/1‑93 (MQ=255)
                                                               gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                           >  1:2331678/1‑93 (MQ=255)
                                                               gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                           <  2:2917921/93‑1 (MQ=255)
                                                               gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                           <  2:1301919/93‑1 (MQ=255)
                                                                cGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGtt                          <  2:8999134/93‑1 (MQ=255)
                                                                cGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGtt                          >  1:4970919/1‑93 (MQ=255)
                                                                 gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc                         <  2:9899876/93‑1 (MQ=255)
                                                                  aTAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTcc                        <  1:8044093/93‑1 (MQ=255)
                                                                   tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa                       >  1:6102731/1‑93 (MQ=255)
                                                                   tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa                       <  1:10016065/93‑1 (MQ=255)
                                                                   tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa                       >  2:7652632/1‑93 (MQ=255)
                                                                    aaTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCat                      >  2:660323/1‑93 (MQ=255)
                                                                     aTGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatg                     >  2:7379725/1‑93 (MQ=255)
                                                                        gtCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgatg                  <  2:6710689/92‑1 (MQ=255)
                                                                          cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                <  1:8460598/93‑1 (MQ=255)
                                                                          cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                <  1:20897/93‑1 (MQ=255)
                                                                          cGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                <  2:9741150/93‑1 (MQ=255)
                                                                            ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              <  2:10035304/93‑1 (MQ=255)
                                                                            ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              >  2:5699186/1‑93 (MQ=255)
                                                                            ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              >  2:2090309/1‑93 (MQ=255)
                                                                              ttggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTACGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt            <  2:87495/90‑1 (MQ=255)
                                                                               tggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtt           >  1:2878967/3‑93 (MQ=255)
                                                                                tgTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt          <  2:4909427/92‑1 (MQ=255)
                                                                                ggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt          <  2:4634533/92‑1 (MQ=255)
                                                                                 gTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt         <  1:4254167/93‑1 (MQ=255)
                                                                                    gTCGCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  2:1809511/93‑1 (MQ=255)
                                                                                    gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  1:6420701/93‑1 (MQ=255)
                                                                                    gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  2:944470/93‑1 (MQ=255)
                                                                                    gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      >  1:9897148/1‑93 (MQ=255)
                                                                                     tccccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCGGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag     <  1:1458904/90‑1 (MQ=255)
                                                                                     tCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag     >  1:353210/1‑93 (MQ=255)
                                                                                     tCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcag     <  1:1393473/93‑1 (MQ=255)
                                                                                      cACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagt    >  1:7230543/1‑93 (MQ=255)
                                                                                       aCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTTCGTTTAcagtt   >  2:6849038/1‑93 (MQ=255)
                                                                                       aCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagtt   >  1:2640329/1‑93 (MQ=255)
                                                                                        ccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttc  >  1:3852413/1‑93 (MQ=255)
                                                                                        ccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttc  >  1:7903238/1‑93 (MQ=255)
                                                                                           |                                                                                         
CCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTC  >  NC_000913/4137128‑4137308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: