Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,323,120 A→T 5.8% D44E (GAT→GAA mrcB ← transglycosylase/transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,323,1200AT 5.8% 264.6 / 4.9 86D44E (GAT→GAAmrcBtransglycosylase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base T (3/2);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TTTCGCGACGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGC  >  NC_002947/5322975‑5323268
                                                                                                                                                 |                                                                                                                                                     
tttCGCGACGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATc                                                                                                                                                     >  7:381855/1‑148 (MQ=255)
       aCGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTaaa                                                                                                                                             <  8:363118/149‑1 (MQ=255)
              ccGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGaa                                                                                                                                      <  5:330400/149‑1 (MQ=255)
               ccccattaTCTAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACcgcg                                                                                                                                        <  1:348508/139‑1 (MQ=255)
                    gcATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATc                                                                                                                                                     >  7:103400/1‑128 (MQ=255)
                    gcATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATc                                                                                                                                                     <  8:103400/128‑1 (MQ=255)
                      aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  4:41431/1‑149 (MQ=255)
                      aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  5:162631/1‑149 (MQ=255)
                      aTCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCaa                                                                                                                              >  5:412471/1‑149 (MQ=255)
                         gAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCTTCGCGGTTTACCTCGATGCCgtcgtcca                                                                                                                           <  7:220017/149‑8 (MQ=255)
                           gCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACcgcg                                                                                                                                        >  2:348508/1‑134 (MQ=255)
                                     gAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCaa                                                                                                               >  6:338737/1‑149 (MQ=255)
                                     gAAGTCGTGCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACTGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCaa                                                                                                               <  7:129252/149‑1 (MQ=255)
                                      aaGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              <  4:62988/149‑1 (MQ=255)
                                       aGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              >  6:154030/1‑148 (MQ=255)
                                       aGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAAc                                                                                                              >  6:111740/1‑148 (MQ=255)
                                           gTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGc                                                                                                         >  3:56789/1‑149 (MQ=255)
                                              cTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTg                                                                                                       <  7:283766/148‑1 (MQ=255)
                                                  cTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:258122/1‑149 (MQ=255)
                                                  cTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:170452/1‑149 (MQ=255)
                                                   tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  8:332094/1‑148 (MQ=255)
                                                   tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  1:16129/1‑148 (MQ=255)
                                                   tCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTc                                                                                                 >  4:272176/1‑149 (MQ=255)
                                                   tCAACTTCTGGCCGGTGAACAGCTCCAGCGATCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCtt                                                                                                  >  1:243751/1‑148 (MQ=255)
                                                   tCAAATTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTGCAGCGGTTGCCGGAGAACTTTTCCGGGACGACAGCTTCGCGTTTAACCGCGCACCCGGCGACGATCACCAGGCCAGCCAGGCTgcgcgt                                                                                                  >  6:327568/1‑146 (MQ=255)
                                                     aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCAc                                                                                                                       >  7:89840/1‑125 (MQ=255)
                                                     aaCTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCAc                                                                                                                       <  8:89840/125‑1 (MQ=255)
                                                          cTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcc                                                                                           <  5:317850/148‑1 (MQ=255)
                                                             gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagc                                                                                       >  4:209194/1‑149 (MQ=255)
                                                             gCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagccca                                                                                         <  6:41951/147‑1 (MQ=255)
                                                               cGGTGCACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCACCGATCACCAGGCCAACCAGCCAGAGCTTCAAcgcccagcc                                                                                      >  2:86583/1‑148 (MQ=255)
                                                               cGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAagcccagcc                                                                                      >  3:181976/1‑148 (MQ=255)
                                                                ggTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAgcccagccca                                                                                    <  5:405875/149‑1 (MQ=255)
                                                                 gTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGcccagcccag                                                                                   >  3:147204/1‑149 (MQ=255)
                                                                         cTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGgcgc                                                                          <  4:147204/149‑1 (MQ=255)
                                                                            cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGcgcg                                                                         <  3:199471/147‑1 (MQ=255)
                                                                            cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  4:181976/148‑1 (MQ=255)
                                                                            cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  6:412471/148‑1 (MQ=255)
                                                                            cAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCgg                                                                        <  5:61292/148‑1 (MQ=255)
                                                                                 gTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCgg                                                                  <  2:238902/149‑1 (MQ=255)
                                                                                  tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  7:170452/149‑1 (MQ=255)
                                                                                  tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  6:162631/149‑1 (MQ=255)
                                                                                  tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 >  2:364542/1‑149 (MQ=255)
                                                                                  tCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGc                                                                 <  3:272176/149‑1 (MQ=255)
                                                                                    gggCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCg                                                               >  8:125376/1‑149 (MQ=255)
                                                                                        aTACACTTTGGCCGGGATGGTCCAGCGCTTGTCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCg                                                           >  3:165360/1‑149 (MQ=255)
                                                                                        aTACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCg                                                           >  4:301207/1‑149 (MQ=255)
                                                                                         tACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCgg                                                          <  5:154030/149‑1 (MQ=255)
                                                                                         tACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCgg                                                          <  2:243751/149‑1 (MQ=255)
                                                                                          acacTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCggg                                                         >  8:144211/1‑149 (MQ=255)
                                                                                            acTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCg                                                       >  7:361852/1‑149 (MQ=255)
                                                                                                tGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGttt                                                    <  3:209194/148‑1 (MQ=255)
                                                                                                     gggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAgg                                              >  5:51661/1‑149 (MQ=255)
                                                                                                      ggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAggg                                             <  3:301207/149‑1 (MQ=255)
                                                                                                      ggATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAggg                                             >  8:125232/1‑149 (MQ=255)
                                                                                                        aTGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGt                                           <  3:41431/149‑1 (MQ=255)
                                                                                                        aTGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGAGAGCGGCCGGTCGGGCGTTTCTGAGGGGt                                           <  3:181187/149‑1 (MQ=255)
                                                                                                          ggTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGc                                         >  2:329150/1‑149 (MQ=255)
                                                                                                              cAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGc                                                                                                                  <  3:110480/73‑1 (MQ=255)
                                                                                                                      gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  5:111740/148‑1 (MQ=255)
                                                                                                                      gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  8:381855/148‑1 (MQ=255)
                                                                                                                      gCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcg                               <  7:332094/147‑1 (MQ=255)
                                                                                                                      gCCGGAGAACTGTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgg                              <  8:269054/148‑1 (MQ=255)
                                                                                                                       ccGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcggg                             <  1:161810/148‑1 (MQ=255)
                                                                                                                        cGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCATCAAAGCCCAGCCCAGCAAAGCGGCGGGACCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggt                            >  1:131354/1‑149 (MQ=255)
                                                                                                                           agaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGGGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                        <  5:338737/149‑1 (MQ=255)
                                                                                                                            gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  1:190714/1‑149 (MQ=255)
                                                                                                                            gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  7:9595/1‑149 (MQ=255)
                                                                                                                            gaACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcgg                        <  1:364542/148‑1 (MQ=255)
                                                                                                                              aCTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCTAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTc                     <  7:258122/149‑1 (MQ=255)
                                                                                                                              aCTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTc                     <  7:144211/149‑1 (MQ=255)
                                                                                                                                ttttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCAt                   >  5:186327/1‑149 (MQ=255)
                                                                                                                                ttttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCAt                   >  2:119147/1‑149 (MQ=255)
                                                                                                                                 tttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATg                  >  7:107970/1‑149 (MQ=255)
                                                                                                                                 tttCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATg                  >  4:312380/1‑149 (MQ=255)
                                                                                                                                     cTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTGAGGTTCGGGTCATggcgg              >  1:138990/1‑149 (MQ=255)
                                                                                                                                       ggACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCggg                                      <  6:186998/123‑1 (MQ=255)
                                                                                                                                       ggACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCggg                                      >  5:186998/1‑123 (MQ=255)
                                                                                                                                        gacgacGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGgcggcg            <  2:354000/148‑1 (MQ=255)
                                                                                                                                        gacgacGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGcggcgg           <  3:312380/149‑1 (MQ=255)
                                                                                                                                         acgacgGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGa          <  7:125232/149‑1 (MQ=255)
                                                                                                                                           gacgGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGAtt        <  1:119147/149‑1 (MQ=255)
                                                                                                                                            ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       <  5:217062/132‑1 (MQ=255)
                                                                                                                                            ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATttcgggttcggg                       >  6:217062/2‑133 (MQ=255)
                                                                                                                                            ccgGCTTCGCGGTTTACCTCGATGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGAGGGGAGTTCGGGTTCGGGTCATGGCGGCGGATta       >  8:351027/2‑149 (MQ=255)
                                                                                                                                                cATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACg   <  4:165360/149‑1 (MQ=255)
                                                                                                                                                 aTCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGc  >  7:68115/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTCGCGACGGTAGCCGAGCGCATCGAGCTCGGTGAGGAAGTCGTTCTTGCTCAACTTCTGGCCGGTGAACAGCTCCAGCGGTCGGGCATACACTTTGGCCGGGATGGTCCAGCGCTTGCCGGAGAACTTTTCCTGGACGACGGCATCGAGGTAAACCGCGAAGCCGGCAACGATCACCAGGCCAACCAGGCTGAGCTTCAAAGCCCAGCCCAGCCAGGCGCGGGAGCGGCCGGTCGGGCGTTTCTGAGGGGTGCGGGGATTTCGGGTTCGGGTCATGGCGGCGGATTATACGC  >  NC_002947/5322975‑5323268

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