Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,210 A→G 73.7% K287E (AAA→GAA)  amiA → N‑acetylmuramoyl‑l‑alanine amidase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,2100AG73.7% 274.2 / 101.0 163K287E (AAA→GAA) amiAN‑acetylmuramoyl‑l‑alanine amidase I
Reads supporting (aligned to +/- strand):  ref base A (19/24);  new base G (59/61);  total (78/85)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 176 total aligned reads displayed.

TTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGG  >  NC_000913/2553114‑2553307
                                                                                                |                                                                                                 
ttAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTAGTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaa                                                                                                <  1:2366577/100‑1 (MQ=255)
 tAGGCATGGCGGCGTTTGGTCAGAATATCGCTATAGCGTTTGCTGAAGGCGTGGTCAGTTATTTCAACTGGTTCGACAACCAGAAAGCACATGCGgaaaa                                                                                               <  1:2044232/100‑1 (MQ=255)
 tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGgaaaa                                                                                               <  2:4300667/100‑1 (MQ=255)
 tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaaaaa                                                                                               >  1:1427955/1‑100 (MQ=255)
 tAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGAGCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGaaaaa                                                                                               >  1:2143269/1‑100 (MQ=255)
  aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAg                                                                                              >  1:3398392/1‑100 (MQ=255)
   ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                             >  1:4359437/1‑100 (MQ=255)
   ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGc                                                                                             <  1:1856583/100‑1 (MQ=255)
   ggCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGc                                                                                             >  2:511843/1‑100 (MQ=255)
    gCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACACTCGGAAAAGCg                                                                                            <  1:1061119/100‑1 (MQ=255)
       cggcggCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATa                                                                                         <  2:214443/100‑1 (MQ=255)
         gcggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATaa                                                                                        <  1:143684/99‑1 (MQ=255)
         gcggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAg                                                                                       >  2:174152/1‑100 (MQ=255)
            gcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTa                                                                                    >  1:2496164/1‑100 (MQ=255)
             cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                   <  1:3445859/100‑1 (MQ=255)
             cgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTAt                                                                                   >  2:646816/1‑100 (MQ=255)
              gTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATg                                                                                  <  2:2650514/100‑1 (MQ=255)
               tttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGa                                                                                 >  2:2344753/1‑100 (MQ=255)
               tttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGa                                                                                 >  2:2898440/1‑100 (MQ=255)
                ttCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGaa                                                                                >  1:3174005/1‑100 (MQ=255)
                  cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                              <  2:4168138/100‑1 (MQ=255)
                  cGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAc                                                                              >  1:3181570/1‑100 (MQ=255)
                   gTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAcc                                                                             <  1:163193/100‑1 (MQ=255)
                    tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAAccc                                                                            <  2:820401/100‑1 (MQ=255)
                    tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                            <  2:2209396/100‑1 (MQ=255)
                    tCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAAccc                                                                            <  2:2209042/100‑1 (MQ=255)
                       gAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGAc                                                                         >  1:4260276/1‑100 (MQ=255)
                          aaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGca                                                                      <  1:1551272/100‑1 (MQ=255)
                            tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                    <  2:865389/100‑1 (MQ=255)
                            tCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGcaca                                                                    <  2:3371078/100‑1 (MQ=255)
                              gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                  >  1:2802774/1‑100 (MQ=255)
                              gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                  >  2:3269043/1‑100 (MQ=255)
                              gCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTAGGAAAAGCGATAAGTTATGAAACCCGACGCACAcc                                                                  >  1:419808/1‑100 (MQ=255)
                                 acaGTGATTGCTGACGGAGTGATCAGTTATTTCAGCTGGTTCGACAACCAGAAAGGACATTCGGAAAAGCGATAAGTTATGAAACCCGAc                                                                         <  1:576747/90‑1 (MQ=255)
                                 acaGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAgg                                                               >  1:3716601/1‑100 (MQ=255)
                                  caGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGt                                                              >  1:3810122/1‑100 (MQ=255)
                                    gCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTa                                                            <  1:281284/100‑1 (MQ=255)
                                     cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                           <  1:1909874/100‑1 (MQ=255)
                                     cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                           <  2:1198515/100‑1 (MQ=255)
                                     cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAACCACAATCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTaa                                                           <  1:497195/100‑1 (MQ=255)
                                         tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAg                                                       <  2:3644952/100‑1 (MQ=255)
                                          gCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGt                                                      >  2:1927556/1‑100 (MQ=255)
                                           cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                     >  2:2701064/1‑100 (MQ=255)
                                           cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                     <  2:4060426/100‑1 (MQ=255)
                                           cTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGtt                                                     >  2:1037671/1‑100 (MQ=255)
                                            tGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGttt                                                    <  2:2435653/100‑1 (MQ=255)
                                                 gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                               <  1:827789/100‑1 (MQ=255)
                                                 gCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCt                                               <  2:2249240/100‑1 (MQ=255)
                                                  cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGc                                                <  1:4064989/98‑1 (MQ=255)
                                                  cGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTc                                              <  2:1517905/100‑1 (MQ=255)
                                                   gTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTcn                                             <  2:41160/100‑2 (MQ=255)
                                                      aTCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAcc                                          <  2:2496177/100‑1 (MQ=255)
                                                      aTCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAAcc                                          <  2:3181588/100‑1 (MQ=255)
                                                         aGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTc                                       >  1:2964235/1‑100 (MQ=255)
                                                          gTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCa                                      <  1:259296/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  1:1021620/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  2:3959589/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCAATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  1:1925897/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  1:511842/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  1:3015764/100‑1 (MQ=255)
                                                            tATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAgg                                    <  1:847253/100‑1 (MQ=255)
                                                             aTTTCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGa                                   <  2:3450818/100‑1 (MQ=255)
                                                               ttCCACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                 <  1:454768/100‑1 (MQ=255)
                                                               ttCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATa                                 >  2:1370346/1‑100 (MQ=255)
                                                                 ccACTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACg                               <  2:355150/100‑1 (MQ=255)
                                                                    cTGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAtt                            >  2:2817654/1‑100 (MQ=255)
                                                                     tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                           >  1:556368/1‑100 (MQ=255)
                                                                     tGGTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGAttt                           >  2:3767227/1‑100 (MQ=255)
                                                                      ggTTCTACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                          <  1:4058796/100‑1 (MQ=255)
                                                                      ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                          <  1:1986597/100‑1 (MQ=255)
                                                                      ggTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTg                          >  2:2953053/1‑100 (MQ=255)
                                                                       gTTCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                         <  1:1529865/100‑1 (MQ=255)
                                                                        ttCGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                        <  2:1241058/100‑1 (MQ=255)
                                                                        ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGt                         <  1:1724873/99‑1 (MQ=255)
                                                                        ttCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTc                        <  1:2621884/100‑1 (MQ=255)
                                                                          cGACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTGTCTGCTCAACCTTCAGGATACgtgtgg                          >  2:648914/1‑91 (MQ=255)
                                                                           gACAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagc                     >  1:3394993/1‑100 (MQ=255)
                                                                           gACAACCAGAAAGCACATTCGAAAGAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCGTCAGGATACGATTTGTcagc                     <  1:4318774/100‑1 (MQ=255)
                                                                             cAACCAGAAAGTACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGACTTGTcagcag                   >  1:2886846/1‑100 (MQ=255)
                                                                             cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                   >  2:4268805/1‑100 (MQ=255)
                                                                             cAACCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                   >  2:367977/1‑100 (MQ=255)
                                                                             cAAACAGAAAGCACATTAGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcag                   <  1:736754/100‑1 (MQ=255)
                                                                              aaCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCagcagc                  <  1:1359390/100‑1 (MQ=255)
                                                                               aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                 >  2:495653/1‑100 (MQ=255)
                                                                               aCCAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCt                 >  1:120859/1‑100 (MQ=255)
                                                                                 cAGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTg                <  1:3707984/99‑1 (MQ=255)
                                                                                  aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc              <  2:31067/100‑1 (MQ=255)
                                                                                  aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc              >  2:4032300/1‑100 (MQ=255)
                                                                                  aGAAAGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAc              <  1:1955482/100‑1 (MQ=255)
                                                                                     aaGCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTcagcaagatcggaa           >  1:4208988/1‑91 (MQ=255)
                                                                                       gCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCCACCTTAAGGATACGATTTGTCAGCAGCTGAccgc           >  1:4366375/1‑98 (MQ=255)
                                                                                       gCACATTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGAccgccg         >  2:1673201/1‑100 (MQ=255)
                                                                                          cactCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg      <  1:3295524/97‑1 (MQ=255)
                                                                                          caTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCg      >  1:3772825/1‑100 (MQ=255)
                                                                                           aTTCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGa     >  2:180726/1‑100 (MQ=255)
                                                                                           aTTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCACCGTCGa     >  2:68516/1‑100 (MQ=255)
                                                                                            ttCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt    <  1:853939/100‑1 (MQ=255)
                                                                                            ttCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGAt    >  1:2741618/1‑100 (MQ=255)
                                                                                             tCGGAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTTACCGCCGTCGATg   >  1:1635324/1‑100 (MQ=255)
                                                                                             tCGAAAAAGCGATAAGTTATGAAACCCGACGCATACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATg   >  1:2917292/1‑100 (MQ=255)
                                                                                                |                                                                                                 
TTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGTTATGAAACCCGACGCACACCAGGTTAAACAGTTTCTGCTCAACCTTCAGGATACGATTTGTCAGCAGCTGACCGCCGTCGATGG  >  NC_000913/2553114‑2553307

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