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,9960GT98.7% 528.5 / ‑2.0 155intergenic (+558/‑301)yiiD/yiiEGNAT family putative N‑acetyltransferase/CopG family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (71/81);  total (72/83)
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.72e-01

Only 100 of 155 total aligned reads displayed.

TAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATAT  >  NC_000913/4078901‑4079091
                                                                                               |                                                                                               
tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTa                                                                                             <  1:1885329/100‑1 (MQ=255)
 aGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATAttca                                                                                            <  1:2847124/100‑3 (MQ=255)
 aGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTaa                                                                                            <  1:272069/100‑1 (MQ=255)
 aGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTaa                                                                                            >  1:407241/1‑100 (MQ=255)
    aaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttt                                                                                         >  1:3142939/1‑100 (MQ=255)
     aaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTc                                                                                        >  1:2632972/1‑100 (MQ=255)
       acaacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAc                                                                                      >  1:4036461/1‑100 (MQ=255)
        caacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                     <  1:1567770/100‑1 (MQ=255)
           cTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTa                                                                                  >  2:1770973/1‑100 (MQ=255)
           cTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTa                                                                                  <  2:1176005/100‑1 (MQ=255)
              tCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAAt                                                                               <  2:2771921/100‑1 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                             >  1:1731442/1‑100 (MQ=255)
                 cGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCt                                                                            >  1:1627144/1‑100 (MQ=255)
                   aaCTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTgg                                                                          <  2:4022440/100‑1 (MQ=255)
                   aaCTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGCTCTGTGATGAAATATTAATTTCACCTTAAATCCTgg                                                                          <  2:346871/100‑1 (MQ=255)
                     cTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTg                                                                        >  2:674998/1‑100 (MQ=255)
                      tGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGc                                                                       <  1:1860772/100‑1 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                      >  1:1116968/1‑100 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                      >  1:494812/1‑100 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                      >  1:1739714/1‑100 (MQ=255)
                         tACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGt                                                                    >  2:2790115/1‑100 (MQ=255)
                          aCCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTAATGAAGTATTAATTGCACCTTAAATCCTGGTGCTGtt                                                                   >  2:2828610/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                  >  1:3994686/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                  <  2:2611742/100‑1 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                  <  1:2281664/100‑1 (MQ=255)
                            caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                 >  1:3488292/1‑100 (MQ=255)
                             aacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTggagatcgga                                                                  >  2:3993201/1‑90 (MQ=255)
                             aacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAgg                                                                >  2:2404437/1‑100 (MQ=255)
                               caaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGAt                                                              >  1:1678114/1‑100 (MQ=255)
                                aaaaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATg                                                             <  1:2684908/100‑1 (MQ=255)
                                aaaaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATg                                                             >  1:3199968/1‑100 (MQ=255)
                                  aaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCt                                                           <  2:3022024/100‑1 (MQ=255)
                                   aTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTg                                                          <  1:1545986/100‑1 (MQ=255)
                                      ataACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACg                                                       <  2:3649881/100‑1 (MQ=255)
                                      ataACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATGAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACg                                                       >  2:2947206/1‑100 (MQ=255)
                                       taACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGt                                                      >  2:574301/1‑100 (MQ=255)
                                       taACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGt                                                      <  1:195307/100‑1 (MQ=255)
                                         aCGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGttt                                                    <  1:1372009/100‑1 (MQ=255)
                                         aCGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGttt                                                    >  2:715206/1‑100 (MQ=255)
                                          cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTCATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACg                                                       >  2:3183801/1‑96 (MQ=255)
                                          cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtttt                                                   <  2:1145087/100‑1 (MQ=255)
                                          cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtttt                                                   >  2:1579420/1‑100 (MQ=255)
                                              gagtatgTAAATACTGGGTGGTGGCTGTATTTACTGGTCTGTGATGAAATATTAATTTTACCTTAGTTCCTGTTGTTTTTAGGATGCTGACGtttt                                                   <  2:2657990/90‑1 (MQ=255)
                                              gAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCgg                                               <  2:1010491/100‑1 (MQ=255)
                                               aGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                              <  2:3339631/100‑1 (MQ=255)
                                                gCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCt                                             <  1:2350217/100‑1 (MQ=255)
                                                gCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCt                                             >  2:180603/1‑100 (MQ=255)
                                                gCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGACTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCt                                             >  2:3338004/1‑100 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            <  1:359411/100‑1 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            >  2:3049351/1‑100 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            <  1:3324660/100‑1 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGGCTGTGATGAAATATTAATTTCACCTTAAATCCTCGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            <  1:6646/100‑1 (MQ=255)
                                                  aaGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTac                                           >  1:1399053/1‑100 (MQ=255)
                                                   aGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTaca                                          <  1:574302/100‑1 (MQ=255)
                                                   aGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTaca                                          >  1:3695892/1‑100 (MQ=255)
                                                    gTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTacac                                         <  1:4104230/100‑1 (MQ=255)
                                                    gTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTacac                                         <  1:3909921/100‑1 (MQ=255)
                                                     tAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACAcc                                        >  1:2208902/1‑100 (MQ=255)
                                                     tAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTAAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTaccc                                         >  2:2681673/1‑97 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTc                                    >  2:3477247/1‑99 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                   >  1:1458805/1‑100 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                   <  1:2488076/100‑1 (MQ=255)
                                                           ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                  <  2:1514697/100‑1 (MQ=255)
                                                           ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                  >  2:3915200/1‑100 (MQ=255)
                                                            cGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGa                                 <  2:3277309/100‑1 (MQ=255)
                                                            cGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGa                                 <  1:1604524/100‑1 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                <  1:3818406/100‑1 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                <  1:3048820/100‑1 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                               <  1:674999/100‑1 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                               <  1:2751810/100‑1 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                <  1:4089122/99‑1 (MQ=255)
                                                                 ggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATg                            <  1:1453903/100‑1 (MQ=255)
                                                                 ggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATg                            >  2:2969320/1‑100 (MQ=255)
                                                                  gtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTAt                             >  1:796246/1‑98 (MQ=255)
                                                                   tggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCa                          <  1:2249520/100‑1 (MQ=255)
                                                                     gCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaa                        >  2:2076439/1‑100 (MQ=255)
                                                                      cAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaat                       <  1:2336356/100‑1 (MQ=255)
                                                                       aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                      <  1:2486197/100‑1 (MQ=255)
                                                                         tACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaa                    <  2:3199954/100‑1 (MQ=255)
                                                                          aCTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                   >  1:2967165/1‑100 (MQ=255)
                                                                           ctgaaaggTCAGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTAACGTTTTTCGTCTACACCATTCAGATTATGCaaacaaat                   <  2:3969573/93‑1 (MQ=255)
                                                                           cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                  >  2:2719443/1‑100 (MQ=255)
                                                                           cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                  <  1:694447/100‑1 (MQ=255)
                                                                           cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                  >  1:3002362/1‑100 (MQ=255)
                                                                            tGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAatat                 >  1:2265112/1‑100 (MQ=255)
                                                                            tGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAatat                 >  2:1423562/1‑100 (MQ=255)
                                                                            tGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAatat                 <  1:2016453/100‑1 (MQ=255)
                                                                             gACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                >  2:162486/1‑100 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg            <  1:1319400/100‑1 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg            <  2:1286614/100‑1 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg            <  1:2176558/100‑1 (MQ=255)
                                                                                    cTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaat         >  2:2125110/1‑100 (MQ=255)
                                                                                    cTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaat         <  1:3492118/100‑1 (MQ=255)
                                                                                      gtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatat       >  2:783489/1‑100 (MQ=255)
                                                                                       tgatgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatata      <  1:3701660/100‑1 (MQ=255)
                                                                                          tgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataata   <  2:529375/100‑1 (MQ=255)
                                                                                          tgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat    <  1:888130/99‑1 (MQ=255)
                                                                                           gaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataatat  >  2:3184643/1‑100 (MQ=255)
                                                                                           gaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataatat  <  2:165622/100‑1 (MQ=255)
                                                                                           gaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataatat  >  1:416366/1‑100 (MQ=255)
                                                                                               |                                                                                               
TAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATAT  >  NC_000913/4078901‑4079091

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