Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,733,600 A→G 7.3% E467G (GAA→GGA)  xylG → xylose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,733,6000AG 7.3% 137.7 / 3.3 63E467G (GAA→GGA) xylGxylose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (34/22);  new base G (2/3);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-02

TGTTACTTAACCCGCGCATTCTCATTCTTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCT  >  NC_000913/3733459‑3733740
                                                                                                                                             |                                                                                                                                            
tgtTACTTAACCCGCGCATTCTCATTCTTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGaa                                                                                                                                             <  1:336193/143‑1 (MQ=255)
     cTTAACCCGCGCATTCTCATTCTTGATGAACCCACCAGGGGTATCGCTATTGGCGCGAACTACGCGATCTACAAATTAATGAACCAACTCGTCCAGCAGGGTATTGCCGTTATTCTCCTCTCTTCCGAATTCCGTGACGTgac                                                                                                                                        >  2:518215/1‑141 (MQ=255)
             gcgcATTCTCATTCTGGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTa                                                                                                                                <  1:90733/143‑1 (MQ=255)
                cattctcattctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACGCGGACAGCAGGGTATTGCCGTGATTTTGCTGTCTTCCGAATAACCTGAAGTGCTCGGCACTAGca                                                                                                                             >  2:1060024/1‑142 (MQ=255)
                 attctcattctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCGGTTGTTGTCATCTTTTCCGAATTACCTAAAGTACTCGATCTTAGcta                                                                                                                            >  1:77402/1‑141 (MQ=255)
                   tctcattctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATc                                                                                                                          <  1:1063374/143‑1 (MQ=255)
                     tcattctTGATGAACCCACCTGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCgt                                                                                                                        <  1:682904/143‑1 (MQ=255)
                     tcattctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTAAAAATTAATTAACCAACTCGTCCAGCAGGGGATTGCCGTTGTTGTCATCTCTTCAGAATTAACTGAAGTGAGCGGACAGAGCGATAgt                                                                                                                        >  1:249802/1‑143 (MQ=255)
                      cattctTGATGAACCCACGAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCGGTTATTGTAATGTATTAAGAATTAACTGAAGTGATGGGCATTAGTGATCgtg                                                                                                                       >  2:388612/1‑143 (MQ=255)
                       attctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGGATTACCTGAAGTGCTCGGCCTTAGCGATCgcgt                                                                                                                      >  1:432777/1‑143 (MQ=255)
                        gtctTGACGAACGCACCAGTGGTATCGATATTGGCGCGAACTACGCGCTCTACAAATTAATTACCCACCTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGTAGTGCTCGGCCTTAGCGATCGTGTa                                                                                                                     <  2:1050505/142‑1 (MQ=255)
                         tctTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTACAGCAGGGGATTGCCGTTCTTGTGATCTCTTCCGAATTCCGTGAAAGGCACGGCCTTGGCGGTCggatac                                                                                                                    >  1:659415/1‑138 (MQ=255)
                            tgatgaACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTACTTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGCACTACCTCAAGTGCTCGCCATTACCCATCgtctcccgg                                                                                                                 >  2:385149/1‑136 (MQ=255)
                             gatgaGCCCTCCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAGCTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGt                                                                                                                <  2:386778/143‑1 (MQ=255)
                               tgaACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTGTTGCCGTTATTGTCCTCTCTTCCGTATTACCTGTAGTGCTCGGCCTTAGCGATCGTGGACTGGTgt                                                                                                              >  1:624376/1‑142 (MQ=255)
                                  aCCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCACCTCGTCCAGCAGGGTATTGCCGTTCTTGTCATCTCTTCCGAATTCCCTGAGGTGCTCGGCCTCAGCCACCATCTACTGGTGAtcc                                                                                                           >  1:544675/1‑141 (MQ=255)
                                     caccaGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGCCTTACCTGAAGTGCACGGCCTTAGCGATCGTGTACTGGTGATGCATg                                                                                                        >  2:366723/1‑143 (MQ=255)
                                             gTATCGATATTGGCGCGAAATACGGTATCTACAAATTAATTAACCATCTCTTCCAGCGGGGTATTGCCGTTATTGTCATCTCTTCCGAGTTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGaaa                                                                                                <  2:140713/143‑1 (MQ=255)
                                             gTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGCATTACCTGAAGTGCTCGGCCTTAGCGATCGTGGACTGGTGCTGCATGAAGGGaca                                                                                                >  1:673931/1‑141 (MQ=255)
                                                tCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGGAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTa                                                                                             >  2:923618/1‑143 (MQ=255)
                                                tCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGGAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTa                                                                                             >  2:363727/1‑143 (MQ=255)
                                                 cGATATGGGCGCGAAATACGAGATCTACAAATTAATTAACCAGCTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTaa                                                                                            <  1:627900/143‑1 (MQ=255)
                                                  gATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCGTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTAGTGATGCATGAGGGGAAAcgaca                                                                                           >  1:12654/1‑139 (MQ=255)
                                                   tccttGGCTCTATATTCGAGTTCTACAACTGTATTTACCATCCTGTCCAGCAGGTTATTTCCGTTAGTGTCCTCTCTGCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGTTGATGCATGAAGGGAAACTaaaa                                                                                          <  2:899676/140‑1 (MQ=255)
                                                       tGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGCTGCATGAAGGGACCCTAAAAGCCa                                                                                      >  2:331879/1‑143 (MQ=255)
                                                         gcgcGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGGAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAAc                                                                                    <  1:363727/143‑1 (MQ=255)
                                                         gcgcGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGGAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAAc                                                                                    <  1:923618/143‑1 (MQ=255)
                                                         gcgcGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTGCTGGTGATGGATGAAGGGAAACTAAAAGCCAAc                                                                                    >  1:960633/1‑143 (MQ=255)
                                                          cgcgAAATACGAGATCTACAAATTCATTAACCCACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAAcc                                                                                   >  2:146060/1‑143 (MQ=255)
                                                          cgcgAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGATAACATCTCCGCGTTGGAATTACCTGAAGTGCTCGGCCTTAGGGTTCGTGTAATGGTTATTCATGACGGACAACCAGAAGCCagac                                                                                   >  2:983023/1‑140 (MQ=255)
                                                            cgAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGACGGGAAACTAAAAGCCGACCTg                                                                                 >  2:672594/1‑143 (MQ=255)
                                                             aaaaTACGAGATCCACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATGGTCATCTCTTCCGAATGACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGa                                                                                <  1:207472/142‑1 (MQ=255)
                                                              aaaTACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTCCCTGAAGTGCTCGGCCTTAGCGCTCGTGTACTGGTCATGCCAGAAGGGAAAATAAAAGCCAAACTGAt                                                                               >  2:403351/1‑143 (MQ=255)
                                                               aaTACGATACCTGCATATTTATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTTAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATa                                                                              <  1:142793/143‑1 (MQ=255)
                                                                        tCTACAAACTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAAc                                                                     <  1:605402/143‑1 (MQ=255)
                                                                          tACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCATAGCGATCGTGGACTGATGATGCAGGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCt                                                                   >  2:519108/1‑143 (MQ=255)
                                                                              aattaattAACCCACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAAGTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACt                                                               <  2:247758/143‑1 (MQ=255)
                                                                                  aattaaCCAACTCGTCCAGCAGGGTATTGCCGTTATTGGCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCAGGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAgg                                                           <  2:422855/143‑1 (MQ=255)
                                                                                   attcaCCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGa                                                          <  1:154692/143‑1 (MQ=255)
                                                                                        ccAACTCGTCCAGCAGGGTATTGCCGTTAATGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAgg                                                     >  1:1061287/1‑143 (MQ=255)
                                                                                          atcTCGTCCAGCAGGGTATTGTCGTTATTGTTATCTCTTCCGAATTACAGGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTg                                                   <  1:549818/141‑1 (MQ=255)
                                                                                           aCTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTCGCGATCCTGTAATGGTGAGGAATGAAGGGAAACTAAAAGCCAACCTGATAACTCCTAAACCGTATCAAGAGCaagcga                                                  >  2:683538/1‑138 (MQ=255)
                                                                                                tCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGCTggca                                             >  1:201325/1‑141 (MQ=255)
                                                                                                tCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGaa                                             >  1:185807/1‑143 (MQ=255)
                                                                                                  cagcagGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGc                                           <  1:594806/143‑1 (MQ=255)
                                                                                                         gtagggCCGTTAGTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTg                                    <  2:863526/138‑1 (MQ=255)
                                                                                                            ttGCCGTCATTGTCATCTCTTCCAAATTACCAGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACATGACCCAGGAGCAGGAGATGGAAGCCGCAttaaag                                 >  2:820161/1‑139 (MQ=255)
                                                                                                             tGCCGTTATTGTCATCTCTTCCGAGTTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTgagga                                <  2:391238/143‑1 (MQ=255)
                                                                                                               ccGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGc                              >  1:377199/1‑143 (MQ=255)
                                                                                                               ccGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGc                              >  1:490537/1‑143 (MQ=255)
                                                                                                               ccGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCGTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCCTTGAGgcgc                              >  2:860069/1‑143 (MQ=255)
                                                                                                                 gTTATTGTCATCTTTTCCGGAGTAGCTGGAGTGCGCGGCCTAAGCGATCGTGTACTGGTGATGCATGAGGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGa                            <  2:627300/143‑1 (MQ=255)
                                                                                                                      tGTCCTCCCTTCCGAATTACCTGAAGTGCTCGGCCTGAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATCACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAAcatc                       <  1:783805/143‑1 (MQ=255)
                                                                                                                         tctctcTGCCGAAGTAGCTGAAGTGCTCGGCCTTAGCGAGCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCCGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATg                    <  2:77402/141‑1 (MQ=255)
                                                                                                                          aTCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACCCAGCAGCAGCTGATGCAAACCGCATTGAGGAGCCAACATCATgc                   >  2:260479/1‑142 (MQ=255)
                                                                                                                           tctcTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGACACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCAGTGAGGAGCGAAcatacggta                  >  2:814893/1‑137 (MQ=255)
                                                                                                                             tctTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCCGGAGGAGGAGAGCGAAACCGCATTGAGGAGCGAACATCATGTCGa                >  2:1029085/1‑143 (MQ=255)
                                                                                                                             gctTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGa                <  1:434954/142‑1 (MQ=255)
                                                                                                                               ttCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGTAACTAAAAGCCAACCTGACAACTCATAACCTGACTAAGAAACAGGTGATGGAAGCCGCAAAGATGAGCGAACATCATGTCGaag              >  2:267758/1‑142 (MQ=255)
                                                                                                                                cccGAACTACCCCAAGTCCTCGGCCTTAGCGACCGTGTACTGGCGATGCATGAAGGGACACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGTAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGaaaa             <  1:69100/142‑1 (MQ=255)
                                                                                                                                  cGAATTACGTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCGTGCCGGGAAACTAAAAGCAAACCGGATAAATCATAAACGGACTCAGGCGCCGGCAATGGTAGCCGCATTGAGGAGCGACGCTCATGTGGagacgc           >  2:105588/1‑138 (MQ=255)
                                                                                                                                        aCCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCg     >  1:738056/1‑143 (MQ=255)
                                                                                                                                         ccTGAAGTGGGCGGCCGTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGt    <  1:814166/143‑1 (MQ=255)
                                                                                                                                           tGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGACAAGCAACACGTCt  >  1:883394/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TGTTACTTAACCCGCGCATTCTCATTCTTGATGAACCCACCAGGGGTATCGATATTGGCGCGAAATACGAGATCTACAAATTAATTAACCAACTCGTCCAGCAGGGTATTGCCGTTATTGTCATCTCTTCCGAATTACCTGAAGTGCTCGGCCTTAGCGATCGTGTACTGGTGATGCATGAAGGGAAACTAAAAGCCAACCTGATAAATCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCT  >  NC_000913/3733459‑3733740

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