Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,268,694 T→A 7.1% L177Q (CTG→CAG)  kdsA → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,268,6940TA 7.1% 298.0 / 11.0 85L177Q (CTG→CAG) kdsA3‑deoxy‑D‑manno‑octulosonate 8‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (40/39);  new base A (3/3);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

TCAAGAAACCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCAC  >  NC_000913/1268574‑1268840
                                                                                                                        |                                                                                                                                                  
tCAAGAAACCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAt                                                                                                                         >  2:941320/1‑148 (MQ=255)
 cAAGAAACCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAtt                                                                                                                        >  2:1124597/1‑148 (MQ=255)
 cAAGAAACCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAtt                                                                                                                        <  1:1124597/148‑1 (MQ=255)
       aCCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCa                                                                                                                          >  3:563977/1‑140 (MQ=255)
       aCCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCa                                                                                                                          <  4:563977/140‑1 (MQ=255)
         cacaGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATgaagaa                                                                                                                >  4:678639/1‑148 (MQ=255)
         cacaGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATgaagaa                                                                                                                <  1:941320/148‑1 (MQ=255)
          acaGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGaaa                                                                                                               >  2:883284/1‑148 (MQ=255)
           caGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAg                                                                                                              >  1:73900/1‑148 (MQ=255)
             gTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGGTAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAgtg                                                                                                            <  4:499317/148‑1 (MQ=255)
             gTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAgtg                                                                                                            <  2:73900/148‑1 (MQ=255)
               ttGTCAGCCCGGGATAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCt                                                                                                                                                    >  4:1093598/1‑106 (MQ=255)
               ttGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCt                                                                                                                                                    <  3:1093598/106‑1 (MQ=255)
                  tCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAgtg                                                                                                            >  1:296239/1‑143 (MQ=255)
                  tCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAgtg                                                                                                            <  2:296239/143‑1 (MQ=255)
                       ccGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTa                                                                                                     <  1:824112/145‑1 (MQ=255)
                       ccGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTa                                                                                                     >  2:824112/1‑145 (MQ=255)
                               gATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTg                                                                                          >  2:635356/1‑148 (MQ=255)
                               gATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTg                                                                                          <  1:883284/148‑1 (MQ=255)
                                    gTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATga                                                                                                                    >  3:639156/1‑117 (MQ=255)
                                    gTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATga                                                                                                                    <  4:639156/117‑1 (MQ=255)
                                    gTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGAtttt                                                                                     <  3:678639/148‑1 (MQ=255)
                                         atCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              >  2:141755/1‑102 (MQ=255)
                                         atCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              <  1:141755/102‑1 (MQ=255)
                                         atCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTa                                                                                                                       >  3:480894/1‑109 (MQ=255)
                                         atCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTa                                                                                                                       <  4:480894/109‑1 (MQ=255)
                                              tGATAAATTCGAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTg                                                                                                        >  1:837289/1‑119 (MQ=255)
                                              tGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTg                                                                                                        <  2:837289/119‑1 (MQ=255)
                                                  aaaTTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATg                                                                                                                                       <  3:382751/84‑1 (MQ=255)
                                                  aaaTTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATg                                                                                                                                       >  4:382751/1‑84 (MQ=255)
                                                   aaTTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGt                                                                                           >  4:513650/1‑127 (MQ=255)
                                                   aaTTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGt                                                                                           <  3:513650/127‑1 (MQ=255)
                                                     ttCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTAt                                                                                                                      <  4:1108789/98‑1 (MQ=255)
                                                     ttCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTAt                                                                                                                      >  3:1108789/1‑98 (MQ=255)
                                                     ttCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCAGGTTGTcat                                                                                                                                           >  4:94814/1‑75 (MQ=255)
                                                     ttCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCAGGTTGTcat                                                                                                                                           <  3:94814/77‑3 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTg                                                                                          <  2:985409/125‑1 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTg                                                                                          >  1:985409/1‑125 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGAt                                                                                        >  1:707101/1‑127 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGAt                                                                                        <  2:707101/127‑1 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACt                                                                   >  3:571588/1‑148 (MQ=255)
                                                      tCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACt                                                                   >  3:438354/1‑148 (MQ=255)
                                                          agaagGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCt                                                                                                                                                    >  2:871197/1‑63 (MQ=255)
                                                          agaagGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCt                                                                                                                                                    <  1:871197/63‑1 (MQ=255)
                                                          agaagGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              <  3:900978/85‑1 (MQ=255)
                                                          agaagGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              >  4:900978/1‑85 (MQ=255)
                                                                  gcAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCAGGTTGTcat                                                                                                                                           >  4:96823/1‑62 (MQ=255)
                                                                  gcAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCAGGTTGTcat                                                                                                                                           <  3:96823/64‑3 (MQ=255)
                                                                   cAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCt                                                                                                         <  2:1005936/97‑1 (MQ=255)
                                                                   cAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCt                                                                                                         >  1:1005936/1‑97 (MQ=255)
                                                                     aCGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              <  3:1139023/74‑1 (MQ=255)
                                                                     aCGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTc                                                                                                                              >  4:1139023/1‑74 (MQ=255)
                                                                     aCGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAt                                                                                                                         >  2:163176/1‑79 (MQ=255)
                                                                     aCGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAt                                                                                                                         <  1:163176/79‑1 (MQ=255)
                                                                        aaaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGAt                                                                                        <  4:685500/109‑1 (MQ=255)
                                                                        aaaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGAt                                                                                        >  3:685500/1‑109 (MQ=255)
                                                                        aaaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGtt                                                 <  4:438354/148‑1 (MQ=255)
                                                                         aaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTc                                                                                                          <  3:850237/90‑1 (MQ=255)
                                                                         aaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTc                                                                                                          >  4:850237/1‑90 (MQ=255)
                                                                         aaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTc                                                                                                          <  1:718019/90‑1 (MQ=255)
                                                                         aaaaGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTc                                                                                                          >  2:718019/1‑90 (MQ=255)
                                                                            aGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAtt                                                                                                                        <  1:415335/73‑1 (MQ=255)
                                                                            aGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCAtt                                                                                                                        >  2:415335/1‑73 (MQ=255)
                                                                                aTTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGa                                                                                  >  4:940765/1‑107 (MQ=255)
                                                                                aTTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGa                                                                                  <  3:940765/107‑1 (MQ=255)
                                                                                     ttGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCTATGCCGCGa                                                       <  3:69499/129‑1 (MQ=255)
                                                                                     ttGCGATCGCGGTGCTAACTTCCGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGa                                                       >  4:69499/1‑129 (MQ=255)
                                                                                        cGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCggt                                 <  1:635356/148‑1 (MQ=255)
                                                                                         gATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTaa                                                                                                    >  2:490999/1‑80 (MQ=255)
                                                                                         gATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTaa                                                                                                    <  1:490999/80‑1 (MQ=255)
                                                                                          aTCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCa                                                                     >  3:649770/1‑110 (MQ=255)
                                                                                          aTCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCa                                                                     <  4:649770/110‑1 (MQ=255)
                                                                                           tCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCggtggt                              >  2:659494/1‑148 (MQ=255)
                                                                                                      aCTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCa                   <  1:659494/148‑1 (MQ=255)
                                                                                                         tCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCtt                                      >  2:1064530/1‑126 (MQ=255)
                                                                                                         tCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCtt                                      <  1:1064530/126‑1 (MQ=255)
                                                                                                               aTGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGt                                                  >  3:168315/1‑108 (MQ=255)
                                                                                                               aTGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGt                                                  <  4:168315/108‑1 (MQ=255)
                                                                                                                tGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTg                                                                              >  2:82848/1‑79 (MQ=255)
                                                                                                                tGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTg                                                                              <  1:82848/79‑1 (MQ=255)
                                                                                                                 gACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCaa                                                               >  3:8923/1‑93 (MQ=255)
                                                                                                                 gACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCaa                                                               <  4:8923/93‑1 (MQ=255)
                                                                                                                 gACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGgtcgt                           <  4:1009997/129‑1 (MQ=255)
                                                                                                                 gACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGgtcgt                           >  3:1009997/1‑129 (MQ=255)
                                                                                                                 gACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGc        >  4:997670/1‑148 (MQ=255)
                                                                                                                   cAACCAGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTgg              <  4:318694/140‑1 (MQ=255)
                                                                                                                   cAACCAGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTgg              >  3:318694/1‑140 (MQ=255)
                                                                                                                    aaCCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTgg     <  4:571588/148‑1 (MQ=255)
                                                                                                                       cTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCAc  >  4:569354/1‑148 (MQ=255)
                                                                                                                        |                                                                                                                                                  
TCAAGAAACCACAGTTTGTCAGCCCGGGACAGATGGGTAATATCGTTGATAAATTCAAAGAAGGCGGCAACGAAAAAGTGATTCTTTGCGATCGCGGTGCTAACTTCGGCTATGACAACCTGGTTGTCGATATGCTGGGCTTCAGCATTATGAAGAAAGTGTCTGGTAACTCGCCGGTGATTTTCGACGTGACCCACGCACTGCAATGCCGCGATCCGTTTGGCGCAGCTTCCGGTGGTCGTCGTGCTCAGGTGGCTGAGCTGGCAC  >  NC_000913/1268574‑1268840

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: