Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,055 G→T 5.5% pseudogene (664/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,0550GT 5.5% 296.7 / 3.7 90pseudogene (664/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base T (2/4);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGA  >  NC_000913/1465908‑1466203
                                                                                                                                                   |                                                                                                                                                    
tctcTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTGCGGGTGa                                                                                                                                                     >  8:72772/1‑149 (MQ=255)
    tGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGa                                                                                                                                                     <  4:196479/145‑1 (MQ=255)
    tGCAAGATAGCAATGGTCGTAAGGCAACGATTAATATTTGGCAGATTGATGAAGCGAAAAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGCTGAGCTTGTTATTACGGGTGATAAt                                                                                                                                                 >  5:260796/1‑149 (MQ=255)
     gCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATg                                                                                                                                                >  5:305188/1‑149 (MQ=255)
      cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATaa                                                                                                                                                  <  4:239368/146‑1 (MQ=255)
      cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATaa                                                                                                                                                  <  7:33395/146‑1 (MQ=255)
      cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATaa                                                                                                                                                  <  2:178936/146‑1 (MQ=255)
      cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGc                                                                                                                                               >  3:266313/1‑149 (MQ=255)
      cAAGATAGCAATGGGCGGAAGGCCACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATa                                                                                                                                                   <  5:391104/145‑1 (MQ=255)
         gATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCcac                                                                                                                                            <  3:381488/149‑1 (MQ=255)
         gATAGCAATGGACGTAATGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCAATTGAAGTGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCca                                                                                                                                             >  8:210119/1‑148 (MQ=255)
            aGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGTGTGATAATGCCACAg                                                                                                                                          >  2:118250/1‑148 (MQ=255)
                 tGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATaa                                                                                                                                                  <  8:231922/135‑1 (MQ=255)
                 tGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATaa                                                                                                                                                  >  7:231922/1‑135 (MQ=255)
                 tGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAaac                                                                                                                                     >  8:136292/1‑148 (MQ=255)
                   gTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAaacaa                                                                                                                                   <  7:121488/148‑1 (MQ=255)
                    tCGTAAGGCAAGGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAacaaca                                                                                                                                 <  6:260796/149‑1 (MQ=255)
                    tCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATCACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACCGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAacaaca                                                                                                                                 <  4:311738/149‑1 (MQ=255)
                     cGTAAGGCAACGATTAATCTTTGGCAGATTGATGCAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAAcaacaa                                                                                                                                <  2:232607/149‑1 (MQ=255)
                      gTAAGGCAACGATTAATCTTTGGCAGATTTATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTTGCAATATAATCACAACGGTTAGCTTGTTATTACGGGTGATAATGCCACAGTAAAcaacaa                                                                                                                                <  7:6200/148‑1 (MQ=255)
                      gTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAAcaacaa                                                                                                                                <  7:208523/148‑1 (MQ=255)
                      gTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGGGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAAt                                                                                                                               >  2:225798/1‑149 (MQ=255)
                       tAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATg                                                                                                                              >  4:146825/1‑149 (MQ=255)
                        aaGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATgg                                                                                                                             >  2:337262/1‑149 (MQ=255)
                         aGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATgg                                                                                                                             <  1:337262/148‑1 (MQ=255)
                          ggCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGc                                                                                                                            >  6:179552/1‑148 (MQ=255)
                           gCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAacaaca                                                                                                                                 <  3:3106/142‑1 (MQ=255)
                           gCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAacaaca                                                                                                                                 >  4:3106/1‑142 (MQ=255)
                           gCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCa                                                                                                                           >  3:130076/1‑148 (MQ=255)
                              aCGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAac                                                                                                                       >  8:298533/1‑149 (MQ=255)
                                     aTCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGtt                                                                                                                <  4:214023/149‑1 (MQ=255)
                                     aTCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACAACCGtt                                                                                                                >  3:161578/1‑149 (MQ=255)
                                     aTCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGATTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGtt                                                                                                                >  3:205408/1‑149 (MQ=255)
                                      tCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGtt                                                                                                                >  3:359091/1‑148 (MQ=255)
                                            gCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGt                                                                                                          >  1:61065/1‑148 (MQ=255)
                                                 ttGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGa                                                                                                     <  8:276038/148‑1 (MQ=255)
                                                    atgaAGCGAATAACACTGTTGCCCTTGAAGGGTTGATAGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTTATAATGCCACAGTAAACAACAATGGCAAAACCACCTTTGACGGTAAAGATTc                                                                                                  <  4:285974/148‑1 (MQ=255)
                                                     tgaAGCGAATAACACTGTTGCGCTTGTAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACTACAATGGCAAAACCACCGTTGACGGTATAGATTCCa                                                                                                <  1:114090/149‑1 (MQ=255)
                                                      gaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTGATTACGGGTGATACTGCCGCCGTAAACAACAATGGCAACACCACCGTTGACGGGAAAGATTcccc                                                                                               >  6:48900/1‑147 (MQ=255)
                                                              tAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACg                                                                                        >  8:187639/1‑148 (MQ=255)
                                                                 cacTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACTAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCCCAGTAAACAACCATGGCAAAACCACCGTTGACGGTAAAGATTCCCCCGGGACGGaaa                                                                                    >  3:240689/1‑149 (MQ=255)
                                                                        gCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTGCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACg                                                                              <  7:72772/148‑1 (MQ=255)
                                                                         cccTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAAGCAACgg                                                                             >  6:284926/1‑148 (MQ=255)
                                                                          ccTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGta                                                                           <  3:146825/149‑1 (MQ=255)
                                                                            ttGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGtaa                                                                          >  8:225734/1‑148 (MQ=255)
                                                                                     tgcgcgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                                                                >  4:229739/4‑149 (MQ=255)
                                                                                     tGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGaaa                                                                <  4:266313/149‑1 (MQ=255)
                                                                                         cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTTATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                                                  >  8:222672/1‑127 (MQ=255)
                                                                                         cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTTATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                                                  <  1:310361/127‑1 (MQ=255)
                                                                                         cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTTATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                                                  <  7:222672/127‑1 (MQ=255)
                                                                                         cgcAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTTATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATc                                                                                  >  2:310361/1‑127 (MQ=255)
                                                                                           cAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGAtt                                                          <  7:210119/149‑1 (MQ=255)
                                                                                             gATGGTGCAACGATGTGGAAATATAATCACCACGGTTAGCTTGTTATTTCGGGTGATAATGCCACAGTAATCAACAATGGCAAAACCTCCGTTGACGGTAATTATTCCACTGGTACGGAAATCATCTGTAATAACGGTAAAGtttttct                                                        <  4:377951/149‑7 (MQ=255)
                                                                                               tGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAgg                                                      <  2:199678/149‑1 (MQ=255)
                                                                                                           gTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTgg                                          <  7:298533/149‑1 (MQ=255)
                                                                                                            tGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGa                                         <  7:275602/149‑1 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAGCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                        <  7:136292/148‑1 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATc                                             <  2:272521/143‑1 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATc                                             >  1:272521/1‑143 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                        <  2:61065/148‑1 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                       >  2:26219/1‑149 (MQ=255)
                                                                                                              gCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATg                                       >  4:172033/1‑149 (MQ=255)
                                                                                                                 atatAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                        >  5:291377/1‑145 (MQ=255)
                                                                                                                 atatAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGAt                                        <  6:291377/145‑1 (MQ=255)
                                                                                                                   ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgc                                  >  1:17427/1‑149 (MQ=255)
                                                                                                                   ataATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGAATGTCAgc                                  >  8:228764/1‑149 (MQ=255)
                                                                                                                    taATCACAACGGTGAGCTTGTTATTACGGGTGATAATCCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgc                                  <  8:143900/148‑1 (MQ=255)
                                                                                                                      aTCTCAATGGTGAGCTTGTTATTACGGGTTATAATGCCACCGTAAAAAACAATGGCATAACCACCGTTTAAGGTAAAGATTCCACCGTTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAgcgg                                <  6:22468/148‑1 (MQ=18)
                                                                                                                        cacaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                             >  1:92710/1‑149 (MQ=255)
                                                                                                                         acaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcgg                             <  5:86872/148‑1 (MQ=255)
                                                                                                                           aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt                                            <  8:177491/131‑1 (MQ=255)
                                                                                                                           aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCt                                            >  7:177491/1‑131 (MQ=255)
                                                                                                                           aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                          >  4:314076/1‑149 (MQ=255)
                                                                                                                           aaCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGcggcggcgg                          >  1:214913/1‑149 (MQ=255)
                                                                                                                            aCGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGt                         >  7:57249/1‑149 (MQ=255)
                                                                                                                                tGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcgg                                >  2:47450/1‑138 (MQ=255)
                                                                                                                                tGAGCTTGGTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACg                     <  1:118093/149‑1 (MQ=255)
                                                                                                                                    cTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                 <  7:225734/149‑1 (MQ=255)
                                                                                                                                     ttGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCAACGTTGACGGTAAAGATTCAAACGGTACGTAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt                >  1:334751/1‑149 (MQ=255)
                                                                                                                                       gTTATTACGGGTGATAATGCCACTGTTAACAACAATGGAAAAACCACGGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCTGGACGGc                                                 <  4:269818/114‑1 (MQ=18)
                                                                                                                                       gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGc                                                 >  3:269818/1‑114 (MQ=255)
                                                                                                                                       gTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAtt                <  7:217034/147‑1 (MQ=255)
                                                                                                                                        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat             <  5:179552/149‑1 (MQ=255)
                                                                                                                                        ttattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat             <  1:26219/149‑1 (MQ=255)
                                                                                                                                         tattaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata            >  4:46191/1‑149 (MQ=255)
                                                                                                                                          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGTAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat             <  4:290937/147‑1 (MQ=255)
                                                                                                                                          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGatat           >  4:164125/1‑149 (MQ=255)
                                                                                                                                          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGata            <  7:187639/148‑1 (MQ=255)
                                                                                                                                          attaCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGat             >  3:290937/1‑147 (MQ=255)
                                                                                                                                                 gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGt    <  6:305188/149‑1 (MQ=255)
                                                                                                                                                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGa  <  8:53228/149‑1 (MQ=255)
                                                                                                                                                   gATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGa  >  8:59885/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGA  >  NC_000913/1465908‑1466203

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