Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,354 A→G 5.3% E272G (GAA→GGA)  waaA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,3540AG 5.3% 307.4 / 6.5 115E272G (GAA→GGA) waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base A (54/52);  new base G (3/3);  total (59/56)
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.31e-01

Only 100 of 117 total aligned reads displayed.

CGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTC  >  NC_000913/3809207‑3809502
                                                                                                                                                   |                                                                                                                                                    
cGCTGCGCCGCCAGTTGGCACGACACCTCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGaa                                                                                                                                                     <  1:332783/149‑1 (MQ=255)
 gCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGAGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAAc                                                                                                                                                    <  1:736824/149‑1 (MQ=255)
   tGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGc                                                                                                                                                  <  1:641879/149‑1 (MQ=255)
      gccgccAGTGGGCACCACACCGCCGGGTATGGATTGCCACCAGCACTCACGACAGCGAAGAGAGTGTGGTGATCGCCGAACATCAGGAAATGTTACAGCAATTCCCGAATTTTTTGCTCATCCTGGTACCCAGTCATCCGGAACGCTTc                                                                                                                                               >  1:578979/1‑149 (MQ=255)
        cgccAGTGGGCACCACCCCGCCCGGTATGGATTGCCCCCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATAGTTACCGCAATTCCCGCATTTATTGGTCATCCTGGTACACCGTCTTCCGGAACGCTTcgc                                                                                                                                             >  1:602082/1‑147 (MQ=255)
        cgccAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGGGGTGCTCGCCGGGCATCAGGCATTGTTACAGCAATTCCCGAGTTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTccc                                                                                                                                             >  1:234619/1‑149 (MQ=255)
        cgccAGTGGGCACCACACCGCCCGGTAGGGATTGCCACCAGCACTCCCGAAGGCGAAGCGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCACTTCCCGAATGTATTGCTCCTCCTGGTACCCCGTCATCCGGAACGCTTccc                                                                                                                                             >  2:1013447/1‑149 (MQ=255)
           cAGTGGGCACCCCACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTGACAGCACTTCCCGAATTTATTGCTCATCCAGGTACCCCGTCATCCGGCTCGCTTCCCggt                                                                                                                                          >  1:31339/1‑148 (MQ=255)
              tGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATga                                                                                                                                       <  1:21785/149‑2 (MQ=255)
              tGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGc                                                                                                                                       <  1:283083/149‑1 (MQ=255)
              ggggCACCACACCGCCCGTTAAGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATga                                                                                                                                       <  2:622081/148‑2 (MQ=255)
                ggCACCCCACCGCCCGGTATGGATTGCCACCAGCACTCACGCAGGCGAAGAGAGTGTGGGGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGATTTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                     <  2:661619/149‑1 (MQ=255)
                ggCACCACACCGCCCGGTATGGATTGCCACCAGCTCTCACGAAGGCGATGAGATTGTGGTGATCGCCGCAGATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCTTCCTGGTAGCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                     <  2:749656/149‑1 (MQ=255)
                ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTAACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                     <  1:106326/149‑1 (MQ=255)
                 gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGACGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATgacgg                                                                                                                                    <  1:16125/149‑5 (MQ=255)
                  caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGATGGCGAAGAGAGTGTGGTGATCGCCGCACAACAGGCATAGTTACAGCAATTCCCGAATTTATTGCTAATCCTGGTACCCCGTCAACAGGAAAGCTTCAAGGAAGCGaat                                                                                                                                   >  1:765250/1‑147 (MQ=255)
                  caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGACGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGAtt                                                                                                                                   >  1:731656/1‑149 (MQ=255)
                    caccacCGCCCGGTATGGCTTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTCCAGCAATTCACTAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTaa                                                                                                                                 <  1:67650/146‑1 (MQ=255)
                         ccGCCCGGTATGGATTGCCACCAGCACTCACGACGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTg                                                                                                                            <  1:458377/149‑1 (MQ=255)
                          cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGATTGTGGTGATCTCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGAGATTTACCTTGt                                                                                                                           >  2:504503/1‑149 (MQ=255)
                          cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGt                                                                                                                           >  1:981317/1‑149 (MQ=255)
                          cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGt                                                                                                                           >  2:480936/1‑149 (MQ=255)
                           gCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTc                                                                                                                          <  2:215994/149‑1 (MQ=255)
                            cccGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGCTCGCCGCGCATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTCTcc                                                                                                                         >  1:135813/1‑149 (MQ=255)
                               ggTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGcc                                                                                                                      <  1:187282/149‑1 (MQ=255)
                                gctaggAATGACACCATCACTCACGATGGCGTAGAGAGTTTCTGGATCGCCGCACAGCAGCCATTGTTACAGCAATTCCCGAATTTCTTGCTCCTCCTGGTACCCCGGCATCGGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCa                                                                                                                     <  2:748281/145‑1 (MQ=255)
                                 tATGGATTGCCACCAGCACTCACGAAGTCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAg                                                                                                                    <  2:697769/149‑1 (MQ=255)
                                     gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCANCCGGAACNCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                                <  1:356532/149‑1 (MQ=255)
                                     gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAGTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCGGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                                >  2:833585/1‑149 (MQ=255)
                                          ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCATCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTa                                                                                                           >  2:949395/1‑149 (MQ=255)
                                            cccaGCACTCACGAAGGCGAAGAGAGTGTGGTGCTCGCGGCACATCCGCCATGGGTACAGCAAGTCCCGACGAGCTTGCGCATCCTGGTCCCCCGTCAGCCGGAACATGTCGGTGTTGCTATTAACCTTGTGCGCGAGGCTGGAGTAAg                                                                                                         >  2:673836/2‑149 (MQ=255)
                                            accaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACACCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAg                                                                                                         <  1:437796/149‑1 (MQ=255)
                                            accaGCACTCACGAAGGCGAAGAGAGTGTGGTGAACGCCGCACATAATGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAg                                                                                                         <  2:882592/149‑1 (MQ=255)
                                             ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                        >  1:153406/1‑149 (MQ=255)
                                             ccaGCACTCACGAAGGCGAAGAGAGTGTGGGTATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                        <  1:728160/149‑1 (MQ=255)
                                              caGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                       <  1:18265/149‑1 (MQ=255)
                                              caGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                       <  1:711708/149‑1 (MQ=255)
                                                gCACTCACGAAGGCGAATAGAGTGTGGTGATGGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGACCGCTTCCCGTATGCGTTTACCCTTGTCCGCCAGGCTGGACTAAGCtat                                                                                                     <  2:868024/149‑1 (MQ=255)
                                                 cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCtata                                                                                                    <  1:984317/149‑1 (MQ=255)
                                                   cTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTTGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATc                                                                                                  <  1:473106/149‑1 (MQ=255)
                                                    tCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATca                                                                                                 >  2:224517/1‑149 (MQ=255)
                                                     cACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcac                                                                                                <  1:886672/149‑1 (MQ=255)
                                                      aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcacc                                                                                               >  2:798993/1‑148 (MQ=255)
                                                      aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATcaca                                                                                               <  2:392619/149‑1 (MQ=255)
                                                        gAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGACCTCGTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTACGCTATATCACACg                                                                                             >  1:745368/1‑149 (MQ=255)
                                                        gAAGGCGAAGAGAGTGGGGTGATCGCCGCACATCAGGCATGGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACg                                                                                             <  2:377037/149‑1 (MQ=255)
                                                             cGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCtt                                                                                        <  2:307909/149‑1 (MQ=255)
                                                              gAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGTTACCCCGTCATCCGGAACGCTGCCCGGCTGCGATTTACCTTGTCCGCCAGGCTGGACTTAGCTATATCACACGCTCTTc                                                                                       >  2:314194/1‑149 (MQ=255)
                                                              gAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGCCTAAGCTATATCACACGCTCTTc                                                                                       >  1:182393/1‑149 (MQ=255)
                                                              gAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTc                                                                                       >  1:809678/1‑149 (MQ=255)
                                                                  tgagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCACTTCCCGAATTTATTGCTCAACCTGGTACCCCGTCATCCGGAACGCTTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTATTCAggg                                                                                   <  2:283156/148‑1 (MQ=255)
                                                                  agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                   <  2:392717/149‑1 (MQ=255)
                                                                  agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                   <  2:372382/149‑1 (MQ=255)
                                                                   gtgtgtGGTGATCGCCGCAGATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCGGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAgggg                                                                                  <  2:983850/147‑1 (MQ=255)
                                                                    agTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGa                                                                                 <  1:523928/149‑1 (MQ=255)
                                                                       gtggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGt                                                                              <  1:543844/149‑1 (MQ=255)
                                                                        tggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGATCGCTTCCCGGATGCGATTAACCTTGTCCTCCAGGCTGGACTAAGCTATATCACACGCTCTACAGGGGAAGTc                                                                             >  1:313575/1‑149 (MQ=255)
                                                                             aTCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTACCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                        <  1:560072/149‑1 (MQ=255)
                                                                             aTCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTACCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                        <  1:154141/149‑1 (MQ=255)
                                                                                gccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCCccac                                                                     <  1:504503/149‑1 (MQ=255)
                                                                                  cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCCTCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTccacca                                                                   <  2:503723/149‑1 (MQ=255)
                                                                                    cacaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCTTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGGCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGc                                                                 <  2:731656/149‑1 (MQ=255)
                                                                                     acaTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCCCACGCTCTTCAGGGGAAGGCCCCTCCACCAGcc                                                                >  2:510867/1‑148 (MQ=255)
                                                                                       aTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTACCCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACg                                                              >  2:419025/1‑149 (MQ=255)
                                                                                       aTCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGCCTAAGCTATATCACGCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACg                                                              >  1:945343/1‑149 (MQ=255)
                                                                                        tCAGGCATTGTTACAGCAATTCCCGAATTTATTGTTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGc                                                             >  1:890920/1‑149 (MQ=255)
                                                                                        tCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGAACCCCGTCATCCGGAACGCTTACCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGGTTTATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGc                                                             >  2:616392/1‑149 (MQ=255)
                                                                                         cAGGCATTGTTACAGCAATACCCGAATTTATAGCTCATCCTGGTACGCCGTCATCCGTAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGct                                                            <  1:52650/149‑2 (MQ=255)
                                                                                           ggCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGTTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCAAACGCTCTTCACGGGATGTCCCCTCCACACGCACGCAgg                                                          >  2:472848/1‑149 (MQ=255)
                                                                                                 gTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCTACAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTgg                                                    >  2:244486/1‑149 (MQ=255)
                                                                                                 gTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTgg                                                    >  2:248194/1‑149 (MQ=255)
                                                                                                  ttACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTACCCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGt                                                   >  2:161009/1‑149 (MQ=255)
                                                                                                   tACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGtt                                                  <  1:480936/149‑1 (MQ=255)
                                                                                                   tACAGCAATTCCCGAATTTATTGCTCATCCTCGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGGTGTGGtt                                                  <  2:66417/149‑1 (MQ=255)
                                                                                                    aCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTg                                                 >  2:813510/1‑149 (MQ=255)
                                                                                                    aCAGCAATTCCCGAATATATTGCTCATCCTGGTACCCCGGCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTg                                                 <  2:534498/149‑1 (MQ=255)
                                                                                                      cgCAAATCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCGTCCCGGATGCGATTAACCTGGTCCGCCAGGCTGGACTACGCTATATCCCACGCTCTTCAGGGGCAGTCCCCTCCCCCCGCCCGCAGGTTGTGGTTGGc                                               <  2:970587/148‑1 (MQ=255)
                                                                                                      aGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGc                                               <  2:118800/149‑1 (MQ=255)
                                                                                                      aGCAATTCCCGAATTTATTGCTCATCCTCGTACCCCGTCATCCGGAACCCTTCCCGGAACAGAATAACCTTGTCCGACAGGCTGGACTAAGCAAAATCACAAGCTCTTCTGCGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGATGGc                                               >  1:589444/1‑149 (MQ=255)
                                                                                                       gCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCCCCAGCACGCAGGTTGGGGTTGGCg                                              >  2:775734/1‑149 (MQ=255)
                                                                                                        cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCGCCTCCACCAGCACGGAGGTTGTGGTTGGCGa                                             >  1:476747/1‑149 (MQ=255)
                                                                                                              ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTGGTGGTTGGCGAGACGAt                                       >  2:636090/1‑149 (MQ=255)
                                                                                                              ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTCTCCGCAAGGCTGGACTAAGCTAGATCACACGCTCTTCAGGGGAAGACCCCTCCAGCAGCACGCAGGGGGTGGTAGGCGAGACGAt                                       >  1:669703/1‑149 (MQ=255)
                                                                                                               cGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATg                                      >  1:838818/1‑149 (MQ=255)
                                                                                                                gAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTACCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGGGGGTGGCGATAGGATgg                                     >  2:464770/1‑149 (MQ=255)
                                                                                                                 aaTTAAATGCGCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGTTACGATggg                                    <  1:1033250/149‑1 (MQ=255)
                                                                                                                         gCTGATCCTGGTACCCCGTCCTCCGGAACGCTTCCCGGATGCGATTAACCTTGACCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGTGAAGTCTCATCCACCAGCTCGGAGGTTGTGGTTGGCGATACGATGTGCGAGTTGa                            <  2:765250/149‑1 (MQ=255)
                                                                                                                         gCTCATCCTGGTACACCGTCATACGGAACGCTTCCCGGATGCGATTACACTTGTCCTCCAGGCTGGACTAAGCGATATAACACGCACTTCAGAGGAAGTCCCCTCAACCAGCACGCAGTGTGTGGTTGGCGATCCGATGGGCTAGTGGa                            >  2:927021/1‑149 (MQ=255)
                                                                                                                          cTCATCCTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGATGAGGTTGGCGATACGATGGGCGAGTTGAt                           >  2:367828/1‑149 (MQ=255)
                                                                                                                            cATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGt                         <  1:126565/149‑1 (MQ=255)
                                                                                                                              tgcTGGGACCCCCGCCTCCGGAACGGTTGCCGGTTGGGGTTAACCTTTTGCGCCATGCTGGTCTTAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTa                       <  2:981317/147‑1 (MQ=255)
                                                                                                                                cTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACt                     <  1:454879/149‑1 (MQ=255)
                                                                                                                                     aCCCCGTCATCCTGTAGGCTTCCCGGATGCGATTAACCTTGTGAGGCAGGCTTGACTAAGGTATATCACTCGCTCTTCAGGGGAATTCACCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATg                <  1:846376/149‑1 (MQ=255)
                                                                                                                                      ccccGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGACCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCTGCACGCAGCTTGTGGTTGGCGATACGATGGGCGCCTTGATGTTACTGTATgg               >  1:844820/1‑149 (MQ=255)
                                                                                                                                         tcgcATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCAt            >  1:868466/4‑149 (MQ=255)
                                                                                                                                          cgcATCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTACTGTATGGCATTg          >  2:777444/3‑149 (MQ=255)
                                                                                                                                             aTCCGGGAAGCGTTCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGcc        >  2:262695/1‑149 (MQ=255)
                                                                                                                                              tCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCCCACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCg       >  1:202809/1‑149 (MQ=255)
                                                                                                                                                cGGAAAGCTTCCCGGATGCGATTATCCTTGAGCGCCCGGCTGGACTAAGCTATAGCACCCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATTCGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAt     <  2:442811/149‑1 (MQ=255)
                                                                                                                                                   aaCGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTACCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGAtctc  <  2:413770/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTC  >  NC_000913/3809207‑3809502

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: