Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,467,165 G→A 5.2% pseudogene (1774/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,467,1650GA 5.2% 238.7 / 2.6 78pseudogene (1774/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (29/44);  new base A (2/2);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAAC  >  NC_000913/1467018‑1467311
                                                                                                                                                   |                                                                                                                                                  
ggCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGc                                                                                                                                                   <  1:360237/149‑1 (MQ=255)
    ttGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACa                                                                                                                                               >  3:21460/1‑149 (MQ=255)
       ttaaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGTACATTGGTGACAACGGCACGGTTAACAACAACGGTGATATTGTTGATTCAGATACTGGATCGATTGGCGTGCTCAGCAACGGTGAGGGGCCAACAGTa                                                                                                                                            >  2:14465/3‑149 (MQ=255)
       tAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTa                                                                                                                                            <  2:263240/149‑1 (MQ=255)
        aaaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTAt                                                                                                                                           <  4:386973/149‑1 (MQ=255)
         aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATc                                                                                                                                          <  1:676736/149‑1 (MQ=255)
         aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATc                                                                                                                                          <  1:354377/149‑1 (MQ=255)
              ggATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATa                                                                                                                                     >  3:341572/1‑149 (MQ=255)
               gATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATAc                                                                                                                                    <  1:206110/149‑1 (MQ=255)
                aTGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACg                                                                                                                                   <  4:638794/149‑1 (MQ=255)
                 tGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACgg                                                                                                                                  <  4:34574/149‑1 (MQ=255)
                   gTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGt                                                                                                                                >  3:548018/1‑149 (MQ=255)
                    tAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTg                                                                                                                               >  2:514877/1‑149 (MQ=255)
                         tGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGtt                                                                                                                          <  1:31717/149‑1 (MQ=255)
                            cTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAAc                                                                                                                       <  4:304843/149‑1 (MQ=255)
                               ttACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGtt                                                                                                                    <  4:114893/149‑1 (MQ=255)
                                tACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTa                                                                                                                   <  1:430978/149‑1 (MQ=255)
                                 aCAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAg                                                                                                                  <  1:396730/149‑1 (MQ=255)
                                               tGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAg                                                                                                    <  3:138444/149‑1 (MQ=255)
                                               tGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAg                                                                                                    <  4:329036/149‑1 (MQ=255)
                                                ggTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAgg                                                                                                   <  3:245978/149‑1 (MQ=255)
                                                 gtgtTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAggg                                                                                                  <  2:168510/149‑1 (MQ=255)
                                                 gtgtTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAggg                                                                                                  <  2:441904/149‑1 (MQ=255)
                                                  tgtTGGGAACATTGGTGACAACGGCACGGTTTACAACAAAGGCGATATTGTTGTTTCCGATACTGGATCGATTTGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGt                                                                                                 >  2:171513/1‑149 (MQ=255)
                                                   gtTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGtt                                                                                                <  4:385557/149‑1 (MQ=255)
                                                    ttGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTc                                                                                               <  1:603948/149‑1 (MQ=255)
                                                        gaACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCa                                                                                           >  1:218569/1‑149 (MQ=255)
                                                          aCATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATc                                                                                         <  2:255363/149‑1 (MQ=255)
                                                               ggTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCAcaac                                                                                    <  4:363816/149‑1 (MQ=255)
                                                                   aCAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGACAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCAcaaccaac                                                                                >  1:283483/1‑149 (MQ=255)
                                                                     aaCGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGaa                                                                                                           <  3:176358/120‑1 (MQ=255)
                                                                     aaCGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGaa                                                                                                           >  4:176358/1‑120 (MQ=255)
                                                                          cACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGACAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTg                                                                            <  3:613762/146‑1 (MQ=255)
                                                                          cACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGACAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTg                                                                            >  4:613762/1‑146 (MQ=255)
                                                                           aCGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGaa                                                                        >  3:423737/1‑149 (MQ=255)
                                                                            aggTTAACAACAAAGGAGATATTGTTGTTTCCGATAATGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAg                                                                       <  4:154202/148‑1 (MQ=255)
                                                                             ggTTAACAACAAAGGAGATATTGTTGTATCCGATACTGGATCGAGTGGCGTGCTCATTAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGTGTTCAGCATCACAACCAACAGTGGGAAgg                                                                      >  2:264352/1‑149 (MQ=255)
                                                                               ttAACAACAAAGGATATATTGTTGTTTCCGATACTAGAGCGATTAGCGTGCGCATCCACGGTGAGGGGGCAACAGTATCCAATACGGGGGATGTTAAGGCGAGCAATGATGCGACAGGGTTCAGCATCGCAACCAACAGTGGGAAgggt                                                                    >  2:266473/1‑147 (MQ=255)
                                                                                 aacaacAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGACAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTc                                                                  <  2:283483/149‑1 (MQ=255)
                                                                                  acaacaAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCg                                                                 >  4:22649/1‑149 (MQ=255)
                                                                                   caacaaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGc                                                                >  4:638437/1‑149 (MQ=255)
                                                                                    aacaaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCt                                                               >  4:161079/1‑149 (MQ=255)
                                                                                    aacaaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCt                                                               >  3:251490/1‑149 (MQ=255)
                                                                                      caaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCAcaaccaac                                                                                <  1:587317/130‑1 (MQ=255)
                                                                                      caaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTgg                                                             >  4:283564/1‑149 (MQ=255)
                                                                                      caaAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGATAACGTTAGCAATGAAGCGACAGGGTTCAGCATCAcaaccaac                                                                                >  2:587317/1‑130 (MQ=255)
                                                                                         aGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTggcag                                                          >  1:349177/1‑149 (MQ=255)
                                                                                           gagaTATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTggcaggc                                                        <  1:267756/149‑1 (MQ=255)
                                                                                                   gttgttTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCa                                                <  2:202085/149‑1 (MQ=255)
                                                                                                   gttgttTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCACCACAACCAACAGGGGGAAGGTTTCGCTGGCAGGCAGTATGCa                                                >  4:97759/1‑149 (MQ=255)
                                                                                                    ttgttTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGCGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAg                                               >  3:389306/1‑149 (MQ=255)
                                                                                                      gttTCCGATACTGGATCGATTGGCGTGCTAATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGt                                             <  1:196500/149‑1 (MQ=255)
                                                                                                       tttCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGtt                                            >  3:453154/1‑149 (MQ=255)
                                                                                                           cGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCg                                        >  3:7341/1‑149 (MQ=255)
                                                                                                            gATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAgg                                                                                                   >  4:656329/1‑89 (MQ=255)
                                                                                                            gATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAgg                                                                                                   <  3:656329/89‑1 (MQ=255)
                                                                                                            gATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGa                                       >  3:417316/1‑149 (MQ=255)
                                                                                                             aTACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCGAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGGGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGAt                                      <  4:501573/149‑1 (MQ=255)
                                                                                                             aTACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGAt                                      <  1:514877/149‑1 (MQ=255)
                                                                                                               aCTGGATAGAATGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGAttt                                    <  2:48734/149‑1 (MQ=255)
                                                                                                                 tGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTtct                                  <  3:643312/149‑1 (MQ=255)
                                                                                                                  ggATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTtctc                                 <  2:419909/149‑1 (MQ=255)
                                                                                                                    aTCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGa                               <  3:381850/149‑1 (MQ=255)
                                                                                                                    aTCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGa                               >  1:419299/1‑149 (MQ=255)
                                                                                                                         ttGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCggg                          >  3:417589/1‑149 (MQ=255)
                                                                                                                          tGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCgggg                         <  2:518014/149‑1 (MQ=255)
                                                                                                                            gCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTa                       >  2:204285/1‑149 (MQ=255)
                                                                                                                                gCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATc                   <  4:417589/149‑1 (MQ=255)
                                                                                                                                gCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATc                   <  3:124527/149‑1 (MQ=255)
                                                                                                                                  tcatcaACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCtt                 >  1:341777/1‑149 (MQ=255)
                                                                                                                                  tcatcaACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCtt                 >  1:694981/1‑149 (MQ=255)
                                                                                                                                     tcaACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAAt              <  2:8111/149‑1 (MQ=255)
                                                                                                                                     tcaACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAAt              <  1:564880/149‑1 (MQ=255)
                                                                                                                                     tcaACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAAt              <  4:417316/149‑1 (MQ=255)
                                                                                                                                      caACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATg             >  2:687716/1‑149 (MQ=255)
                                                                                                                                          ggTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGcaa         <  4:317045/149‑1 (MQ=255)
                                                                                                                                          ggTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGcaa         <  4:119862/149‑1 (MQ=255)
                                                                                                                                           gTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGcaac        <  2:639743/149‑1 (MQ=255)
                                                                                                                                               gggggCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATa    >  4:91400/1‑149 (MQ=255)
                                                                                                                                               gggggCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCCTAATGGCAACAATa    >  2:116869/1‑149 (MQ=255)
                                                                                                                                                ggggCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATaa   >  4:358554/1‑149 (MQ=255)
                                                                                                                                                ggggCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATaa   >  1:138953/1‑149 (MQ=255)
                                                                                                                                                 gggCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAAc  <  3:91400/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGACCGGGGTAGATCTTAATGGCAACAATAAC  >  NC_000913/1467018‑1467311

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