Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,144 G→T 10.7% pseudogene (753/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,1440GT10.7% 134.5 / 2.7 66pseudogene (753/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (21/37);  new base T (4/4);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT  >  NC_000913/1465996‑1466290
                                                                                                                                                    |                                                                                                                                                  
gcgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCCACGGTAATAACGGGAAAGTg                                                                                                                                                    >  2:232340/1‑149 (MQ=255)
 cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGa                                                                                                                                                   >  1:46643/1‑149 (MQ=255)
   cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                                                                                                                                                 <  1:732038/149‑1 (MQ=255)
   cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                                                                                                                                                 <  2:826455/149‑1 (MQ=255)
   cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGAAAATGCCACAGTAAACAACAATGGCAAAACCACCGGTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                                                                                                                                                 <  2:492404/149‑1 (MQ=255)
     gATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGCCGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCa                                                                                                                                               >  2:855293/1‑149 (MQ=255)
     gATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCa                                                                                                                                               >  2:192455/1‑149 (MQ=255)
      aTGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGTTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGTAAATCAACGGTAATAACGGGAAAGTGATTCAg                                                                                                                                              <  2:334327/149‑1 (MQ=255)
      aTGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAg                                                                                                                                              <  2:165973/149‑1 (MQ=255)
           gcAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTATAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGAGATTCAGGACgg                                                                                                                                         >  1:618479/1‑149 (MQ=255)
                gAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATc                                                                                                                                    <  2:291793/149‑1 (MQ=255)
                 aaGTGGCAATATAATCACACCGGTGAGCTTGTTATTACGGGTGATAAAGCCACAGTAAACAACAATGGCAAAACCAACGTTGACGGTAAAGATCCCACCGGTACGGATAACAACGGTAATAACGGGAAAGTTATTCAGGACAGCGATcg                                                                                                                                   >  1:963366/1‑148 (MQ=255)
                   gTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                                                                                                                 <  1:881286/149‑1 (MQ=255)
                     ggCAATATAATCACAACGGTTAGCTTGTTATTACGGGTGATAATTCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                                                                                                               <  1:464435/149‑1 (MQ=255)
                     ggCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGCCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                                                                                                               <  1:541269/149‑1 (MQ=255)
                      gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                                                                                                              >  1:396665/1‑149 (MQ=255)
                         atatAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCa                                                                                                                           <  2:373571/149‑1 (MQ=255)
                           ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATACCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgc                                                                                                                         >  1:200780/1‑149 (MQ=255)
                            taATCACAAGGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTAAGGAATTCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAgcg                                                                                                                        >  1:543874/1‑149 (MQ=255)
                             aaTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                                                                                                       <  1:515204/149‑1 (MQ=255)
                                 acaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcggc                                                                                                                   >  2:515345/1‑149 (MQ=255)
                                  caACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTTCGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcggcg                                                                                                                  >  1:309027/1‑149 (MQ=255)
                                      ggTGAGGTTGGTAATACGGGTGATAATGCCACACTAAACAAGAAGCGCAAACACACGGGTGACGGTAAAGATTCCACCGGTACTGAAATCAACGGTAATACCGGGTAAGTGATTCAGTACGGCTATCTGGATGTCAGCGGCGGCGGTCa                                                                                                              <  1:199637/149‑1 (MQ=17)
                                        tGAGCTTGTTATTACGGGTGATAATGCCACAGTAACCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGCCGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                            >  2:209995/1‑149 (MQ=255)
                                            gttGTTATTACGTGTGATAATGCCACAGTAAACAACAATTGCAAAACCACAGTTGACGGATAAGATTCCTCCGGTTCGGAAATCAACGGTAATAACGGGAAACTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                                                                                        <  2:64687/148‑1 (MQ=17)
                                              tGTTATTACGTTTGATAATGCCACAGAAAACAACAATGGCAAAACCACCGTTGACGGAAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                                                                                      <  1:145846/149‑1 (MQ=255)
                                               gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                                                                                <  2:911644/138‑1 (MQ=255)
                                               gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                                                                                                                >  1:911644/1‑138 (MQ=255)
                                               gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGa                                                                                                     >  1:617063/1‑149 (MQ=255)
                                                         gTGATAATGCCACAGTAAACAACAATGGCAAACCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCCACGGTAATAACGGGAAAGTGATTCAGGTCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgct                                                                                           >  2:255947/1‑147 (MQ=255)
                                                         gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGt                                                                                           <  1:192455/149‑1 (MQ=255)
                                                             tAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACa                                                                                       <  1:763959/149‑1 (MQ=255)
                                                                 gCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgcgc                                                                                   >  1:391338/1‑149 (MQ=255)
                                                                  ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                            <  1:912899/123‑1 (MQ=255)
                                                                  ccACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                                                                                                            >  2:912899/1‑123 (MQ=255)
                                                                   cacaGTAAACCACAATGGCAAAACCACCGTGGACGGTAAAGATTCCACCGGTACGTAAATCAACGGTAATAACCGGAACGTGATTCAGGACGGCGATCTGGATCTCAGCGGCGGCGGTCACGGGATTGATATCACCGGTGACAGCGCGa                                                                                 <  1:478015/149‑1 (MQ=17)
                                                                     caGTAAACAACAATGGCAAAACCACCGTTGACGGTAACGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                                                                                                                   >  1:879330/1‑97 (MQ=255)
                                                                     caGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCt                                                                                                                                   <  2:879330/97‑1 (MQ=255)
                                                                      aGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACgg                                                                              <  2:983143/149‑1 (MQ=255)
                                                                       gTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                             <  2:7062/149‑1 (MQ=255)
                                                                           acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGCCAGCGCGACGGTGGAt                                                                         >  2:96009/1‑149 (MQ=255)
                                                                           acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                         >  1:536820/1‑149 (MQ=255)
                                                                           acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAGCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                         <  2:930625/149‑1 (MQ=255)
                                                                                aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGCAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt                                                                                                       <  1:817035/114‑1 (MQ=255)
                                                                                aaTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt                                                                                                       >  2:817035/1‑114 (MQ=12)
                                                                                 caggCACACCCACAGTTGACGGACAAGAAGGCACCGTCACGGAAATCAAAGGTATTAACGGGATAGTGATTCAGCACCGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                   <  2:697472/147‑1 (MQ=14)
                                                                                  tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                  >  2:84540/1‑149 (MQ=255)
                                                                                      aaaaCCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAg                                                                                      >  1:221880/1‑125 (MQ=17)
                                                                                             ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAcc                                                       <  1:124873/149‑1 (MQ=255)
                                                                                             ccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGAcc                                                       <  1:499203/149‑1 (MQ=255)
                                                                                                  gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAc                                                  <  1:515345/149‑1 (MQ=255)
                                                                                                     ggTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCAcccc                                               >  1:609392/1‑147 (MQ=255)
                                                                                                      gTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAt                                              <  2:46643/149‑1 (MQ=255)
                                                                                                         aaGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                           <  2:445040/149‑1 (MQ=40)
                                                                                                           gATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGa                                         <  2:369785/149‑1 (MQ=39)
                                                                                                            aTTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCCCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAg                                        >  1:776964/1‑149 (MQ=255)
                                                                                                              tCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGACATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                      <  2:617749/149‑1 (MQ=255)
                                                                                                                   cGGTACGGAAATCAACGGTACTAACGGGCAAGTGATTCAGGACGGCGAGCTGCATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGCGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGACACCGATCCGGAGGCCATgg                                 <  2:618479/149‑1 (MQ=11)
                                                                                                                   cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                 <  1:320612/149‑1 (MQ=40)
                                                                                                                   cGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                 <  1:106147/149‑1 (MQ=255)
                                                                                                                         ggAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCTGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATcc                           <  1:24899/149‑1 (MQ=255)
                                                                                                                         ggAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGAGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATcc                           >  2:929829/1‑149 (MQ=255)
                                                                                                                                    gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGt                <  1:1018565/149‑1 (MQ=255)
                                                                                                                                       ataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGAc             <  2:292014/149‑1 (MQ=255)
                                                                                                                                             ggAAAGGGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCTCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGcc       <  1:619524/149‑1 (MQ=255)
                                                                                                                                                acgTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATc    <  1:1020417/147‑1 (MQ=255)
                                                                                                                                                  gTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGt  <  1:232340/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGT  >  NC_000913/1465996‑1466290

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: