Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,809,361 C→T 6.8% F274F (TTC→TTT) ‡ 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,3610CT 6.8% 264.1 / 7.0 109F274F (TTC→TTT) ‡waaA3‑deoxy‑D‑manno‑octulosonic‑acid transferase (KDO transferase)
Reads supporting (aligned to +/- strand):  ref base C (49/50);  new base T (3/4);  total (54/55)
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.80e-01

Only 100 of 115 total aligned reads displayed.

GCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTG  >  NC_000913/3809213‑3809509
                                                                                                                                                    |                                                                                                                                                    
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)
        ggggCACCACACCGCCCGTTAAGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGGAAGCGTTCCGGATga                                                                                                                                              <  2:622081/148‑2 (MQ=255)
          ggCACCACACCGCCCGGTATGGATTGCCACCAGCTCTCACGAAGGCGATGAGATTGTGGTGATCGCCGCAGATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCTTCCTGGTAGCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                            <  2:749656/149‑1 (MQ=255)
          ggCACCACACCGCCCGGTATGGATTGCCACCAGCACTAACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGa                                                                                                                                            <  1:106326/149‑1 (MQ=255)
            caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGATGGCGAAGAGAGTGTGGTGATCGCCGCACAACAGGCATAGTTACAGCAATTCCCGAATTTATTGCTAATCCTGGTACCCCGTCAACAGGAAAGCTTCAAGGAAGCGaat                                                                                                                                          >  1:765250/1‑147 (MQ=255)
            caccacACCGCCCGGTATGGATTGCCACCAGCACTCACGACGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGAtt                                                                                                                                          >  1:731656/1‑149 (MQ=255)
                  acCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCtt                                                                                                                                    >  2:472830/1‑149 (MQ=255)
                    cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGATTGTGGTGATCTCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGAGATTTACCTTGt                                                                                                                                  >  2:504503/1‑149 (MQ=255)
                    cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGt                                                                                                                                  >  2:480936/1‑149 (MQ=255)
                    cGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGt                                                                                                                                  >  1:981317/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)
                              ggATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAGCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTAATCCTCAAACCCCGTCAACCATAACGATTCACGAATTCTATTAACCTTGACAGCCacaaa                                                                                                                        >  2:680107/1‑145 (MQ=255)
                               gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCANCCGGAACNCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                                       <  1:356532/149‑1 (MQ=255)
                               gATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAGTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCGGGATGCGATTAACCTTGTCCGCCAGGCTg                                                                                                                       >  2:833585/1‑149 (MQ=255)
                                    ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTa                                                                                                                  >  1:494328/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)
                                      accaGCACTACCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGCATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTATCCTTGTCCGCCAGGCTGGACTAAg                                                                                                                >  2:777319/1‑149 (MQ=255)
                                       ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               <  2:269217/149‑1 (MQ=255)
                                       ccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               >  1:153406/1‑149 (MQ=255)
                                       ccaGCACTCACGAAGGCGAAGAGAGTGTGGGTATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGc                                                                                                               <  1:728160/149‑1 (MQ=255)
                                        caGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCt                                                                                                              >  1:413770/1‑149 (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)
                                                  gAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACg                                                                                                    >  1:221137/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)
                                                            agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAggg                                                                                          <  2:372382/149‑1 (MQ=255)
                                                             gtgtgtGGTGATCGCCGCAGATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCGGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAgggg                                                                                         <  2:983850/147‑1 (MQ=255)
                                                              agTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGa                                                                                        <  1:523928/149‑1 (MQ=255)
                                                                  tggtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGATCGCTTCCCGGATGCGATTAACCTTGTCCTCCAGGCTGGACTAAGCTATATCACACGCTCTACAGGGGAAGTc                                                                                    >  1:313575/1‑149 (MQ=255)
                                                                    gtgATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTccc                                                                                  <  2:12606/149‑1 (MQ=255)
                                                                      gATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTccccg                                                                                >  2:884894/1‑148 (MQ=255)
                                                                       aTCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTACCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                               <  1:560072/149‑1 (MQ=255)
                                                                       aTCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTACCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTc                                                                               <  1:154141/149‑1 (MQ=255)
                                                                          gccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACACCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCCccac                                                                            <  1:504503/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)
                                                                                             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)
                                                                                                        ccGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCCACAGCACGCAGGTTGTGGTGGGCGATACGAt                                              >  2:382512/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)
                                                                                                                   gCTCATCCTGGTACCCCGTCATCGGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGGTGTGGTTGGCGATACGATGGGCGAGTTGa                                   <  1:814382/149‑1 (MQ=255)
                                                                                                                   gCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGa                                   <  1:161009/149‑1 (MQ=255)
                                                                                                                   gCTCATCCTGGTACACCGTCATACGGAACGCTTCCCGGATGCGATTACACTTGTCCTCCAGGCTGGACTAAGCGATATAACACGCACTTCAGAGGAAGTCCCCTCAACCAGCACGCAGTGTGTGGTTGGCGATCCGATGGGCTAGTGGa                                   >  2:927021/1‑149 (MQ=255)
                                                                                                                    cTCATCCTGGTACCCCGTCATCCGGCACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTATGCTATATCACACGCTCTTCAGGGGAAGTCCCCTGCACCATCACGCAGGTTGTGGTGGGCGATACGATGTGAGAGTTGAt                                  >  1:583607/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)
                                                                                                                                  ccGTCATCCGGAACGCTTGCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATTCGATGGGCGAGTTGATGTTACTGTATGGCa                    >  1:717165/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)
                                                                                                                                                    cccGGATGCGATTTACCTTGTCCGCCAGGCTGGACTAAGCTATATCACCCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCTTTGCCGATCTCGCCTttg  >  1:609837/1‑149 (MQ=255)
                                                                                                                                                    cccGGATGCGATTAACCTTGTCCGCCCGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCATCAACCAGCACGCAGGTTGTGGTTGGCGATACGATGGTCGAGTTGATGTTACTGTATGGCTTTGCCGATCTCGCCTttg  <  2:153406/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTG  >  NC_000913/3809213‑3809509

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: