Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 1,670,899 C→T Q50* (CAG→TAG)  U712_08610 → Flagellar biosynthesis protein flhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,670,8990CT98.6% 230.4 / ‑4.0 70Q50* (CAG→TAG) U712_08610Flagellar biosynthesis protein flhA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (40/28);  total (42/28)
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.66e-01

CGTTTTAATAAGTGTTGTCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGCAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTAAAGTCGTCGAGACTTTCG  >  CP006881/1670772‑1671031
                                                                                                                               |                                                                                                                                    
cGTTTTAATAAGTGTTGTCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACAACAATGAATATGTAGGAAGCGCt                                                                                                                            >  7:93695/1‑138 (MQ=255)
   tttAATAAGTGTTGTCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAg                                                                                                                        >  6:57062/1‑139 (MQ=255)
            tgttgtCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGGTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGAtt                                                                                                               >  4:232328/1‑139 (MQ=255)
            tgttgtCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGAtt                                                                                                               >  5:363228/1‑139 (MQ=255)
            tgttgtCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGAtt                                                                                                               >  8:261936/1‑139 (MQ=255)
             gttgtCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGAttt                                                                                                              >  2:315962/1‑139 (MQ=255)
                gtCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTc                                                                                                           >  8:123697/1‑139 (MQ=255)
                 tCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTcc                                                                                                          >  5:16389/1‑139 (MQ=255)
                   cTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGt                                                                                                        >  4:143381/1‑139 (MQ=255)
                      aTTGTGGCGATGCTGGTGATTCCATTTCCTCCAAGCCGGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAAAGAAAATGTAGGAAGCGCTTCAGTTTTCGATTATTCCGTCAc                                                                                                     >  4:226149/1‑139 (MQ=255)
                           ggCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTtg                                                                                                >  2:29046/1‑139 (MQ=255)
                                tGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTtgct                                                                                              <  6:429583/136‑1 (MQ=255)
                                tGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTtgct                                                                                              >  5:429583/1‑136 (MQ=255)
                                tGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATCTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTgctgct                                                                                           >  4:206558/1‑139 (MQ=255)
                                tGCTGGTCATTCCAATTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTAGATCGTGCTACTCACCACAATGAATATGTAGGAAGCGCTTCAGATTTCGATTTTTCCGTCACTGTtgctgca                                                                                           >  8:352504/1‑138 (MQ=255)
                                  cTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCttt                                                                                         <  7:21550/139‑1 (MQ=255)
                                     gTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTg                                                                                       <  8:93695/138‑1 (MQ=255)
                                     gTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGa                                                                                      >  8:355679/1‑139 (MQ=255)
                                          tCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTa                                                                                 >  4:106196/1‑139 (MQ=255)
                                          tCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTa                                                                                 >  6:392833/1‑139 (MQ=255)
                                          tCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTa                                                                                 >  1:174244/1‑139 (MQ=255)
                                                            gTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTg                                                                                       >  7:334855/1‑115 (MQ=255)
                                                            gTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTg                                                                                       <  8:334855/115‑1 (MQ=255)
                                                             ttAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTc                                                                                                       <  3:30996/98‑1 (MQ=255)
                                                             ttAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTc                                                                                                       >  4:30996/1‑98 (MQ=255)
                                                             ttAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACg                                                              >  6:226284/1‑139 (MQ=255)
                                                             ttAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACg                                                              >  8:412212/1‑139 (MQ=255)
                                                             ttAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACg                                                              >  8:400377/1‑139 (MQ=255)
                                                               aaGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGt                                                             <  3:80640/138‑1 (MQ=255)
                                                                    ttttGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCCCTGTTGGTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCaa                                                        >  6:308545/1‑138 (MQ=255)
                                                                    ttttGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGATTGATCGTCTCAAc                                                       >  6:256493/1‑139 (MQ=255)
                                                                       tGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTAc                                                    >  2:43899/1‑139 (MQ=255)
                                                                       tGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCGGCGTTTGACCCTAGTTCGTTTGGGCTTGAACGTCTCAACTAc                                                    >  1:414878/1‑139 (MQ=255)
                                                                        gATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTAc                                                    <  5:57062/138‑1 (MQ=255)
                                                                        gATCATTATTAATATTTCTCTTGCGTTGATCGTGCGTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTAc                                                    <  4:318217/138‑1 (MQ=255)
                                                                          tCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCg                                                 <  7:355679/139‑1 (MQ=255)
                                                                          tCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCg                                                 <  5:392833/139‑1 (MQ=255)
                                                                                ttaATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCaa                                            <  1:43899/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACTATGAATTTGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTGTCAACTACCCGTTCAATTCTTTca                                    <  7:352504/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTcac                                   >  7:61525/1‑139 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTcac                                   >  5:287933/1‑139 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    <  7:123697/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    >  7:67893/1‑138 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    <  1:29046/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    <  6:16389/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    >  3:119212/1‑138 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    <  3:106196/138‑1 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    >  2:98351/1‑138 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGCAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTca                                    >  6:75981/1‑138 (MQ=255)
                                                                                        tctcTTGCGTTGATCGTGCTTCTCACCACAATGAATATGAAGCAAGCGCTTCAGTTTACGATTTTTCCGACACTGTTGCTTCTTTTGACGCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATACTTTca                                    >  4:319569/1‑138 (MQ=255)
                                                                                         ctctTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTcac                                   <  1:83009/138‑1 (MQ=255)
                                                                                           ctTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAAc                                                               >  3:89681/1‑108 (MQ=255)
                                                                                           ctTGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAAc                                                               <  4:89681/108‑1 (MQ=255)
                                                                                            ttGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTTCCCGTTCAATTCTTTcaca                                  <  3:259936/136‑1 (MQ=255)
                                                                                            ttGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTcaca                                  <  2:174244/136‑1 (MQ=255)
                                                                                            ttGCGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTcaca                                  <  4:24152/136‑1 (MQ=255)
                                                                                               cGTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTga                            <  7:400377/139‑1 (MQ=255)
                                                                                                gTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTgag                           <  8:67893/139‑1 (MQ=255)
                                                                                                gTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTgag                           <  5:283455/139‑1 (MQ=255)
                                                                                                gTTGATCGTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTgag                           <  3:143381/139‑1 (MQ=255)
                                                                                                       gTGCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTa                    <  6:363228/139‑1 (MQ=255)
                                                                                                         gCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTaaa                  >  3:37146/1‑139 (MQ=255)
                                                                                                         gCTTCTCACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTa                    >  3:74387/1‑137 (MQ=255)
                                                                                                            tctcACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGt                                                             >  3:19828/1‑93 (MQ=255)
                                                                                                            tctcACCACAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGt                                                             <  4:19828/93‑1 (MQ=255)
                                                                                                              tccccacAATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTAAcgtcgt             >  3:148326/4‑139 (MQ=255)
                                                                                                                  cacaATGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTAAAGTCGTCgag         >  3:286756/1‑139 (MQ=255)
                                                                                                                       tGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGtt                                               <  7:433790/96‑1 (MQ=255)
                                                                                                                       tGAATATGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGtt                                               >  8:433790/1‑96 (MQ=255)
                                                                                                                          atatGTAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTAAAGTCGTCGAGACTTTCg  >  3:451144/1‑138 (MQ=255)
                                                                                                                               |                                                                                                                                    
CGTTTTAATAAGTGTTGTCCTCATTGTGGCGATGCTGGTGATTCCATTTCCTCCATGGCTGTTAAGTATTTTGATCATTATTAATATTTCTCTTGCGTTGATCGTGCTTCTCACCACAATGAATATGCAGGAAGCGCTTCAGTTTTCGATTTTTCCGTCACTGTTGCTGCTTTTGACCCTATTTCGTTTGGGCTTGAACGTCTCAACTACCCGTTCAATTCTTTCACACGGTGAGGGAGGTAAAGTCGTCGAGACTTTCG  >  CP006881/1670772‑1671031

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