Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,152 (G)5→6 coding (1641/1659 nt) nrfE → putative cytochrome c‑type biogenesis protein NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,1471.G98.3% 211.9 / ‑3.6 60coding (1636/1659 nt)nrfEputative cytochrome c‑type biogenesis protein NrfE
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (33/26);  total (34/26)
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

CCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTGCCAG  >  NC_000913/4293011‑4293281
                                                                                                                                         |                                                                                                                                      
ccATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCggg                                                                                                                                     <  2:826137/140‑1 (MQ=255)
       tGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCggggggaaga                                                                                                                              >  2:435225/1‑135 (MQ=255)
       tGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCggggggaaga                                                                                                                              >  2:258610/1‑135 (MQ=255)
       tGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCggggggaaga                                                                                                                              <  1:258610/140‑6 (MQ=255)
          tATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAgc                                                                                                                            >  1:226738/1‑140 (MQ=255)
              cGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGa                                                                                                                        <  2:931474/140‑1 (MQ=255)
               ggTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGAt                                                                                                                       <  2:13307/140‑1 (MQ=255)
                gTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGTGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATg                                                                                                                      >  2:781470/1‑140 (MQ=255)
                 taaTGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGa                                                                                                                     >  2:414027/3‑140 (MQ=255)
                         gagaAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAggg                                                                                                             <  2:25125/140‑1 (MQ=255)
                         gagaAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAggg                                                                                                             <  1:322906/140‑1 (MQ=255)
                             aaaCTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTtc                                                                                                         >  1:619066/1‑140 (MQ=255)
                              aaCTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTtct                                                                                                        <  2:1089359/140‑1 (MQ=255)
                              aaCTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTtct                                                                                                        <  2:580337/140‑1 (MQ=255)
                                cTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCa                                                                                                      <  2:344032/140‑1 (MQ=255)
                                   ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCa                                                                                                      <  2:787519/137‑1 (MQ=255)
                                   ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCa                                                                                                      >  1:787519/1‑137 (MQ=255)
                                   ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                                                                                   >  1:494137/1‑140 (MQ=255)
                                   ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                                                                                   >  1:433261/1‑140 (MQ=255)
                                   ggCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGc                                                                                                   <  1:414027/140‑1 (MQ=255)
                                        gATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTActg                                                                                              <  2:333905/140‑1 (MQ=255)
                                              tACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTa                                                                                        <  2:725526/140‑1 (MQ=255)
                                                   ttttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTAc                                                                                   >  2:878757/1‑140 (MQ=255)
                                                   ttttCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTAc                                                                                   >  1:486671/1‑140 (MQ=255)
                                                           tgtatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTg                                                                           >  1:970984/1‑140 (MQ=255)
                                                             tatgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTgcg                                                                         <  2:361808/140‑1 (MQ=255)
                                                               tgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCa                                                                       <  2:880665/140‑1 (MQ=255)
                                                               tgtaCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCa                                                                       >  2:51247/1‑140 (MQ=255)
                                                                 taCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACg                                                                     >  1:1139159/1‑140 (MQ=255)
                                                                 taCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACg                                                                     >  2:84099/1‑140 (MQ=255)
                                                                  aCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGc                                                                    <  2:970984/140‑1 (MQ=255)
                                                                  aCAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGc                                                                    >  1:241019/1‑140 (MQ=255)
                                                                      aGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAg                                                                >  2:997111/1‑140 (MQ=255)
                                                                      aGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAg                                                                <  2:241019/140‑1 (MQ=255)
                                                                      aGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAg                                                                >  1:255830/1‑140 (MQ=255)
                                                                      aGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAg                                                                <  2:1139159/140‑1 (MQ=255)
                                                                       gCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAgg                                                               >  2:1063773/1‑140 (MQ=255)
                                                                             tgCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCa                                                                       <  2:1118877/126‑1 (MQ=255)
                                                                             tgCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCa                                                                       >  1:1118877/1‑126 (MQ=255)
                                                                                  tGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACAcc                                                    <  2:787521/140‑1 (MQ=255)
                                                                                     aTCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTgg                                                 >  1:247520/1‑140 (MQ=255)
                                                                                               gaggaTTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCa                                       >  1:177248/1‑140 (MQ=255)
                                                                                               gaggaTTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCa                                       >  2:563393/1‑140 (MQ=255)
                                                                                                aggaTTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCaa                                      >  2:1142509/1‑140 (MQ=255)
                                                                                                aggaTTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCaa                                      >  1:965955/1‑140 (MQ=255)
                                                                                                    ttgttgATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCg                                  >  2:985578/1‑140 (MQ=255)
                                                                                                       ttgATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGcaa                               >  1:703501/1‑140 (MQ=255)
                                                                                                       ttgATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGcaa                               >  2:321385/1‑140 (MQ=255)
                                                                                                            gatTGCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAacaaca                          >  1:530410/1‑140 (MQ=255)
                                                                                                                gCGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAAcaacagcaa                      <  2:226738/140‑1 (MQ=255)
                                                                                                                 cGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAAcaacagcaac                     <  2:255830/140‑1 (MQ=255)
                                                                                                                 cGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAAcaacagcaac                     <  1:1063773/140‑1 (MQ=255)
                                                                                                                 cGGGCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATAAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAAcaacagcaac                     >  2:994955/1‑140 (MQ=255)
                                                                                                                   ggCGCATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAAcaacagcaacag                   <  1:563393/140‑1 (MQ=255)
                                                                                                                       cATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGt               <  2:177248/140‑1 (MQ=255)
                                                                                                                       cATTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGt               >  2:814668/1‑140 (MQ=255)
                                                                                                                        aTTGTTAAGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGtt              <  1:878757/140‑1 (MQ=255)
                                                                                                                         ttGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTaa            >  2:282340/1‑140 (MQ=255)
                                                                                                                              aaGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAatat        >  1:1152923/1‑140 (MQ=255)
                                                                                                                              aaGCGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAatat        >  1:1093541/1‑140 (MQ=255)
                                                                                                                                 cGGATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTgc     <  1:985578/140‑1 (MQ=255)
                                                                                                                                   gATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTgcca   >  2:832184/1‑140 (MQ=255)
                                                                                                                                   gATGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTgcca   <  2:247520/140‑1 (MQ=255)
                                                                                                                                    aTGGCGGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTgccag  >  2:1061319/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CCATGACTGGTATGCGGTCATGGGGGAGAAAACTGGGCCGGATCGTTACGCTTTTCGTTTGTATGTACAAAGCGGTGTGCGCTGGATCTGGGGGGGAGGATTGTTGATGATTGCGGGCGCATTGTTAAGCGGATGGCGGGGGAAGAAGCGCGATGAATAAAGGGCTTCTCACGCTACTGCTGTTATTTACCTGTTTTGCGCACGCTCAGGTCGTAGACACCTGGCAATTCGCCAATCCGCAACAACAGCAACAGGCGTTAAATATTGCCAG  >  NC_000913/4293011‑4293281

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