Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 981,315 A→G Y90C (TAT→TGT)  ldtD → murein L,D‑transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913981,3150AG95.8% 314.3 / ‑5.3 96Y90C (TAT→TGT) ldtDmurein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (41/50);  minor base T (4/0);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 4.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTATCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTT  >  NC_000913/981169‑981459
                                                                                                                                                    |                                                                                                                                                 
ggTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCCGCCACTGCCTGAAAGGGCGGGAGAAAAAACCCGCACCCCAAACGAATTTCCATTACCCGCAGGGTACAAACCGGGTTGTCtt                                                                                                                                              >  2:48684/1‑151 (MQ=255)
 gTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGGCCCGCCAATGCCTGAAGGGGGGGCAGAAAAAACCCCCACGCAAAACGAATCTCAATTAACCGCAGGGTACAAGCCGGGTTGTCtt                                                                                                                                              >  2:373842/1‑150 (MQ=255)
  tGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCCATTCCTGAAGGGGCGGCAGAAAAAAACCGCACCCCAAACGAATCTCAATTAACCGCCGGGTACAAACCGGTTTGTCttt                                                                                                                                             >  4:596989/1‑149 (MQ=255)
   ggCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTaa                                                                                                                                            <  2:850708/150‑1 (MQ=255)
     cTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCCGCCAGTGCCTGAAGGGGCGGGAGAAAAAACCCGCACGCAAATCGCATCTCAAATAACCCCAGGGTTAAAACCCGGTTTTCTTAAcc                                                                                                                                          >  4:1060612/1‑150 (MQ=255)
      tGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTGCAGCCATTGCCTGAAGGGGCGGCAGAAAAAACCCGCACGCCAAATGAATCTCAATTACCCGCCGGGTAACAACCGGGTTGTCTTAACCa                                                                                                                                         >  4:1167571/1‑150 (MQ=255)
      tGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCAGCCAATGCCTCAAAGGGCGGCAGAAAAAACCCGCCCGCAAAACGAATATCCATTACCCCCCGGGTACAAAACGGTTTTTCTTaaac                                                                                                                                          >  1:387909/1‑147 (MQ=255)
       gTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCa                                                                                                                                         <  1:370915/149‑1 (MQ=255)
         cAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGGCCCACCAATGCCTGAAGGAGCGGGCGGAAAAAACCCCCCCCAAAACGAAACCCAATTAACCCCAGGGTGCAAAACGGGTTGTCtttaacaacct                                                                                                                                     >  1:142065/1‑142 (MQ=255)
         cAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCt                                                                                                                                      <  3:230728/150‑1 (MQ=255)
             gAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAGGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAAc                                                                                                                                 <  4:1109027/151‑1 (MQ=255)
              aaCAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGACCCGCCACTTCCTGAAGGGGCGGGAGAAAAAACCCGCCCGCAAATAGAATTTCAATTACCCCCAGGGTACAAGCCGGGTTTTTTTAACCAGCTTCCAc                                                                                                                                 >  1:394450/1‑150 (MQ=255)
               aCAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCCTTTCCTTAAGGTGCGGCAGAAAAAGCCCGCACGCCAATCGAATCTCAATTACCCCCCGGGTACAAACCCGGTTGTCTTTACCAGCTTCaa                                                                                                                                  >  3:133934/1‑148 (MQ=255)
               aCAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACtg                                                                                                                               <  3:350807/151‑1 (MQ=255)
                taGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCAGCCATTGCCTGAAGGGGCGGCAGAAAAAACCCGCACGCAAAACGAATCTCAATTACCCGCCGGGTACAAACCGGGTTTTCTTAAACAGATTTAACtgt                                                                                                                              >  2:718288/2‑151 (MQ=255)
                cAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACtgt                                                                                                                              <  2:653900/151‑1 (MQ=255)
                 aGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCCACCCTTTCCTTAAAGAGGGGGCGAAAAAGACCCCCCGCAAAACGAAACTCAATTATCCGCAGGGTACAAACCGGGTTTTCTTAACCAGCTTCAACtgtt                                                                                                                             >  4:703974/1‑151 (MQ=255)
                  gggCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACtgttg                                                                                                                            <  1:430505/151‑1 (MQ=255)
                  gggCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACtgttg                                                                                                                            <  4:1439959/151‑1 (MQ=255)
                        ggCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATg                                                                                                                        <  3:202363/149‑1 (MQ=255)
                         gCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAAGGGCGGCAGAAAAAACCCGCACGCCAATCGAATCTCAATTACCCGCAGGGTACAAGCCGGGTTGTCTTAACCAGCTTCAACTGTTGTATGcc                                                                                                                      >  4:183015/1‑150 (MQ=255)
                           aCTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCAACCACTGCCTGAAGGTGGGGCCGAAAAAACCCCCACGCCAAACGAATCTCAATTACCCCCCGGTTACAAACCCGGTTGTCTTAACCAACTTCAAATGTTGTTTGCCGc                                                                                                                    >  2:301145/1‑150 (MQ=255)
                            cTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCa                                                                                                                   <  4:235687/150‑1 (MQ=255)
                             tGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAACCATTGCCTGAAAGGGCGGCAGAAAAAACCCGCACGCAAATAGAAACTCAATTAACCGCAGGGTACAAACCCGGTTGTCTTTAACAACCTTCACTGTTGTATGcccc                                                                                                                    >  4:1111239/1‑146 (MQ=255)
                              gCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGGTCCAGCCCCTGCCTGAAAGGGGGGCGGAAAAAACCCCCACGCAAAACGAATCTCAATTACCCCCCGGGTACAAACCGGGTTGTTTTTAACAACTTTAACTGTTTTTTGCCggcc                                                                                                                  >  2:1305610/1‑146 (MQ=255)
                              gCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCAcgc                                                                                                                <  4:34419/151‑1 (MQ=255)
                              gCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCAcgc                                                                                                                <  2:616914/151‑1 (MQ=255)
                                cGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGa                                                                                                              <  4:731243/151‑1 (MQ=255)
                                    ggCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATg                                                                                                          <  2:183196/151‑1 (MQ=255)
                                        cAAGCCACGGCAATAATGGCGGGGGTCCAACCATTTCCTGAAGGGGCGGCCGAAAAAACCCGCCCGCAAAACGAATCACAATTACCCGCAGGGTAAAAGCCGGGTTGTCTTAACCAACTTTAACTGTTGTATGCCGCACGCGGTATGCaa                                                                                                       >  4:59274/1‑150 (MQ=255)
                                         aaGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAAcc                                                                                                     <  1:98214/151‑1 (MQ=255)
                                         aaGCCACGGCAATAATGGCGGGGAGCCAGCCACTGCCCTAATAAGGGGGCGAAAAAGACCGCACGCAAAACGAATCTCCATTACCCGCAGGGTAAAAAACGGGTTGTGTTAACAAGCTTCAAATTTTGTATGCCGCACGCGATATCAAAcc                                                                                                     >  2:1641842/1‑151 (MQ=255)
                                            ccACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCAt                                                                                                  <  1:428272/151‑1 (MQ=255)
                                              aCGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCAtgt                                                                                                <  2:1279896/151‑1 (MQ=255)
                                               cGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCAtgt                                                                                                <  2:704543/150‑1 (MQ=255)
                                                 gCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTggg                                                                                             <  2:1001745/151‑1 (MQ=255)
                                                  cAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGa                                                                                            <  1:366770/151‑1 (MQ=255)
                                                   aataatGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGaa                                                                                           <  2:995132/151‑1 (MQ=255)
                                                   aataatGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGaa                                                                                           <  2:775427/151‑1 (MQ=255)
                                                     taatGGCGGGGATCCAGCCACTGCCTTAAGGGGGGGCCGAAAAAACCCCCACGCCAAACGAAACTCAATTAACCGCCGGGTTCAAACCGGGTTGTCTTAACCAACCTCAACTGTTGTATGCCGCACGCCGTATGCAACCCCAGTGGGgaa                                                                                          >  1:22559/1‑150 (MQ=255)
                                                      aatGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGaaaa                                                                                         <  4:1303687/150‑1 (MQ=255)
                                                         ggCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCCCACGCAATTCGAATCCCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAATTGTTGTAGGCCGCACGCGATATGCAACCCATGGGGGAAAACCg                                                                                      >  4:1128525/1‑150 (MQ=255)
                                                          gCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTg                                                                                    >  1:127365/1‑151 (MQ=255)
                                                           cGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGa                                                                                   <  2:433113/151‑1 (MQ=255)
                                                             gggATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATg                                                                                 >  4:1233155/1‑151 (MQ=255)
                                                               gATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCt                                                                               <  3:103850/151‑1 (MQ=255)
                                                                  ccAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGtt                                                                            >  4:995820/1‑151 (MQ=255)
                                                                   cAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGtt                                                                            >  2:819785/1‑150 (MQ=255)
                                                                    aGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTaa                                                                          <  3:409892/151‑1 (MQ=255)
                                                                     gCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTaaa                                                                         <  4:821378/151‑1 (MQ=255)
                                                                               gAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAg                                                                <  2:1117263/150‑1 (MQ=255)
                                                                               gAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGc                                                               >  1:107866/1‑151 (MQ=255)
                                                                                 aGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCa                                                              >  1:30470/1‑150 (MQ=255)
                                                                                      cGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCt                                                        >  2:1431905/1‑151 (MQ=255)
                                                                                               aaaGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAgg                                                >  2:1490936/1‑151 (MQ=255)
                                                                                                aaGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGt                                               >  2:483628/1‑151 (MQ=255)
                                                                                                 aGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTg                                              >  4:160038/1‑151 (MQ=255)
                                                                                                   cccGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGc                                            <  2:1569801/151‑1 (MQ=255)
                                                                                                    ccGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCg                                           >  4:1325243/1‑151 (MQ=255)
                                                                                                     cgcacgcaAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCg                                           >  2:1482945/1‑150 (MQ=255)
                                                                                                     cgcacgcaAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCg                                           >  2:1156537/1‑150 (MQ=255)
                                                                                                       cacgcaAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGAtt                                        <  2:1260588/151‑1 (MQ=255)
                                                                                                           caAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGcca                                    <  2:106238/151‑2 (MQ=255)
                                                                                                            aaaTCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCgg                                   >  1:148708/1‑151 (MQ=255)
                                                                                                             aaTCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGt                                  <  3:287170/151‑1 (MQ=255)
                                                                                                                cGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTc                               <  3:106095/151‑1 (MQ=255)
                                                                                                                 gAATCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTcc                              >  4:1163676/1‑151 (MQ=255)
                                                                                                                  aaTCTCAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCa                             <  4:1590579/151‑1 (MQ=255)
                                                                                                                      tcAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCg                         >  2:80852/1‑151 (MQ=255)
                                                                                                                       cAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAg                            <  2:857259/147‑1 (MQ=255)
                                                                                                                       cAATTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGcc                          <  3:85917/149‑1 (MQ=255)
                                                                                                                        aaTTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGc                        >  2:144960/1‑150 (MQ=255)
                                                                                                                        aaTTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCa                       <  2:1323590/151‑1 (MQ=255)
                                                                                                                        aaTTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCa                       <  3:265713/151‑1 (MQ=255)
                                                                                                                         aTTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGc                        >  4:559415/1‑149 (MQ=255)
                                                                                                                         aTTACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCAGGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGc                           <  4:1653428/146‑1 (MQ=255)
                                                                                                                          ttACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGt                     >  2:1246708/1‑151 (MQ=255)
                                                                                                                           tACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGtt                    >  4:408332/1‑151 (MQ=255)
                                                                                                                           tACCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGt                     <  4:1535714/150‑1 (MQ=255)
                                                                                                                            aCCCGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGttt                   >  3:295375/1‑151 (MQ=255)
                                                                                                                             cccGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTta                  <  2:596167/151‑1 (MQ=255)
                                                                                                                             cccGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTta                  <  4:629593/151‑1 (MQ=255)
                                                                                                                             cccGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTta                  <  4:1168570/151‑1 (MQ=255)
                                                                                                                              ccGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTtaa                 <  1:358664/151‑1 (MQ=255)
                                                                                                                               cGCAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTtaat                >  1:337068/1‑151 (MQ=255)
                                                                                                                                 cAGGTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTtaataa              <  4:306244/151‑1 (MQ=255)
                                                                                                                                    ggTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAAt            >  3:324921/1‑151 (MQ=255)
                                                                                                                                    ggTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAAt            <  2:1450756/151‑1 (MQ=255)
                                                                                                                                     gTTACAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATg           <  4:48779/151‑1 (MQ=255)
                                                                                                                                        aCAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATggg         <  2:1621756/150‑1 (MQ=255)
                                                                                                                                        aCAAGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGt        >  2:1184724/1‑151 (MQ=255)
                                                                                                                                           aGCCGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTaga     <  1:353823/151‑1 (MQ=255)
                                                                                                                                             ccGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGt   >  1:89369/1‑151 (MQ=255)
                                                                                                                                             ccGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGt   <  2:1109833/151‑1 (MQ=255)
                                                                                                                                              cGGTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGtt  >  2:1450343/1‑151 (MQ=255)
                                                                                                                                               ggTTTGTCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGtt  >  2:964526/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTGGCTGTCAGTGAACAGGGCGAGGCACTGCCGCAGGCGCAAGCCACGGCAATAATGGCGGGGATCCAGCCATTGCCTGAAGGTGCGGCAGAAAAAGCCCGCACGCAAATCGAATCTCAATTACCCGCAGGTTACAAGCCGGTTTATCTTAACCAGCTTCAACTGTTGTATGCCGCACGCGATATGCAACCCATGTGGGAAAACCGTGATGCTGTTAAAGCCTTCCAGCAACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTT  >  NC_000913/981169‑981459

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