Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,946,757 T→C Y141C (TAC→TGC)  mltA ← membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,946,7570TC97.5% 296.2 / ‑2.4 81Y141C (TAC→TGC) mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (44/35);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

CACTACCCTGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGGCGCATA  >  NC_000913/2946626‑2946891
                                                                                                                                   |                                                                                                                                      
cACTACCCTGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAgg                                                                                                                                 >  2:467459/1‑139 (MQ=255)
    aCCCTGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATAc                                                                                                                             <  1:872803/139‑1 (MQ=255)
        tGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGa                                                                                                                         >  2:594595/1‑139 (MQ=255)
        tGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGa                                                                                                                         >  2:295808/1‑139 (MQ=255)
        tGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGa                                                                                                                         >  1:798538/1‑139 (MQ=255)
         gCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGaa                                                                                                                        >  2:924555/1‑139 (MQ=255)
             aTCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCg                                                                                                                    >  2:550018/1‑139 (MQ=255)
               ccATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGcc                                                                                                                  <  2:741487/139‑1 (MQ=255)
               ccATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGcc                                                                                                                  >  1:486088/1‑139 (MQ=255)
               ccATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGcc                                                                                                                  >  1:425672/1‑139 (MQ=255)
                    aTGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGc                                                                                                             >  2:437011/1‑139 (MQ=255)
                    aTGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGc                                                                                                             <  2:425672/139‑1 (MQ=255)
                         gTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTa                                                                                                        >  2:1171182/1‑139 (MQ=255)
                                aTCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGgcgcg                                                                                                 <  1:924555/139‑1 (MQ=255)
                                   aGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCt                                                                                              >  2:1000636/1‑139 (MQ=255)
                                          tACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGcc                                                                                                                  <  1:685757/112‑1 (MQ=255)
                                          tACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGcc                                                                                                                  >  2:685757/1‑112 (MQ=255)
                                          tACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTAc                                                                                       <  2:798538/139‑1 (MQ=255)
                                              gTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGc                                                                                   <  1:550018/139‑1 (MQ=255)
                                                 aGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTAc                                                                                       <  2:1034500/132‑1 (MQ=255)
                                                 aGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTAc                                                                                       >  1:1034500/1‑132 (MQ=255)
                                                 aGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCgtg                                                                                >  1:758204/1‑139 (MQ=255)
                                                   cGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGta                                                                              >  2:1035971/1‑139 (MQ=255)
                                                     agaATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaat                                                                            <  1:1035971/139‑1 (MQ=255)
                                                     agaATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaat                                                                            <  2:588873/139‑1 (MQ=255)
                                                                 tCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTa                                                                                        <  1:1106145/115‑1 (MQ=255)
                                                                 tCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTa                                                                                        >  2:1106145/1‑115 (MQ=255)
                                                                  cACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTg                                                               >  1:276156/1‑139 (MQ=255)
                                                                   aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCACATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGc                                                                                   >  1:81295/1‑118 (MQ=255)
                                                                   aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGc                                                                                   <  2:81295/118‑1 (MQ=255)
                                                                    cTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCa                                                             <  1:467459/139‑1 (MQ=255)
                                                                    cTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAAGAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCa                                                             <  1:437011/139‑1 (MQ=255)
                                                                     tCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGCCCACGTTTTGGCGGCATACGGCAAATCGGATACTGGAACTCGCACTGGCGGGTATGGCGCGCCTGAATTATCGGTGTGTAATAACCGGTAAACTGCAc                                                            >  1:1080946/1‑139 (MQ=255)
                                                                     tCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCAc                                                            >  1:699118/1‑139 (MQ=255)
                                                                     tCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCAc                                                            >  1:63090/1‑139 (MQ=255)
                                                                       aaTGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGt                                                          >  1:120538/1‑139 (MQ=255)
                                                                        aTGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt                                                         >  2:1177200/1‑139 (MQ=255)
                                                                        aTGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGtt                                                         >  2:637892/1‑139 (MQ=255)
                                                                         tGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGgcgcg                                                                                                 <  1:437146/98‑1 (MQ=255)
                                                                         tGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGgcgcg                                                                                                 >  2:437146/1‑98 (MQ=255)
                                                                         tGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTa                                                        >  1:173907/1‑139 (MQ=255)
                                                                         tGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTa                                                        <  2:63090/139‑1 (MQ=255)
                                                                              ccGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATa                                                   <  2:486088/139‑1 (MQ=255)
                                                                              ccGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATa                                                   <  1:500850/139‑1 (MQ=255)
                                                                               cgcgTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAgg                                                                                                                                 <  2:428998/60‑1 (MQ=255)
                                                                               cgcgTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAgg                                                                                                                                 >  1:428998/1‑60 (MQ=255)
                                                                               cgcgTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAg                                                  >  1:1111801/1‑139 (MQ=255)
                                                                                 cgTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGtt                                                >  1:969612/1‑139 (MQ=255)
                                                                                  gTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTg                                               <  1:1171182/139‑1 (MQ=255)
                                                                                   tAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGt                                              >  1:517626/1‑139 (MQ=255)
                                                                                      aTCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCgg                                           <  2:120538/139‑1 (MQ=255)
                                                                                          ccGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCAcc                                       >  1:339361/1‑139 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                                                                     <  1:1096309/76‑1 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATgg                                                                                                     >  2:1096309/1‑76 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTa                                                                                        >  1:766278/1‑89 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTa                                                                                        <  2:766278/89‑1 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCgg                                                                                    <  2:55637/93‑1 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCgg                                                                                    >  1:55637/1‑93 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaata                                                                           >  2:632612/1‑102 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaata                                                                           <  1:632612/102‑1 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCa                                                             >  2:721423/1‑116 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCa                                                             <  1:721423/116‑1 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCt                                      >  1:180876/1‑139 (MQ=255)
                                                                                           cGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCt                                      >  1:902001/1‑139 (MQ=255)
                                                                                             cACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaat                                                                            <  1:1060785/99‑1 (MQ=255)
                                                                                             cACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGtaat                                                                            >  2:1060785/1‑99 (MQ=255)
                                                                                                  gaCGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCgg                                           >  2:268329/1‑127 (MQ=255)
                                                                                                  gaCGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCgg                                           <  1:268329/127‑1 (MQ=255)
                                                                                                   aCGGCAGACGACCACGTTTTGGCGGCATACGTCAAATAGGAGACTGGAACGCGCGCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCAGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTg                              <  2:1080946/139‑1 (MQ=255)
                                                                                                   aCGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACt                                                                                                                            >  2:336663/1‑45 (MQ=255)
                                                                                                   aCGGCAGACGACCACGTTTTGGCGGCATACGGCAAATAGGATACt                                                                                                                            <  1:336663/45‑1 (MQ=255)
                                                                                                          acgacCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCa                       >  1:1017956/1‑139 (MQ=255)
                                                                                                           cgacCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCAt                      <  2:517626/139‑1 (MQ=255)
                                                                                                             acCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCt                                      >  2:950366/1‑121 (MQ=255)
                                                                                                             acCACGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCt                                      <  1:950366/121‑1 (MQ=255)
                                                                                                                aCGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATg                 <  2:249367/139‑1 (MQ=255)
                                                                                                                 cGTTTTGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGc                >  1:428055/1‑139 (MQ=255)
                                                                                                                    tttGGCGGCATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGa             <  2:969612/139‑1 (MQ=255)
                                                                                                                            cATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGaa            <  2:285629/132‑1 (MQ=255)
                                                                                                                            cATACGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGaa            >  1:285629/1‑132 (MQ=255)
                                                                                                                               aCGGCAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGGCGCata  <  2:428055/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
CACTACCCTGCACATCCATAATGAAGTTATCCATCAGGGAGTTACTGTAAGCGAGAATATATTTATCACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGCGGGTATGGCGCGCCTGAATTACCGGCGTGTAATAACCGGTAAACTGCACGTTACCATAGTTGTCGGCACCTTCCATCTGCCAGGCATCAATGCCGAACTGGCGCATA  >  NC_000913/2946626‑2946891

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: