Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,078,996 G→T intergenic (+558/‑301) yiiD → / → yiiE GNAT family putative N‑acetyltransferase/CopG family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,078,9960GT90.8% 275.6 / 12.7 87intergenic (+558/‑301)yiiD/yiiEGNAT family putative N‑acetyltransferase/CopG family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (2/6);  new base T (39/40);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATATA  >  NC_000913/4078901‑4079092
                                                                                               |                                                                                                
tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTa                                                                                              >  1:1103884/1‑100 (MQ=255)
    aaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttt                                                                                          >  2:1342659/1‑100 (MQ=255)
      aacaacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCa                                                                                        <  2:1508178/100‑1 (MQ=255)
        caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                      <  1:1342653/100‑1 (MQ=255)
        caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                      <  2:1333354/100‑1 (MQ=255)
        caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                      <  2:606682/100‑1 (MQ=255)
              tCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGCCCGGTCTGTGATGAAATATTAATTTCACCTTAAAt                                                                                >  2:191524/1‑100 (MQ=255)
              tCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAAt                                                                                <  2:157676/100‑1 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAAt                                                                                >  2:791006/1‑99 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATc                                                                               >  1:277557/1‑100 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATc                                                                               <  1:895191/100‑1 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATc                                                                               >  2:1318525/1‑100 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  1:983185/1‑100 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  2:1056794/1‑100 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  2:281135/1‑100 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              <  2:983188/100‑1 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  2:2135479/1‑100 (MQ=255)
                  gAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTg                                                                            >  1:1716592/1‑100 (MQ=255)
                    aCTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGt                                                                          <  1:281135/100‑1 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCt                                                                       >  2:550456/1‑100 (MQ=255)
                         tACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGt                                                                     <  2:1716602/100‑1 (MQ=255)
                          aCCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGtt                                                                    >  2:441648/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                   <  1:1721292/100‑1 (MQ=255)
                            caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                  >  2:1121912/1‑100 (MQ=255)
                            caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                  >  2:261070/1‑100 (MQ=255)
                               caaCATTATAACGATGAGCAAGTAAACAACACTTGGTCTCAGTACCGCACGGTCTGTGATCAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAgg                                                                 <  1:2380614/98‑1 (MQ=255)
                               caaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGAt                                                               >  2:282353/1‑100 (MQ=255)
                                aaaaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATg                                                              >  2:530222/1‑100 (MQ=255)
                                 aaaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGc                                                             <  1:737655/100‑1 (MQ=255)
                                   aTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTg                                                           <  1:598149/100‑1 (MQ=255)
                                   aTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTg                                                           <  1:1318519/100‑1 (MQ=255)
                                      ataACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACg                                                        >  2:2129015/1‑100 (MQ=255)
                                      ataACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACg                                                        <  1:1384165/100‑1 (MQ=255)
                                       taACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGt                                                       >  2:1229977/1‑100 (MQ=255)
                                            aTGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTc                                                  >  2:840931/1‑100 (MQ=255)
                                             tGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCg                                                 <  1:1056791/100‑1 (MQ=255)
                                              gAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCgg                                                <  2:66320/100‑1 (MQ=255)
                                               aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                               <  1:648314/100‑1 (MQ=255)
                                               aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                               <  1:540454/100‑1 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                             >  2:842313/1‑100 (MQ=255)
                                                   aGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTaca                                           <  2:1103887/100‑1 (MQ=255)
                                                    gTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTacac                                          <  2:99631/100‑1 (MQ=255)
                                                     tAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACAcc                                         >  1:2463114/1‑100 (MQ=255)
                                                      aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATCTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCa                                        >  2:1912524/1‑100 (MQ=255)
                                                      aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCa                                        >  1:1381792/1‑100 (MQ=255)
                                                       aaGACCGGTTAGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCAt                                       <  2:1292840/100‑1 (MQ=255)
                                                        aGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCAtt                                      >  2:2400708/1‑100 (MQ=255)
                                                         gACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTc                                     >  2:962030/1‑100 (MQ=255)
                                                         gACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTc                                     >  2:1170243/1‑100 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                    <  2:2507605/100‑1 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                    <  2:245414/100‑1 (MQ=255)
                                                           ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                   >  1:836751/1‑100 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                 <  2:459356/100‑1 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                                <  1:8346/100‑1 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                                <  1:550455/100‑1 (MQ=255)
                                                               ttGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTa                               <  1:2378010/100‑1 (MQ=255)
                                                                tggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTAt                              <  2:2597333/100‑1 (MQ=255)
                                                                 ggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATg                             <  1:1404783/100‑1 (MQ=255)
                                                                 ggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATg                             >  1:162478/1‑100 (MQ=255)
                                                                  gtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGc                            >  2:563758/1‑100 (MQ=255)
                                                                   tggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCa                           >  2:1863068/1‑100 (MQ=255)
                                                                    ggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaa                          >  1:675541/1‑100 (MQ=255)
                                                                       aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                       >  2:2398120/1‑100 (MQ=255)
                                                                         tACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaa                     <  2:230044/100‑1 (MQ=255)
                                                                          aCTGACCGGTCTGTGATGAAATGTTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                    <  1:267106/100‑1 (MQ=255)
                                                                          aCTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                    <  1:261071/100‑1 (MQ=255)
                                                                          aCTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                    >  2:1606887/1‑100 (MQ=255)
                                                                           cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                   <  2:836752/100‑1 (MQ=255)
                                                                            tGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAatat                  >  2:2359860/1‑100 (MQ=255)
                                                                             ttccgatCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                 <  2:1360752/94‑1 (MQ=255)
                                                                              aactgtatgtgATGAAAGATTAATGTCACCTTAAATCCTGGTGCTTTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtct                <  1:2449001/93‑1 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             >  1:977453/1‑100 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             >  1:1816102/1‑100 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             <  2:575298/100‑1 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             >  2:332928/1‑100 (MQ=255)
                                                                                  gTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGa            >  2:2481305/1‑100 (MQ=255)
                                                                                  gTCTGTGATGAAATATGAGTTTCACCTTAAAGCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             <  1:2297311/99‑1 (MQ=255)
                                                                                    tcttgatgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaat          <  1:1835867/97‑1 (MQ=255)
                                                                                     tgtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaata         <  1:2400690/100‑1 (MQ=255)
                                                                                      gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtcagatcggaa        >  1:1360746/1‑91 (MQ=255)
                                                                                      gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatat        <  2:277557/100‑1 (MQ=255)
                                                                                      gggatgatAGATTAATTTCACTTTATGTCCTGGTTCTGTTATGATGGTGACGTTTTTCGGCTACACGATTCAGATTATGCGAATAAATATCTCTGa            <  1:2150985/94‑1 (MQ=255)
                                                                                       tgatgaAATATTAATTTCACCTTAAATCCTGGTGCTATTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaataga       >  2:1835877/1‑98 (MQ=255)
                                                                                        gatgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataa      <  1:441647/100‑1 (MQ=255)
                                                                                         atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat     <  1:962027/100‑1 (MQ=255)
                                                                                         atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTAACGTTTTTCGGCTACACCATTCAGAGTATGCAAATAAATATCTCTGaatataa      <  1:1047001/99‑1 (MQ=255)
                                                                                            aaaTATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAatataatata  <  2:675542/100‑1 (MQ=255)
                                                                                               |                                                                                                
TAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATATA  >  NC_000913/4078901‑4079092

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