Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,078,996 G→T 100% 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,9960GT99.4% 575.2 / ‑4.5 164intergenic (+558/‑301)yiiD/yiiEGNAT family putative N‑acetyltransferase/CopG family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (79/84);  total (80/84)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 165 total aligned reads displayed.

TTTAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATATA  >  NC_000913/4078899‑4079092
                                                                                                 |                                                                                                
tttAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAatat                                                                                                <  1:3330064/100‑1 (MQ=255)
  tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATAtt                                                                                               <  2:403533/99‑1 (MQ=255)
  tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTa                                                                                              >  1:3986092/1‑100 (MQ=255)
  tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTa                                                                                              <  2:2337382/100‑1 (MQ=255)
   aaaaaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATAtt                                                                                               <  2:451211/96‑1 (MQ=255)
    gAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAt                                                                                            >  1:1526905/1‑100 (MQ=255)
     aaaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAtt                                                                                           <  2:549016/100‑1 (MQ=255)
     aaaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAtt                                                                                           >  2:3791717/1‑100 (MQ=255)
      aaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttt                                                                                          >  2:380717/1‑100 (MQ=255)
        aacaacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCa                                                                                        <  1:2928163/100‑1 (MQ=255)
         acaacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAc                                                                                       >  2:1220283/1‑100 (MQ=255)
          caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                      <  2:4235061/100‑1 (MQ=255)
          caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                      <  2:1887153/100‑1 (MQ=255)
           aacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCt                                                                                     <  2:804199/100‑1 (MQ=255)
           aacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCt                                                                                     >  1:1623358/1‑100 (MQ=255)
            acTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCtt                                                                                    >  1:3904190/1‑100 (MQ=255)
            acTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCtt                                                                                    >  2:2044926/1‑100 (MQ=255)
                tCCCGAACTGTTACCAACAAAATTATAACGATGATCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAAt                                                                                <  2:293824/100‑1 (MQ=255)
                  ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  1:2517327/1‑100 (MQ=255)
                  ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  2:874517/1‑100 (MQ=255)
                  ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                              >  1:3634462/1‑100 (MQ=255)
                     aaCTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTgg                                                                           <  2:1979736/100‑1 (MQ=255)
                       cTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTg                                                                         <  2:1623089/100‑1 (MQ=255)
                         gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                       >  2:732633/1‑100 (MQ=255)
                          ttACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTg                                                                      <  2:2258969/100‑1 (MQ=255)
                           tACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGt                                                                     >  2:4094693/1‑100 (MQ=255)
                            aCCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGtt                                                                    >  1:1430624/1‑100 (MQ=255)
                             ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                   <  2:2521821/100‑1 (MQ=255)
                             ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTAACCGGTCTGTGATGAAATATTAATTTCACCTGAAATCGTGGTGCTGTTa                                                                   >  1:3113625/1‑100 (MQ=255)
                              caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                   >  1:3422266/1‑99 (MQ=255)
                              caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                  <  1:298213/100‑1 (MQ=255)
                              caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                  >  2:2084354/1‑100 (MQ=255)
                              caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGATGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTACTGTTAg                                                                  >  1:91416/1‑100 (MQ=255)
                               aacaaAATTATAACGATGAGCAAGTAAAGACCTGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAgg                                                                 <  2:4060698/100‑1 (MQ=255)
                               aacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAgg                                                                 <  2:4358426/100‑1 (MQ=255)
                               aacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGATCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTATTa                                                                   >  1:2622183/1‑98 (MQ=255)
                                 caaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGAt                                                               >  2:966486/1‑100 (MQ=255)
                                      ttATAACGATGAGCGAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGa                                                          <  1:795040/100‑1 (MQ=255)
                                       tataACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGAc                                                         <  2:3986119/100‑1 (MQ=255)
                                         taACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGt                                                       >  1:3271731/1‑100 (MQ=255)
                                            cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtttt                                                    <  1:506445/100‑1 (MQ=255)
                                               tGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCg                                                 <  1:2044918/100‑1 (MQ=255)
                                                gAGCAAGTAAAGCCCGGTTGATGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCg                                                 <  2:2850388/99‑1 (MQ=255)
                                                gAGCAAGTAAAGACCGGTTTGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCgg                                                <  1:1414541/100‑1 (MQ=255)
                                                 aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                               <  2:3562208/100‑1 (MQ=255)
                                                 aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                               <  1:4147966/100‑1 (MQ=255)
                                                 aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                               <  1:1765155/100‑1 (MQ=255)
                                                       tAAAGACCGGTTGGTCGCAGTACTGACCGGTCTGTTATGAAATATTAATTTCACCTTAAATCCTGGTTCTGTTTTTATGCTGACGTTTTTCGGCTaca                                           <  1:3718802/98‑1 (MQ=255)
                                                        aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACAcc                                         >  1:1946691/1‑99 (MQ=255)
                                                         aaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCAt                                       <  2:3414125/100‑1 (MQ=255)
                                                          aGACCGGTTTGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCAtt                                      >  1:4093183/1‑100 (MQ=255)
                                                           gACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTc                                     <  2:1526910/100‑1 (MQ=255)
                                                            aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                    <  1:1309851/100‑1 (MQ=255)
                                                            aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                    <  2:3138366/100‑1 (MQ=255)
                                                            aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                    <  2:1339589/100‑1 (MQ=255)
                                                             ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                   >  2:1598495/1‑100 (MQ=255)
                                                              cGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGa                                  >  1:790211/1‑100 (MQ=255)
                                                               ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                 <  2:1695473/100‑1 (MQ=255)
                                                               ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                 >  2:3665231/1‑100 (MQ=255)
                                                               ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                 >  1:1796459/1‑100 (MQ=255)
                                                                gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                 <  2:2157210/99‑1 (MQ=255)
                                                                 ttGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTa                               <  1:856474/100‑1 (MQ=255)
                                                                  tggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTAt                              <  2:1339798/100‑1 (MQ=255)
                                                                  tggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTAt                              <  2:3154847/100‑1 (MQ=255)
                                                                    gtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGc                            >  2:270782/1‑100 (MQ=255)
                                                                    gtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGc                            <  2:4105766/100‑1 (MQ=255)
                                                                    gtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGc                            >  2:1810892/1‑100 (MQ=255)
                                                                     tggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCa                           >  1:3131748/1‑100 (MQ=255)
                                                                      ggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaa                          >  1:1026966/1‑100 (MQ=255)
                                                                      ggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCa                           >  2:1197289/1‑99 (MQ=255)
                                                                         ggctggctccGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                       <  1:566216/92‑1 (MQ=255)
                                                                         aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTTACGTTTTTCGGCTACACCATTCAGATTATGCaaata                       <  2:3885511/100‑1 (MQ=255)
                                                                         aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTTTTCTTTTTCATTCAGATTATGCaaata                       >  1:1973974/1‑100 (MQ=255)
                                                                         aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                       >  1:1105681/1‑100 (MQ=255)
                                                                         aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                       <  2:2794287/100‑1 (MQ=255)
                                                                          gTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataa                      >  1:2466698/1‑100 (MQ=255)
                                                                           tACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaa                     <  1:380716/100‑1 (MQ=255)
                                                                            aCTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                    >  2:1013577/1‑100 (MQ=255)
                                                                             cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                   <  1:1371998/100‑1 (MQ=255)
                                                                               gACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                 >  2:586903/1‑100 (MQ=255)
                                                                               gACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                 >  2:3827053/1‑100 (MQ=255)
                                                                                aCCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtct                <  2:1843718/100‑1 (MQ=255)
                                                                                 ccGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtctc               <  1:1220278/100‑1 (MQ=255)
                                                                                 ccGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtctc               <  1:1839060/100‑1 (MQ=255)
                                                                                   ngTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             >  2:1821273/2‑100 (MQ=255)
                                                                                   ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             <  2:1061371/100‑1 (MQ=255)
                                                                                   ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             >  2:1259162/1‑100 (MQ=255)
                                                                                   agTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg             <  1:3105998/99‑1 (MQ=255)
                                                                                       tgtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaata         <  2:3239428/100‑1 (MQ=255)
                                                                                       tgtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaata         <  1:2000750/100‑1 (MQ=255)
                                                                                        gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATGAATATCTCTGaatat        >  2:4326183/1‑100 (MQ=255)
                                                                                        gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatat        >  2:2048748/1‑100 (MQ=255)
                                                                                        gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaa           >  1:966783/1‑97 (MQ=255)
                                                                                        gtgATGAAATATTAATTTCACCTTAAATCCTGGGGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGAAAATAAATATCTCTGaata         >  1:3861267/1‑99 (MQ=255)
                                                                                           atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat     >  1:3942587/1‑100 (MQ=255)
                                                                                           atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat     >  1:3519000/1‑100 (MQ=255)
                                                                                           atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat     <  2:3921761/100‑1 (MQ=255)
                                                                                           atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat     >  2:2190255/1‑100 (MQ=255)
                                                                                             gaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataatat   >  1:1672527/1‑100 (MQ=255)
                                                                                              aaaTATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAatataatata  <  1:874513/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TTTAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATATA  >  NC_000913/4078899‑4079092

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