Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 553,885 G→A 7.0% T19I (ACC→ATC)  lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913553,8850GA 7.0% 319.4 / 9.9 99T19I (ACC→ATC) lpxHUDP‑2,3‑diacylglucosamine pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (51/41);  new base A (4/3);  total (55/44)
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.98e-01

Only 100 of 101 total aligned reads displayed.

TCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGA  >  NC_000913/553738‑554031
                                                                                                                                                   |                                                                                                                                                  
tCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGt                                                                                                                                                   <  3:748679/149‑1 (MQ=255)
   gACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                <  3:6333/149‑1 (MQ=255)
      acCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGc                                                                                                                                             >  3:521413/1‑149 (MQ=255)
           tttGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGt                                                                                                                                         <  1:722872/148‑1 (MQ=255)
             tGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGttct                                                                                                                                      <  4:307333/149‑1 (MQ=255)
                tCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTcc                                                                                                                                   >  3:145402/1‑149 (MQ=255)
                tCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTcc                                                                                                                                   <  4:490468/149‑1 (MQ=255)
                  gccgccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGc                                                                                                                                             >  3:73906/1‑137 (MQ=255)
                  gccgccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGc                                                                                                                                             <  4:73906/137‑1 (MQ=255)
                     gccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGc                                                                                                                               <  1:825517/148‑1 (MQ=255)
                          cATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGAt                                                                                                                         >  4:817069/1‑149 (MQ=255)
                             cTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGaa                                                                                                                      <  3:780688/149‑1 (MQ=255)
                              ttGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                <  3:616870/122‑1 (MQ=255)
                              ttGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGAt                                                                                                                                                >  4:616870/1‑122 (MQ=255)
                                 cGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                           <  4:545940/124‑1 (MQ=255)
                                 cGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCg                                                                                                                                           >  3:545940/1‑124 (MQ=255)
                                   aTGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTg                                                                                                                <  3:817069/149‑1 (MQ=255)
                                    tGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTg                                                                                                                <  2:860456/148‑1 (MQ=255)
                                      gAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                   >  2:899677/1‑143 (MQ=255)
                                      gAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGAt                                                                                                                   <  1:899677/143‑1 (MQ=255)
                                          gggTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaaa                                                                                                         <  3:359081/149‑1 (MQ=255)
                                            gTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCaa                                                                                                             >  1:532583/1‑143 (MQ=255)
                                            gTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCaa                                                                                                             <  2:532583/143‑1 (MQ=255)
                                            gTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                       <  4:659134/149‑1 (MQ=255)
                                             ttGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAagag                                                                                                      >  4:714071/1‑149 (MQ=255)
                                             ttGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAagag                                                                                                      <  3:612461/149‑1 (MQ=255)
                                                      tcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCa                                                                                              <  2:177974/148‑1 (MQ=255)
                                                        gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACg                                                                                            <  2:887481/148‑1 (MQ=255)
                                                         tcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGtt                                                                                                                                        <  4:172134/103‑1 (MQ=255)
                                                         tcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGATCACCGCCGGtt                                                                                                                                        >  3:172134/1‑103 (MQ=255)
                                                          cgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGAtt                                                                                         <  3:596303/149‑1 (MQ=255)
                                                            ccAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTa                                                                                        >  1:166845/1‑148 (MQ=255)
                                                             cAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACt                                                                                      >  4:53800/1‑149 (MQ=255)
                                                                   aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaa                                                                                                          >  3:564908/1‑123 (MQ=255)
                                                                   aTGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATaa                                                                                                          <  4:564908/123‑1 (MQ=255)
                                                                    tGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAAc                                                                               <  4:145402/149‑1 (MQ=255)
                                                                       ttCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGt                                                                            >  4:381559/1‑149 (MQ=255)
                                                                       ttCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGt                                                                            >  4:832816/1‑149 (MQ=255)
                                                                          aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                         >  3:642436/1‑149 (MQ=255)
                                                                          aaaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCAc                                                                         >  3:40177/1‑149 (MQ=255)
                                                                           aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                   <  2:745240/122‑1 (MQ=255)
                                                                           aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAgtgt                                                                                                   >  1:745240/1‑122 (MQ=255)
                                                                           aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACg                                                                        >  4:702771/1‑149 (MQ=255)
                                                                           aaCAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACg                                                                        >  3:133781/1‑149 (MQ=255)
                                                                               gATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCttt                                                                    <  4:87/149‑1 (MQ=255)
                                                                                   gCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAAt                                                                <  4:133781/149‑1 (MQ=255)
                                                                                     cAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATg                                                               >  1:124564/1‑148 (MQ=255)
                                                                                      aaGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGAt                                                             >  4:373835/1‑149 (MQ=255)
                                                                                        gAATATACAGCGCGTCGGCCTTGCGGGCTTTCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGTGTGTCGCCACGATTACTCTCTaa                                                                                <  2:742444/128‑1 (MQ=255)
                                                                                        gAATATACAGCGCGTCGGCCTTGCGGGCTTTCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTaa                                                                                >  1:742444/1‑128 (MQ=255)
                                                                                           tataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTCCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTa                                                                                 <  2:260056/124‑1 (MQ=255)
                                                                                           tataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTCCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTa                                                                                 >  1:260056/1‑124 (MQ=255)
                                                                                            ataCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGtc                                                       >  3:406978/1‑149 (MQ=255)
                                                                                             taCAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATTATaaaggc                                                       >  1:596315/1‑144 (MQ=255)
                                                                                               cAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAATATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGtcttc                                                    <  4:521413/149‑1 (MQ=255)
                                                                                                  cgcgTCGGCCTTGGGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTcttctt                                                  >  1:685225/1‑148 (MQ=255)
                                                                                                  cgcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATg                                                               >  1:789546/1‑135 (MQ=255)
                                                                                                  cgcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATg                                                               <  2:789546/135‑1 (MQ=255)
                                                                                                   gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGCTTACGTCTTCtttt                                                <  3:67167/149‑1 (MQ=255)
                                                                                                   gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                <  3:447568/149‑1 (MQ=255)
                                                                                                   gcgTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                <  3:381559/149‑1 (MQ=255)
                                                                                                         gCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACg                                          >  4:809834/1‑149 (MQ=255)
                                                                                                           cTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGATATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTc                                        >  4:522618/1‑149 (MQ=255)
                                                                                                            ttGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTcc                                       >  3:512465/1‑149 (MQ=255)
                                                                                                            ttGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTcc                                       >  3:88349/1‑149 (MQ=255)
                                                                                                                 ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCt                                      >  2:198949/1‑145 (MQ=255)
                                                                                                                 ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCt                                      <  1:198949/145‑1 (MQ=255)
                                                                                                                 ggCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTg                                  <  3:665475/149‑1 (MQ=255)
                                                                                                                   cTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGc                                 <  2:685225/148‑1 (MQ=255)
                                                                                                                       ccccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATaac                            >  3:469396/1‑147 (MQ=255)
                                                                                                                         ccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                           >  2:552415/1‑148 (MQ=255)
                                                                                                                         ccGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCg                           >  2:320062/1‑148 (MQ=255)
                                                                                                                          cGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGc                          >  2:530514/1‑148 (MQ=255)
                                                                                                                          cGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGc                          >  2:68307/1‑148 (MQ=255)
                                                                                                                          cGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGAGGCAGATCTTCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGc                          >  2:601685/1‑148 (MQ=255)
                                                                                                                           gCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                        >  4:147856/1‑149 (MQ=255)
                                                                                                                            cTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACg                                                                        >  1:469632/1‑100 (MQ=255)
                                                                                                                            cTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACg                                                                        <  2:469632/100‑1 (MQ=255)
                                                                                                                            cTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTa                        >  2:867572/1‑148 (MQ=255)
                                                                                                                              aaaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACg                      >  2:615999/1‑148 (MQ=255)
                                                                                                                               aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                       <  4:73198/66‑1 (MQ=255)
                                                                                                                               aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAaga                                                                                                       >  3:73198/1‑66 (MQ=255)
                                                                                                                               aaaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGa                     >  2:846122/1‑148 (MQ=255)
                                                                                                                                aaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                <  1:303186/120‑1 (MQ=255)
                                                                                                                                aaaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCtttt                                                >  2:303186/1‑120 (MQ=255)
                                                                                                                                 aaaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCg                  >  4:509657/1‑149 (MQ=255)
                                                                                                                                  aaCGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCg                  >  2:849580/1‑148 (MQ=255)
                                                                                                                                    cGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGt                >  2:685761/1‑148 (MQ=255)
                                                                                                                                      cAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTg              <  2:236632/148‑1 (MQ=255)
                                                                                                                                      cAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTg              <  1:530514/148‑1 (MQ=255)
                                                                                                                                       gaagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCaa                                                                 >  4:41390/5‑96 (MQ=38)
                                                                                                                                       gaagaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCaa                                                                 <  3:41390/92‑1 (MQ=38)
                                                                                                                                       aGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa            >  4:5515/1‑149 (MQ=255)
                                                                                                                                        gAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCaa           >  3:893706/1‑149 (MQ=255)
                                                                                                                                        gAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCa            <  1:849580/148‑1 (MQ=255)
                                                                                                                                          gaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAaca         >  4:140937/2‑149 (MQ=255)
                                                                                                                                          gaaCCGGCGATCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAaca         >  3:552483/2‑149 (MQ=255)
                                                                                                                                          aaaCCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAac          >  1:466029/1‑148 (MQ=255)
                                                                                                                                              cggcgaTCACCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCt      >  2:295420/1‑148 (MQ=255)
                                                                                                                                                 cggTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGa  <  3:832816/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAAAAACGCAGAAAACCGGCGGTGATCGCCGGTTCTTCCACGCAGAGATGAAGATCTGCAATAAAGAGTGTCGCCACGATTACTCGCTAACGGTCACGCTTTCAATGATAACGTCTTCTTTTGGCACGTCCTGGTGCATACCGCTACGACCGGTTGCAACACCTTTGA  >  NC_000913/553738‑554031

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