Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,819,067 C→A 36.6% intergenic (‑132/‑124) mltB ← / → srlA membrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,819,0670CA36.6% 381.0 / 260.9 257intergenic (‑132/‑124)mltB/srlAmembrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (70/93);  new base A (42/52);  total (112/145)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 266 total aligned reads displayed.

ATAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTG  >  CP009273/2818918‑2819215
                                                                                                                                                     |                                                                                                                                                    
aTAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTc                                                                                                                                                      >  2:4111189/1‑150 (MQ=255)
aTAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTc                                                                                                                                                      <  2:2084888/150‑1 (MQ=255)
aTAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTc                                                                                                                                                      <  1:848172/150‑1 (MQ=255)
  accgacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCTCCCAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTTCATAACTCAACTCATTCCCCTCGCTATCCTTTTATTaaa                                                                                                                                                    >  1:3190052/3‑150 (MQ=255)
  aaCGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTaaa                                                                                                                                                    >  1:884749/1‑150 (MQ=255)
     gacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACt                                                                                                                                                 <  1:4177207/150‑1 (MQ=255)
     gacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACt                                                                                                                                                 >  1:1472768/1‑150 (MQ=255)
      acgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACtt                                                                                                                                                >  2:860299/1‑150 (MQ=255)
      acgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACtt                                                                                                                                                >  2:3601208/1‑150 (MQ=255)
         cTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCa                                                                                                                                             >  2:2497116/1‑150 (MQ=255)
           tGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCaaa                                                                                                                                           >  1:4188916/1‑150 (MQ=255)
            gAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAAt                                                                                                                                          <  2:4411841/150‑1 (MQ=255)
              aCATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTa                                                                                                                                        <  2:3229951/150‑1 (MQ=255)
              aCATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTa                                                                                                                                        <  2:3475907/150‑1 (MQ=255)
              aCATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTa                                                                                                                                        >  1:1008517/1‑150 (MQ=255)
                  gAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAat                                                                                                                                    >  2:4507983/1‑150 (MQ=255)
                   aGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAata                                                                                                                                   >  2:1052347/1‑150 (MQ=255)
                     gggTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATAtt                                                                                                                                 <  1:4205684/150‑1 (MQ=255)
                     gggTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCAGTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATAtt                                                                                                                                 <  2:1742168/150‑1 (MQ=255)
                      ggTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATAttt                                                                                                                                <  1:4355614/150‑1 (MQ=255)
                       gTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTa                                                                                                                               >  2:3501468/1‑150 (MQ=255)
                       gTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTa                                                                                                                               <  1:3446022/150‑1 (MQ=255)
                        tacATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTAt                                                                                                                              <  2:3901802/148‑1 (MQ=255)
                        tCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTAt                                                                                                                              <  2:1286240/150‑1 (MQ=255)
                           aTTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCtt                                                                                                                           <  1:128881/150‑1 (MQ=255)
                            tttAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCttt                                                                                                                          >  2:1405809/1‑150 (MQ=255)
                            tttAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCttt                                                                                                                          <  1:3106501/150‑1 (MQ=255)
                               aaCAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAt                                                                                                                       <  2:3438582/150‑1 (MQ=255)
                                 cAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAttt                                                                                                                     <  2:3068861/150‑1 (MQ=255)
                                  aGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAtttt                                                                                                                    >  2:3000925/1‑150 (MQ=255)
                                   gATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTg                                                                                                                   >  1:3975526/1‑150 (MQ=255)
                                    aTTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCAtttttc                                                                                                                  <  2:2948890/150‑3 (MQ=255)
                                     ttCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCCTTTTTCg                                                                                                                 <  1:2908571/150‑1 (MQ=255)
                                        aaCCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATc                                                                                                              <  1:940264/150‑1 (MQ=255)
                                        aaCCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATc                                                                                                              <  1:1999422/150‑1 (MQ=255)
                                          ccAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCaa                                                                                                            >  1:1883711/1‑150 (MQ=255)
                                             ggggCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAAt                                                                                                         <  1:70893/150‑1 (MQ=255)
                                             ggggCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAAt                                                                                                         >  1:2719926/1‑150 (MQ=255)
                                                        ggTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTaa                                                                                              >  2:2523571/1‑150 (MQ=255)
                                                         gTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTaaa                                                                                             <  2:1488455/150‑1 (MQ=255)
                                                          tAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAAt                                                                                            >  2:21638/1‑150 (MQ=255)
                                                             aGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCtt                                                                                         <  2:3367616/150‑1 (MQ=255)
                                                             aGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCtt                                                                                         <  1:1702915/150‑1 (MQ=255)
                                                              gCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                                                        <  2:732188/150‑1 (MQ=255)
                                                              gCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                                                        <  1:4440146/150‑1 (MQ=255)
                                                              gCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                                                        <  1:3343911/150‑1 (MQ=255)
                                                               cATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTc                                                                                       <  2:3294179/150‑1 (MQ=255)
                                                                aTCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCa                                                                                      <  1:4476316/150‑1 (MQ=255)
                                                                aTCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCa                                                                                      >  2:1291532/1‑150 (MQ=255)
                                                                  cACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAAt                                                                                    >  1:1404572/1‑150 (MQ=255)
                                                                   aCGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATc                                                                                   <  1:3519210/150‑1 (MQ=255)
                                                                        ccGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTgatt                                                                              <  1:2391774/150‑1 (MQ=255)
                                                                            acaAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTgattagat                                                                          <  1:3938711/150‑1 (MQ=255)
                                                                              aaGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgattagatta                                                                        >  1:2522567/1‑150 (MQ=255)
                                                                                ggAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAgg                                                                      >  1:2770746/1‑150 (MQ=255)
                                                                                ggAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAgg                                                                      >  1:2770839/1‑150 (MQ=255)
                                                                                     cGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGcc                                                                 >  2:4448719/1‑150 (MQ=255)
                                                                                     cGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGcc                                                                 <  1:518815/150‑1 (MQ=255)
                                                                                      ggTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCg                                                                <  1:485058/150‑1 (MQ=255)
                                                                                        tAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGtt                                                              <  1:1529799/150‑1 (MQ=255)
                                                                                         aGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGttt                                                             <  1:642010/150‑1 (MQ=255)
                                                                                          gTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTg                                                            <  2:2275167/150‑1 (MQ=255)
                                                                                          gTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTg                                                            <  2:1574951/150‑1 (MQ=255)
                                                                                            cACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGt                                                          >  1:2215440/1‑150 (MQ=255)
                                                                                               tGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGtaat                                                       >  1:3057474/1‑150 (MQ=255)
                                                                                                gCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGtaata                                                      >  2:3088331/1‑150 (MQ=255)
                                                                                                 cccGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGtaataa                                                     >  1:2590728/1‑150 (MQ=255)
                                                                                                   cGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACATTTAAATCTTTCAATCTGATTAGATTATGTTGCCGTTTGTTAATaaaa                                                   >  1:3616116/1‑150 (MQ=255)
                                                                                                     aTACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACa                                                 <  2:1755784/150‑1 (MQ=255)
                                                                                                     aTACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACa                                                 <  1:2791726/150‑1 (MQ=255)
                                                                                                        cGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATa                                              <  2:385158/150‑1 (MQ=255)
                                                                                                            cTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATc                                          <  2:1626795/150‑1 (MQ=255)
                                                                                                                     actcaactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgag                                 <  2:2599758/150‑1 (MQ=255)
                                                                                                                       tcaactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgagag                               <  2:1883711/150‑1 (MQ=255)
                                                                                                                       tcaactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgagag                               >  2:3127136/1‑150 (MQ=255)
                                                                                                                        caactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgagaga                              >  2:3847043/1‑150 (MQ=255)
                                                                                                                        caactcaTTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgagaga                              <  1:3543935/150‑1 (MQ=255)
                                                                                                                         aactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGaa                             >  2:418125/1‑150 (MQ=255)
                                                                                                                         aactcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGaa                             >  1:3031321/1‑150 (MQ=255)
                                                                                                                         aactcaTTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGaa                             <  2:3972666/150‑1 (MQ=255)
                                                                                                                           ctcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACa                           >  1:1099636/1‑150 (MQ=255)
                                                                                                                           ctcaTTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATATTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAAGAAATCCTGAAGGAGAGAACa                           >  1:1340937/1‑150 (MQ=255)
                                                                                                                            tcaTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACaa                          <  1:3333435/150‑1 (MQ=255)
                                                                                                                              aTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATg                        <  1:1214637/150‑1 (MQ=255)
                                                                                                                                tCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGAt                      >  2:3660085/1‑150 (MQ=255)
                                                                                                                                  cccTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAg                    <  1:1121521/150‑1 (MQ=255)
                                                                                                                                  cccTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTCATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAg                    <  1:614744/150‑1 (MQ=255)
                                                                                                                                    cTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAAACTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGaa                  >  2:2047076/1‑150 (MQ=255)
                                                                                                                                     tCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGaaa                 >  1:2172051/1‑150 (MQ=255)
                                                                                                                                       gCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAAcc               <  2:2831707/150‑1 (MQ=255)
                                                                                                                                       gCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAAcc               <  1:3216912/150‑1 (MQ=255)
                                                                                                                                         tATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCAt             >  2:54517/1‑150 (MQ=255)
                                                                                                                                         tATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCAt             >  1:146066/1‑150 (MQ=255)
                                                                                                                                          aTCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCAtt            <  1:4318380/150‑1 (MQ=255)
                                                                                                                                             cTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACt         >  2:4498233/1‑150 (MQ=255)
                                                                                                                                               tttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCa       >  2:1179430/1‑150 (MQ=255)
                                                                                                                                                ttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCAt      <  1:1304500/150‑1 (MQ=255)
                                                                                                                                                 tATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATg     <  1:4066887/150‑1 (MQ=255)
                                                                                                                                                 tATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATg     >  2:2646429/1‑150 (MQ=255)
                                                                                                                                                  aTTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATgg    >  2:2877204/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTG  >  CP009273/2818918‑2819215

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