Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,288,554 G→C pseudogene (361/2525 nt) yejO ← pseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,288,5540GC97.0% 226.7 / ‑0.4 67pseudogene (361/2525 nt)yejOpseudogene, autotransporter outer membrane homology,putative transport, Not classified, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (32/33);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑GCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688
                                                                                                                                       |                                                                                                                                       
gaCGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑Ccccc                                                                                                                                     >  1:220115/1‑139 (MQ=255)
 aCGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑cccccc                                                                                                                                    <  1:969110/139‑1 (MQ=255)
              aaTACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGAC‑‑GCCCCC‑CCATCTTTTTCAAc                                                                                                                     >  2:501486/1‑139 (MQ=255)
                   cTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGAAG‑CCCCCC‑CCATCTTTTTCCCCGAtt                                                                                                                 >  2:97520/1‑139 (MQ=255)
                          gCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  1:550302/139‑1 (MQ=255)
                               cgTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCt                                                                                                     >  2:1142134/1‑139 (MQ=255)
                                gTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCtt                                                                                                    >  2:840002/1‑139 (MQ=255)
                                     tgttTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATc                                                                                               <  1:97520/139‑1 (MQ=255)
                                                cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCAtt                                                                                    >  1:208463/1‑139 (MQ=255)
                                                cACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCAtt                                                                                    >  2:236691/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  2:259591/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  2:803447/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  2:198442/1‑139 (MQ=255)
                                                 aCCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTc                                                                                   >  2:1003719/1‑139 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:998157/1‑115 (MQ=255)
                                                  ccGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:998157/115‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCGGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGCATCATTCGCTg                                                                               <  2:874428/139‑1 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  1:406213/1‑139 (MQ=255)
                                                     ccTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:1003719/139‑1 (MQ=255)
                                                        gCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTgcgc                                                                            <  1:840002/139‑1 (MQ=255)
                                                         cTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  2:827269/1‑108 (MQ=255)
                                                         cTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  1:827269/108‑1 (MQ=255)
                                                          tttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCa                                                                          >  1:1045194/1‑139 (MQ=255)
                                                           ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  2:208463/139‑1 (MQ=255)
                                                           ttCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCAt                                                                         <  2:1206686/139‑1 (MQ=255)
                                                            tCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                                                        <  1:621413/139‑1 (MQ=255)
                                                             cTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATaa                                                                       <  1:803447/139‑1 (MQ=255)
                                                                  cTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  <  2:220115/139‑1 (MQ=255)
                                                                   tACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTg                                                                  >  2:1031063/1‑139 (MQ=255)
                                                                    aCCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTa                                                                >  2:1431434/1‑139 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              >  1:642343/1‑108 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCg                                                                                              <  2:642343/108‑1 (MQ=255)
                                                                     ccGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAg                                                               <  2:406213/139‑1 (MQ=255)
                                                                                aatCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACGCCCCCCC‑CCATCTTTTTCAACGATTAACCCCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTTTAGCACACCatta                                                     >  1:1474534/1‑139 (MQ=255)
                                                                                         tctTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAg                                           <  1:229742/139‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  1:1229142/1‑66 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  2:1229142/66‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          >  2:1389589/1‑66 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgcc                                                                                                          <  1:1389589/66‑1 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                >  1:273612/1‑76 (MQ=255)
                                                                                                   ccAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAAt                                                                                                <  2:273612/76‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGttt                                   <  2:72174/136‑1 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGttt                                   >  1:72174/1‑136 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  1:1137752/1‑139 (MQ=255)
                                                                                                    cAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGcc                                >  2:520920/1‑139 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           <  2:320488/111‑1 (MQ=255)
                                                                                                     agtagtCCTCCCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGcaca                                                           >  1:320488/1‑111 (MQ=255)
                                                                                                          tcctccCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCa                          <  1:146629/139‑1 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAcc                        >  1:1461790/1‑139 (MQ=255)
                                                                                                            ctccCTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATAATTCGCTGCGCCATAAAGCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCAcc                        >  2:970476/1‑139 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgccc                                                                                                         >  1:1451771/1‑56 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGcgccc                                                                                                         <  2:1451771/56‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  2:1442688/1‑82 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  1:1442688/82‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                                                        <  1:780328/89‑1 (MQ=255)
                                                                                                              cccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATa                                                                        >  2:780328/1‑89 (MQ=255)
                                                                                                               ccTTTTCGATAGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCAt                     <  2:1461790/139‑1 (MQ=255)
                                                                                                                    tCGATAGCGACAAAGAC‑‑GCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAAc               <  1:501486/139‑1 (MQ=255)
                                                                                                                     cGATAGCGACAAGGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAAc               >  2:188458/1‑139 (MQ=255)
                                                                                                                     cGATAGCCGCAAAGAAG‑CCCCCCACCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGaa                <  2:1474534/139‑1 (MQ=255)
                                                                                                                         aGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               >  1:359850/1‑71 (MQ=255)
                                                                                                                         aGCGACAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTg                                                                               <  2:359850/71‑1 (MQ=255)
                                                                                                                             aCAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGAt       <  1:1431434/139‑1 (MQ=255)
                                                                                                                              cAAAGACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATa      <  1:246089/139‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:672746/139‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:452610/139‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  2:1137752/139‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  1:495174/139‑1 (MQ=255)
                                                                                                                                  gACG‑CCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCg  <  1:236691/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GACGGTTTGTTCCGAATACCTCCATGGCCCGCGTGGTTGTTTTAATAGCACCGCCTGCTTTCTGATCTACCGCAAATGCAAATCCCCCCTCTTTAACCTCCAGTAGTCCTCCCTTTTCGATAGCGACAAAGACG‑GCCCCC‑CCATCTTTTTCAACGATTAAGCGCCCGCCTTTAATCGTGGTATCATTCGCTGCGCCATAAACCTGTAGCACACCATTATTGAGCGTAGTCTGGTTTGCCGCGCCACCCATTTGAACAATTTGATAACCG  >  NC_000913/2288420‑2288688

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