Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 858,673 A→C 5.1% L130R (CTC→CGC) ‡ ybiU ← DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913858,6730AC 5.1% 229.3 / 4.4 79L130R (CTC→CGC) ‡ybiUDUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base C (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

ATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  NC_000913/858525‑858818
                                                                                                                                                    |                                                                                                                                                 
aTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                                   <  7:16733/149‑1 (MQ=255)
   cGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCAGAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTTCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                                 <  3:207027/148‑1 (MQ=255)
   cGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                                 <  2:302082/148‑1 (MQ=255)
       ttAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTccc                                                                                                                                            <  5:302555/149‑1 (MQ=255)
            ccATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaag                                                                                                                                        >  6:212381/1‑148 (MQ=255)
             cATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaaga                                                                                                                                       >  1:418732/1‑148 (MQ=255)
             cATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaaga                                                                                                                                       >  5:207462/1‑148 (MQ=255)
                   tttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATc                                                                                                                                >  8:146660/1‑149 (MQ=255)
                    ttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                          >  7:170988/1‑138 (MQ=255)
                    ttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCga                                                                                                                                          <  8:170988/138‑1 (MQ=255)
                       cATCACTTGAAAATGTCCACAGCCGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGTTGTGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCg                                                                                                                            <  6:67412/149‑1 (MQ=255)
                        aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgg                                                                                                                           >  3:131052/1‑149 (MQ=255)
                        aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgg                                                                                                                           >  8:439916/1‑149 (MQ=255)
                             tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      >  3:215731/1‑149 (MQ=255)
                             tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      <  6:55163/149‑1 (MQ=255)
                             tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      >  3:36984/1‑149 (MQ=255)
                             tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagat                                                                                                                                       <  1:501502/132‑4 (MQ=255)
                             tttCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagat                                                                                                                                       >  2:501502/1‑129 (MQ=255)
                               tCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                      <  7:146660/147‑1 (MQ=255)
                                 aaaTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTg                                                                                                                   >  1:71716/1‑148 (MQ=255)
                                       ccACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCt                                                                                                            >  4:413042/1‑149 (MQ=255)
                                         acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                          >  6:481859/1‑149 (MQ=255)
                                         acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                          >  7:487952/1‑149 (MQ=255)
                                          caGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGt                                                                                                         <  4:215731/149‑1 (MQ=255)
                                           aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTCGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                        <  2:227467/149‑1 (MQ=255)
                                           aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGt                                                                                                         >  8:32326/1‑148 (MQ=255)
                                            gACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                        >  7:381157/1‑148 (MQ=255)
                                                    agaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGtt                                                                                               <  5:481859/149‑1 (MQ=255)
                                                    agaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGtt                                                                                               <  7:439916/149‑1 (MQ=255)
                                                     gaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTg                                                                                              >  1:129683/1‑149 (MQ=255)
                                                          gACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                         <  5:212381/149‑1 (MQ=255)
                                                           aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCTACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTcc                                                                                        <  4:36984/149‑1 (MQ=255)
                                                           aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTcc                                                                                        <  8:381157/149‑1 (MQ=255)
                                                           aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                         >  7:247120/1‑148 (MQ=255)
                                                               cgcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     <  3:413042/148‑1 (MQ=255)
                                                                gcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGAGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     >  6:413216/1‑147 (MQ=255)
                                                                      gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATa                                                                                   <  2:185493/143‑1 (MQ=255)
                                                                      gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATa                                                                                   >  1:185493/1‑143 (MQ=255)
                                                                       ccATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAg                                                                             >  1:260294/1‑148 (MQ=255)
                                                                         aTTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCAt                                                                          <  6:207462/149‑1 (MQ=255)
                                                                         aTTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCAGGGTCATTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCGAGCAt                                                                          >  5:235179/1‑149 (MQ=255)
                                                                            tcttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAATTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATAGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCAAGCATCGa                                                                       >  1:113757/1‑149 (MQ=255)
                                                                             cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      >  7:377089/1‑149 (MQ=255)
                                                                              ttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                     >  5:58862/1‑149 (MQ=255)
                                                                                cACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      <  4:131052/146‑1 (MQ=255)
                                                                                  ctctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATc                                                                 >  5:106749/1‑149 (MQ=255)
                                                                                       gCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCa                                                              <  2:71716/147‑1 (MQ=255)
                                                                                           gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      <  1:417705/135‑1 (MQ=255)
                                                                                           gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                                      >  2:417705/1‑135 (MQ=255)
                                                                                             cTGCATTTGCGCCTGCGCCCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                      >  7:243815/1‑149 (MQ=255)
                                                                                             cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                             >  8:187408/1‑110 (MQ=255)
                                                                                             cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                             <  7:187408/110‑1 (MQ=255)
                                                                                             cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                      >  3:87777/1‑149 (MQ=255)
                                                                                                  tttGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 >  7:292618/1‑149 (MQ=255)
                                                                                                        ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCgg                                           <  7:32326/149‑1 (MQ=255)
                                                                                                         cTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                          >  7:460266/1‑149 (MQ=255)
                                                                                                            cGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggggg                                        <  2:418732/148‑1 (MQ=255)
                                                                                                                cAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATg                                   <  8:487952/149‑1 (MQ=255)
                                                                                                                  gTAGATGGGGTAAATCTCGTGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGACGGGGCCTTTTTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAATGATCCCAGCCTAGCGCTTGTACGCGGGGGCACTGGc                                 >  6:476506/1‑149 (MQ=255)
                                                                                                                  gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGCCCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTa                                                         >  2:118886/1‑125 (MQ=255)
                                                                                                                  gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTa                                                         <  1:118886/125‑1 (MQ=255)
                                                                                                                     gATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCt                               <  8:247120/148‑1 (MQ=255)
                                                                                                                      aTGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGATATTATCGCCGTGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCTACTGATCCCAGCCTAGCGCTTGt                                                 >  1:37701/1‑129 (MQ=255)
                                                                                                                      aTGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                                 <  2:37701/129‑1 (MQ=255)
                                                                                                                           gTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGCCTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCa                        >  3:36522/1‑149 (MQ=255)
                                                                                                                            tAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc                       <  5:413216/149‑1 (MQ=255)
                                                                                                                              aaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCg                     <  2:260294/149‑1 (MQ=255)
                                                                                                                               aTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCg                     >  7:14712/1‑148 (MQ=255)
                                                                                                                                  tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    >  7:75094/1‑130 (MQ=255)
                                                                                                                                  tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                                    <  8:75094/130‑1 (MQ=255)
                                                                                                                                    gggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAc               >  7:396284/1‑149 (MQ=255)
                                                                                                                                            aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     >  4:9202/1‑71 (MQ=255)
                                                                                                                                            aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                                     <  3:9202/71‑1 (MQ=255)
                                                                                                                                            aaGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt       >  4:178452/1‑149 (MQ=255)
                                                                                                                                              gcgcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGtt      >  1:182091/1‑148 (MQ=255)
                                                                                                                                               cgcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTa    <  6:106749/149‑1 (MQ=255)
                                                                                                                                                gcTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTaa   >  5:60475/1‑149 (MQ=255)
                                                                                                                                                 cggAGCGGCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCCGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  5:476506/147‑1 (MQ=255)
                                                                                                                                                 cTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  8:292618/149‑1 (MQ=255)
                                                                                                                                                 cTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAt  <  8:377089/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGTTTAAT  >  NC_000913/858525‑858818

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