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,9960GT96.7% 316.6 / ‑0.4 93intergenic (+558/‑301)yiiD/yiiEGNAT family putative N‑acetyltransferase/CopG family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (40/48);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGAATATAATAT  >  NC_000913/4078901‑4079091
                                                                                               |                                                                                               
tAGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTa                                                                                             <  1:1820834/100‑1 (MQ=255)
 aGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTaa                                                                                            >  1:322949/1‑100 (MQ=255)
 aGAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTaa                                                                                            >  1:206312/1‑100 (MQ=255)
 aGAAAAACAACTTTCCCGAAATGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTaa                                                                                            >  2:1409785/1‑100 (MQ=255)
  gAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAt                                                                                           <  2:2034932/100‑1 (MQ=255)
  gAAAAACAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAt                                                                                           <  1:153441/100‑1 (MQ=255)
   tccgatctcTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAtt                                                                                          <  2:1448328/92‑1 (MQ=255)
   aaaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAtt                                                                                          <  2:200494/100‑1 (MQ=255)
   aaaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAtt                                                                                          <  2:1748375/100‑1 (MQ=255)
    aaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttt                                                                                         <  2:1074326/100‑1 (MQ=255)
    aaaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttt                                                                                         <  1:2069859/100‑1 (MQ=255)
     aaaCAACTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTc                                                                                        >  1:1781296/1‑100 (MQ=255)
      aacaacTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCa                                                                                       <  1:776832/100‑1 (MQ=255)
        tctcTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAcc                                                                                     <  1:246141/97‑1 (MQ=255)
         ctcTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCt                                                                                    <  1:1801805/98‑1 (MQ=255)
           cTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAAttagatcgga                                                                                  >  1:1448332/1‑92 (MQ=255)
           cTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCAccaga                                                                                  >  2:246142/1‑97 (MQ=255)
           cTTTCCCGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCtag                                                                                  >  2:1801799/1‑98 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATc                                                                              >  1:1831567/1‑100 (MQ=255)
               cccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATc                                                                              >  1:1438618/1‑100 (MQ=255)
               cccGAACTGTTACCAACAAAATTATACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                             >  2:775531/1‑100 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                             <  1:1422536/100‑1 (MQ=255)
                ccGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATcc                                                                             >  2:582936/1‑100 (MQ=255)
                 cGAACTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCt                                                                            <  2:2205573/100‑1 (MQ=255)
                    aCTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGt                                                                         <  1:1487794/100‑1 (MQ=255)
                     cTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCTGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTg                                                                        <  2:1622301/100‑1 (MQ=255)
                     cTGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTg                                                                        <  1:582935/100‑1 (MQ=255)
                      tGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGc                                                                       <  1:186498/100‑1 (MQ=255)
                      tGTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGc                                                                       <  1:1727496/100‑1 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGca                                                                      >  2:1727491/1‑99 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                      >  1:545385/1‑100 (MQ=255)
                       gTTACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCt                                                                      >  1:548420/1‑100 (MQ=255)
                         tACCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGt                                                                    <  2:453077/100‑1 (MQ=255)
                          aCCAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGtt                                                                   >  1:547663/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGTTGTTa                                                                  >  2:575818/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                  >  2:576268/1‑100 (MQ=255)
                           ccAACAAAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTa                                                                  <  2:128433/100‑1 (MQ=255)
                            caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                 >  2:1847862/1‑100 (MQ=255)
                            caacaaAATTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAg                                                                 >  1:150422/1‑100 (MQ=255)
                              acaaAATTATAACTATGAGCAAGTAAAGACCGGTTGGTTGCAGTACTGACCGGTCTGAGATGAAAGATTAATTCCACCTGAAATCCTGGTGCTGTTa                                                                  >  2:1426091/1‑97 (MQ=255)
                                 aaaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGc                                                            >  1:1957130/1‑100 (MQ=255)
                                  aaTTATAACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAACATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCt                                                           <  2:1414694/100‑1 (MQ=255)
                                       taACGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGt                                                      >  2:1697899/1‑100 (MQ=255)
                                        aaCGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtt                                                     <  1:789595/100‑1 (MQ=255)
                                        aaCGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtt                                                     >  1:1348157/1‑100 (MQ=255)
                                         aCGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGttt                                                    >  1:2021427/1‑100 (MQ=255)
                                          cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtttt                                                   <  2:1702795/100‑1 (MQ=255)
                                          cGATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGtttt                                                   <  2:236347/100‑1 (MQ=255)
                                           gATGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGttttt                                                  >  2:1697561/1‑100 (MQ=255)
                                            aTGAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTc                                                 <  1:878839/100‑1 (MQ=255)
                                              gAGCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCgg                                               <  2:206313/100‑1 (MQ=255)
                                               aGCAAGTAAAGACAGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGc                                              >  2:328308/1‑100 (MQ=255)
                                                gCAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCt                                             <  2:33839/100‑1 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            >  1:990327/1‑100 (MQ=255)
                                                 cAAGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTa                                            >  1:1743666/1‑100 (MQ=255)
                                                  aaGTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTac                                           <  1:775533/100‑1 (MQ=255)
                                                    gTAAAGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTacac                                         <  1:1951131/100‑1 (MQ=255)
                                                      aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCa                                       >  1:2128373/1‑100 (MQ=255)
                                                      aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCa                                       <  1:1261138/100‑1 (MQ=255)
                                                      aaaGACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAACATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCa                                       >  2:1261135/1‑100 (MQ=255)
                                                         gACCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTc                                    <  2:1770267/100‑1 (MQ=255)
                                                          tctggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                   <  2:503629/97‑1 (MQ=255)
                                                          aCCGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCa                                   <  2:868275/100‑1 (MQ=255)
                                                           ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                  <  1:260637/100‑1 (MQ=255)
                                                           ccGGTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                  >  2:478290/1‑100 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCaaga                                >  1:503627/1‑97 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAg                                  >  2:794782/1‑98 (MQ=255)
                                                             ggTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAt                                >  2:1099965/1‑100 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                               >  2:1155878/1‑100 (MQ=255)
                                                              gTTGGTGGCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGAtt                               <  2:1438614/100‑1 (MQ=255)
                                                                 ggtggCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATg                            <  1:972084/100‑1 (MQ=255)
                                                                     gCAGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaa                        >  2:348485/1‑100 (MQ=255)
                                                                       gatcttGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                      <  2:632067/95‑1 (MQ=255)
                                                                       aGTACTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaata                      <  2:1348155/100‑1 (MQ=255)
                                                                       aGTACTGACCGGTCTGTGATGAAAGATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACCCCCTTCAGCTTATGCaaat                       >  2:2050873/1‑99 (MQ=255)
                                                                          aCTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataaat                   >  2:1935321/1‑100 (MQ=255)
                                                                           cTGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAaataaata                  <  1:794784/100‑1 (MQ=255)
                                                                            tGACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCaaataagatc                 >  1:632068/1‑96 (MQ=255)
                                                                             gACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                >  2:1501248/1‑100 (MQ=255)
                                                                             gACCGGTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATAtc                <  1:328306/100‑1 (MQ=255)
                                                                                cGGTCTGTGATGAAATATTAATTTAACCCTAAAACCTGGTTCTGTTAGGATGCTGACATTTTTCGGCTACACCATTCATATTATGCAAATAAATATcact             >  2:1374382/1‑100 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg            <  1:791018/100‑1 (MQ=255)
                                                                                 ggTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTg            <  2:547664/100‑1 (MQ=255)
                                                                                  gTCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGa           >  2:2049161/1‑100 (MQ=255)
                                                                                   tCTGTGATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaa          <  1:855772/100‑1 (MQ=255)
                                                                                    ntgtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaat         >  2:1239445/2‑100 (MQ=255)
                                                                                     tgtgATGAAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaata        <  2:2128364/100‑1 (MQ=255)
                                                                                         atgaAATATTAATTTCTTCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACGATTCAGATTAGGCAAATAAATATCTCTGaatataat    <  2:548421/100‑1 (MQ=255)
                                                                                         atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCTTTCAGATTATGCAAATAAATATCTCTGaatataat    <  1:1099967/100‑1 (MQ=255)
                                                                                         atgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataat    <  2:545386/100‑1 (MQ=255)
                                                                                          tgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataata   >  1:2051800/1‑100 (MQ=255)
                                                                                          tgaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataata   <  1:1155881/100‑1 (MQ=255)
                                                                                           gaAATATTAATTTCACCTTAAATCCTGGTGCTGTTAGGATGCTGACGTTTTTCGGCTACACCATTCAGATTATGCAAATAAATATCTCTGaatataatat  <  1:1409788/100‑1 (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: