Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,104,750 T→C H112R (CAC→CGC)  wbbJ ← putative lipopolysaccharide biosynthesis O‑acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,104,7500TC98.5% 228.3 / ‑4.2 65H112R (CAC→CGC) wbbJputative lipopolysaccharide biosynthesis O‑acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (22/42);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGTGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGTCC  >  NC_000913/2104602‑2104895
                                                                                                                                                    |                                                                                                                                                 
ggCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGGCTGGAGGTATATTTGGCGAACTCATTGGATCAGAgcgc                                                                                                                                                 <  1:487064/151‑1 (MQ=255)
         tACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGa                                                                                                                                         >  2:667114/1‑150 (MQ=255)
           cGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAAcc                                                                                                                                      <  1:464071/151‑1 (MQ=255)
           cGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAAcc                                                                                                                                      <  1:388020/151‑1 (MQ=255)
                ccATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGAt                                                                                                                                 <  2:378495/151‑1 (MQ=255)
                    tACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATg                                                                                                                             >  2:744278/1‑151 (MQ=255)
                      ccAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGAt                                                                                                                           >  1:588832/1‑151 (MQ=255)
                       cAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATc                                                                                                                          <  1:702646/151‑1 (MQ=255)
                           aatTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGt                                                                                                                       <  2:2794/150‑1 (MQ=255)
                               gTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGtaata                                                                                                                   >  1:389584/1‑150 (MQ=255)
                                 tCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTaataaa                                                                                                                 <  1:200906/150‑1 (MQ=255)
                                  ccAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTaataaata                                                                                                               <  2:686608/151‑1 (MQ=255)
                                    aGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtact                                                                                                             >  2:290481/1‑151 (MQ=255)
                                      gCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtacttt                                                                                                           >  2:746849/1‑151 (MQ=255)
                                          aaCCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAAtactttactt                                                                                                       <  2:754903/151‑1 (MQ=255)
                                                cacaTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAAt                                                                                                  >  2:768858/1‑150 (MQ=255)
                                                 acaTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATaa                                                                                                <  1:735022/151‑1 (MQ=255)
                                                 acaTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATaa                                                                                                <  1:667641/151‑1 (MQ=255)
                                                 acaTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATaa                                                                                                <  1:245195/151‑1 (MQ=255)
                                                 acaTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATa                                                                                                 >  1:457241/1‑150 (MQ=255)
                                                   aTTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGc                                                                                              <  2:232425/151‑1 (MQ=255)
                                                        cACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTAt                                                                                          <  1:371511/150‑1 (MQ=255)
                                                          cccAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATccc                                                                                       <  1:6543/151‑1 (MQ=255)
                                                            caaccaaACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCg                                                                                      >  1:398839/1‑150 (MQ=255)
                                                            caaccaaACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGa                                                                                     <  1:584977/151‑1 (MQ=255)
                                                                 aaaCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTAt                                                                                >  2:428004/1‑151 (MQ=255)
                                                                  aaCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGGGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGACCCGTGATTATGAACGGTAATAAATACTTTACTTGCAATAAGCGTATCCTGACCTAt                                                                                >  1:471305/1‑150 (MQ=255)
                                                                   aCCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCg                                                                              <  1:398147/151‑1 (MQ=255)
                                                                               aTTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCaa                                                                  <  2:694638/151‑1 (MQ=255)
                                                                                  acaacaGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAAtt                                                                <  2:604613/150‑1 (MQ=255)
                                                                                      caGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGc                                                           >  1:68383/1‑151 (MQ=255)
                                                                                          tGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGat                                                        >  2:35971/1‑150 (MQ=255)
                                                                                          tGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGat                                                        >  2:425210/1‑150 (MQ=255)
                                                                                           gAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGatat                                                      >  1:422636/1‑151 (MQ=255)
                                                                                           gAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGatat                                                      <  1:367640/151‑1 (MQ=255)
                                                                                           gAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGatat                                                      <  2:136352/151‑1 (MQ=255)
                                                                                           gAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGatat                                                      <  1:251113/151‑1 (MQ=255)
                                                                                           gAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGata                                                       <  2:434973/150‑1 (MQ=255)
                                                                                            aaGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGa                                                              >  2:433885/1‑142 (MQ=255)
                                                                                            aaGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGa                                                              <  1:433834/142‑1 (MQ=255)
                                                                                                   cAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAg                                              <  1:619341/151‑1 (MQ=255)
                                                                                                     aCGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTc                                            <  2:621060/151‑1 (MQ=255)
                                                                                                         gcgcATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTa                                        <  2:290856/151‑1 (MQ=255)
                                                                                                           gcATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTaa                                       <  1:621461/150‑1 (MQ=255)
                                                                                                             aTGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACtt                                    <  2:63723/151‑1 (MQ=255)
                                                                                                                   ggAGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACAt                              <  1:455024/151‑1 (MQ=255)
                                                                                                                     aGGTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTa                            <  2:636061/151‑1 (MQ=255)
                                                                                                                      ggTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTAt                           <  2:598703/151‑1 (MQ=255)
                                                                                                                      ggTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTAt                           <  1:640381/151‑1 (MQ=255)
                                                                                                                       gTATATTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATc                          <  2:402097/151‑1 (MQ=255)
                                                                                                                           atTTGGCGAACTCATTGGATCAGAGTGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGa                       <  2:203688/150‑1 (MQ=255)
                                                                                                                           atTTGGCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGa                       >  1:517645/1‑150 (MQ=255)
                                                                                                                               ggCGAACTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGaaaaaa                  <  1:425566/151‑1 (MQ=255)
                                                                                                                                   aaCTCATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAAt              >  2:525454/1‑151 (MQ=255)
                                                                                                                                       cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCa            <  2:414292/149‑1 (MQ=255)
                                                                                                                                       cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCa            <  2:573451/149‑1 (MQ=255)
                                                                                                                                       cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACg          <  2:92703/151‑1 (MQ=255)
                                                                                                                                       cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACg          >  1:109693/1‑151 (MQ=255)
                                                                                                                                       cATTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACg          >  1:173456/1‑151 (MQ=255)
                                                                                                                                        aTTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACg          <  1:540754/150‑1 (MQ=255)
                                                                                                                                        aTTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGc         <  2:375632/151‑1 (MQ=255)
                                                                                                                                        aTTGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGc         <  1:252122/151‑1 (MQ=255)
                                                                                                                                          tGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGcc        <  1:236438/150‑1 (MQ=255)
                                                                                                                                          tGGATCAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGcc        >  1:135377/1‑150 (MQ=255)
                                                                                                                                               cAGAGCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGTcc  <  2:5529/151‑1 (MQ=255)
                                                                                                                                                agagCGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGTcc  >  1:28420/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGCCGACTACGACTCCATTACCAATAATTGTTCCAGGCAAAACCGTCACATTCTCACCCAACCAAACCCTCTGGCCAATTACAACAGCTGAAGATTCCAACGTGCGCATGTCTGGAGGTATATTTGGCGAACTCATTGGATCAGAGTGCTTAAAGGAACCGTGATTATGATCGGTAATAAATACTTTACTTGCAATAAGCGTATCCCGACCTATCGTAACGCTCTCAATTGAGGCGATATGAACATAGTCGTTAACTTGCACATTATCGGAAAAAAAAATCACGCCACGTCC  >  NC_000913/2104602‑2104895

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