Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,349,163 A→G 6.7% V216A (GTG→GCG) ‡ mtgA ← biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,349,1630AG 6.7% 293.6 / 8.5 90V216A (GTG→GCG) ‡mtgAbiosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (41/43);  new base G (2/4);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

ATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAA  >  NC_000913/3349016‑3349310
                                                                                                                                                   |                                                                                                                                                   
aTTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCa                                                                                                                                                     <  2:811366/148‑1 (MQ=255)
    cccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTa                                                                                                                                                 >  1:17984/1‑148 (MQ=255)
     ccTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAg                                                                                                                                                >  2:742030/1‑148 (MQ=255)
      cTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGcc                                                                                                                                              >  4:574175/1‑149 (MQ=255)
       tGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGccc                                                                                                                                             <  3:776132/149‑1 (MQ=255)
          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtca                                                                                                                                          <  3:521520/149‑5 (MQ=255)
          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtc                                                                                                                                           <  1:494534/148‑4 (MQ=255)
          aGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                                                          <  3:122323/149‑1 (MQ=255)
           gAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtcag                                                                                                                                         <  3:423888/149‑6 (MQ=255)
           gAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                                                          >  2:60172/1‑148 (MQ=255)
            aaaGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATgg                                                                                                                                        >  3:350588/1‑149 (MQ=255)
            aaaGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATgg                                                                                                                                        >  3:843503/1‑149 (MQ=255)
              aGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCg                                                                                                                                      <  4:259856/149‑1 (MQ=255)
                aTCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgag                                                                                                                                    <  4:350588/149‑1 (MQ=255)
                aTCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCga                                                                                                                                     <  2:542324/148‑1 (MQ=255)
                 tCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgag                                                                                                                                    >  2:713135/1‑148 (MQ=255)
                  cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                                                                  <  3:327234/149‑1 (MQ=255)
                  cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                                                                  <  3:728531/149‑1 (MQ=255)
                  cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagga                                                                                                                                  <  4:525311/149‑1 (MQ=255)
                  cTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgagg                                                                                                                                   <  1:60172/148‑1 (MQ=255)
                    gCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCgaggag                                                                                                                                 <  1:521274/148‑1 (MQ=255)
                      aaCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGAc                                                                                                                               <  2:141651/148‑1 (MQ=255)
                       aCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCg                                                                                                                                      >  1:724208/1‑140 (MQ=255)
                       aCCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCg                                                                                                                                      <  2:724208/140‑1 (MQ=255)
                        ccTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                                                                                                            <  4:860153/149‑1 (MQ=255)
                         cTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACttt                                                                                                                            <  1:742030/148‑1 (MQ=255)
                          tCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGa                                                                                                                          >  4:476175/1‑149 (MQ=255)
                             tGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                                                                                                             >  2:84769/1‑143 (MQ=255)
                             tGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                                                                                                             <  1:84769/143‑1 (MQ=255)
                                     aaaCGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGa                                                                                                               >  3:209246/1‑149 (MQ=255)
                                           tGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtca                                                                                                                                          <  4:637135/116‑5 (MQ=255)
                                           tGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGTGCGTAGccgtca                                                                                                                                          >  3:637135/1‑112 (MQ=255)
                                           tGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAgg                                                                                                          <  2:401775/148‑1 (MQ=255)
                                             gCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTaa                                                                                                       <  3:390550/149‑1 (MQ=255)
                                               ccTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATa                                                                                                     <  3:635066/149‑1 (MQ=255)
                                                cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                                                                                                             <  3:372799/124‑1 (MQ=255)
                                                cTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACtt                                                                                                                             >  4:372799/1‑124 (MQ=255)
                                                    cTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCga                                                                                                                                     >  2:361914/1‑112 (MQ=255)
                                                    cTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCga                                                                                                                                     <  1:361914/112‑1 (MQ=255)
                                                    cTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCt                                                                                                <  3:541002/149‑1 (MQ=255)
                                                         cgcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTCACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCa                                                                                            >  2:155773/1‑148 (MQ=255)
                                                          gcgAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGt                                                                                          >  4:373743/1‑149 (MQ=255)
                                                               gCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGc                                                                                                                                       <  1:14398/99‑1 (MQ=255)
                                                               gCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGc                                                                                                                                       >  2:14398/1‑99 (MQ=255)
                                                                  tAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGc                                                                                  <  3:574175/149‑1 (MQ=255)
                                                                     tCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTc                                                                               >  4:97472/1‑149 (MQ=255)
                                                                           tGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGAcc                                                                          <  2:513869/148‑1 (MQ=255)
                                                                              tgttgtTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGt                                                                      <  3:373743/149‑1 (MQ=255)
                                                                                 tgttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTaa                                                                    >  2:809271/1‑148 (MQ=255)
                                                                                  gttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCgg                                                                        <  4:59349/143‑1 (MQ=255)
                                                                                  gttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCgg                                                                        >  3:59349/1‑143 (MQ=255)
                                                                                   ttgCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGtt                                                                 >  4:546709/1‑149 (MQ=255)
                                                                                      cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTg                                                               >  2:499928/1‑148 (MQ=255)
                                                                                      cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGc                                                              >  3:502625/1‑149 (MQ=255)
                                                                                      cATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGc                                                              >  3:382290/1‑149 (MQ=255)
                                                                                       aTAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGc                                                              <  1:276485/148‑1 (MQ=255)
                                                                                       aTAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                             >  4:222557/1‑149 (MQ=255)
                                                                                        tAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                             >  2:448655/1‑148 (MQ=255)
                                                                                        tAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCt                                                             >  1:674018/1‑148 (MQ=255)
                                                                                        tAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTc                                                            >  4:689561/1‑149 (MQ=255)
                                                                                          aaCGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                                                          >  3:120773/1‑69 (MQ=255)
                                                                                          aaCGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGAt                                                                                                                                          <  4:120773/69‑1 (MQ=255)
                                                                                                tCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGtt                                                     <  1:809271/148‑1 (MQ=255)
                                                                                                  aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATa                                                                                                     >  3:27699/1‑98 (MQ=255)
                                                                                                  aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATa                                                                                                     <  4:27699/98‑1 (MQ=255)
                                                                                                  aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                  <  3:476175/149‑1 (MQ=255)
                                                                                                  aCCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtg                                                   <  1:129310/148‑1 (MQ=255)
                                                                                                   ccGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                  >  2:47600/1‑148 (MQ=255)
                                                                                                        tAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGaaa                                             >  1:554823/1‑148 (MQ=255)
                                                                                                           cTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATaa                                          >  2:388507/1‑148 (MQ=255)
                                                                                                               tACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                  >  2:559001/1‑136 (MQ=255)
                                                                                                               tACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTtgt                                                  <  1:559001/136‑1 (MQ=255)
                                                                                                                aCATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTtgt                                    <  4:843503/149‑1 (MQ=255)
                                                                                                                   tCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGcc                                 >  4:245909/1‑149 (MQ=255)
                                                                                                                   tCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGcc                                 >  4:94333/1‑149 (MQ=255)
                                                                                                                   tCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGcc                                 <  4:382290/149‑1 (MQ=255)
                                                                                                                   tCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGc                                  >  1:609565/1‑148 (MQ=255)
                                                                                                                     tGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCgg                                                        <  4:395356/124‑1 (MQ=255)
                                                                                                                     tGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCgg                                                        >  3:395356/1‑124 (MQ=255)
                                                                                                                        cGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAg                             <  2:554823/148‑1 (MQ=255)
                                                                                                                          tAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCt                           <  1:499928/148‑1 (MQ=255)
                                                                                                                              aTCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGAc                      <  3:190919/149‑1 (MQ=255)
                                                                                                                                ccACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGttt                                                                <  1:139559/105‑1 (MQ=255)
                                                                                                                                ccACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGttt                                                                >  2:139559/1‑105 (MQ=255)
                                                                                                                                ccACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACg                     >  1:133345/1‑148 (MQ=255)
                                                                                                                                  aCGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGAATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGcc                   >  1:597416/1‑148 (MQ=255)
                                                                                                                                     ccTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                                                    >  1:531665/1‑112 (MQ=255)
                                                                                                                                     ccTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGttt                                                    <  2:531665/112‑1 (MQ=255)
                                                                                                                                           ggCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCg         >  3:451544/1‑149 (MQ=255)
                                                                                                                                              tACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCg      <  4:566124/149‑1 (MQ=255)
                                                                                                                                                  cACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCa   <  1:155773/148‑1 (MQ=255)
                                                                                                                                                   aCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCaa  <  1:388507/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTGCCCTGCAGAAAGATCTGCAACCTCGTGAGAGCAAAACGCTGGCCCTCACTTCGCGCGAAGCTTAATCCAGCTGGTGTTGTTGCATAAACGGCTCACCGCCTAACTGATACATCTGCCGTAAAATCCACGCCTGACGGCTACGCACGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTCGCCAA  >  NC_000913/3349016‑3349310

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