Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,342,169 C→A 5.1% H83Q (CAC→CAA)  pyrF → orotidine‑5'‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,342,1690CA 5.1% 288.5 / 5.0 79H83Q (CAC→CAA) pyrForotidine‑5'‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base A (2/2);  total (36/43)
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.69e-01

ACAAGATCGACCCACGCGATTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGC  >  NC_000913/1342021‑1342316
                                                                                                                                                    |                                                                                                                                                   
aCAAGATCGACCCACGCGATTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCAc                                                                                                                                                     >  4:460436/1‑149 (MQ=255)
      tCGACCCACGCGATTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTc                                                                                                                                               <  2:480534/149‑1 (MQ=255)
              cgcgATTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCtgctgc                                                                                                                                       >  2:621178/1‑149 (MQ=255)
                  aTTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCt                                                                                                                                   >  3:56597/1‑149 (MQ=255)
                     gtcgtcTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGAc                                                                                                                                <  4:225774/149‑1 (MQ=255)
                      tcgtcTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACt                                                                                                                               >  4:205581/1‑149 (MQ=255)
                         tcTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAg                                                                                                                            >  3:66823/1‑149 (MQ=255)
                            gAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCg                                                                                                                         >  3:303979/1‑149 (MQ=255)
                             aaGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCgt                                                                                                                        <  1:149291/149‑1 (MQ=255)
                               ggTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTa                                                                                                                             <  4:184741/142‑1 (MQ=255)
                               ggTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTa                                                                                                                             >  3:184741/1‑142 (MQ=255)
                               ggTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCgtgt                                                                                                                      <  1:227957/149‑1 (MQ=255)
                                gTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGtgtg                                                                                                                     <  1:559686/149‑1 (MQ=255)
                                  cGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGa                                                                                                                   <  2:108709/149‑1 (MQ=255)
                                   ggCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGAt                                                                                                                  <  4:280380/149‑1 (MQ=255)
                                     cAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATgg                                                                                                                <  1:101990/149‑1 (MQ=255)
                                       aaGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTg                                                                                                              <  1:465826/149‑1 (MQ=255)
                                       aaGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTg                                                                                                              >  3:501910/1‑149 (MQ=255)
                                          agaTGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAAt                                                                                                           <  4:585225/149‑1 (MQ=255)
                                            aTGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGt                                                                                                         <  4:200981/149‑1 (MQ=255)
                                               tttACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCa                                                                                                      >  2:137259/1‑149 (MQ=255)
                                                ttACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCAt                                                                                                     <  4:518588/149‑1 (MQ=255)
                                                ttACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCAt                                                                                                     <  3:205581/149‑1 (MQ=255)
                                                     ttGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCctc                                                                                                >  3:542954/1‑149 (MQ=255)
                                                       gTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCtg                                                                                              >  1:136389/1‑149 (MQ=255)
                                                           gggCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCtggtgg                                                                                          >  3:132230/1‑149 (MQ=255)
                                                           gggCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCtggtgg                                                                                          <  1:621178/149‑1 (MQ=255)
                                                             gCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTgggg                                                                                        <  4:66823/149‑1 (MQ=255)
                                                                acaGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGgcgc                                                                                     >  3:380042/1‑149 (MQ=255)
                                                                acaGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGgcgc                                                                                     <  4:250177/149‑1 (MQ=255)
                                                                    tttGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTat                                                                                 <  2:112398/149‑1 (MQ=255)
                                                                        tgCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCAAGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGc                                                                                                   >  4:336166/1‑127 (MQ=255)
                                                                        tgCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCAAGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGc                                                                                                   <  3:620170/127‑1 (MQ=255)
                                                                        tgCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCAAGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGc                                                                                                   <  3:336166/127‑1 (MQ=255)
                                                                        tgCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCAAGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGc                                                                                                   >  4:620170/1‑127 (MQ=255)
                                                                               aCTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAg                                                                      >  4:135832/1‑149 (MQ=255)
                                                                                cTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGAc                                                                           >  1:332042/1‑143 (MQ=255)
                                                                                cTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGAc                                                                           <  2:332042/143‑1 (MQ=255)
                                                                                 ttCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAATTTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAgcg                                                                    >  4:269644/1‑149 (MQ=255)
                                                                                 ttCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAgcg                                                                    >  3:107940/1‑149 (MQ=255)
                                                                                  tCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAgcgc                                                                   <  1:618063/149‑1 (MQ=255)
                                                                                    aaCAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTAtgatga                                                                            <  4:545815/138‑1 (MQ=255)
                                                                                    aaCAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTAtgatga                                                                            >  3:545815/1‑138 (MQ=255)
                                                                                     aCAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTg                                                                >  4:493835/1‑149 (MQ=255)
                                                                                       aGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAg                                                              <  1:137259/149‑1 (MQ=255)
                                                                                                  gATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTcc                                                   <  2:338225/149‑1 (MQ=255)
                                                                                                  gATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTcc                                                   <  4:542954/149‑1 (MQ=255)
                                                                                                    tatCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGt                                                 <  1:116940/149‑1 (MQ=255)
                                                                                                     atCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGtt                                                >  3:425258/1‑149 (MQ=255)
                                                                                                     atCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGtt                                                <  4:56597/149‑1 (MQ=255)
                                                                                                     atCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGtt                                                <  3:135832/149‑1 (MQ=255)
                                                                                                      tctttcttGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGttt                                               <  3:460436/149‑1 (MQ=255)
                                                                                                       ctttcttGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTg                                              >  1:346743/1‑149 (MQ=255)
                                                                                                            ttGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCaaa                                         >  3:190822/1‑149 (MQ=255)
                                                                                                            ttGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCaaa                                         <  2:378907/149‑1 (MQ=255)
                                                                                                                   gAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCAc                                  >  1:519465/1‑149 (MQ=255)
                                                                                                                   gAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGCGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCAc                                  >  4:217691/1‑149 (MQ=255)
                                                                                                                      aTTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGc                               >  3:33349/1‑149 (MQ=255)
                                                                                                                       ttCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCAc                                                          >  3:201042/1‑121 (MQ=255)
                                                                                                                       ttCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCAc                                                          <  4:201042/121‑1 (MQ=255)
                                                                                                                       ttCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCt                              >  1:19685/1‑149 (MQ=255)
                                                                                                                          cACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCtttt                           <  4:190822/149‑1 (MQ=255)
                                                                                                                          cACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATCACCCCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACTGCtttt                           <  2:310557/149‑1 (MQ=255)
                                                                                                                           aCGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTg                          >  1:135133/1‑149 (MQ=255)
                                                                                                                                 tCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCt                    <  4:33349/149‑1 (MQ=255)
                                                                                                                                   cccAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCtgt                  >  3:104261/1‑149 (MQ=255)
                                                                                                                                    ccAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCtgtg                 >  4:276585/1‑149 (MQ=255)
                                                                                                                                     cAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGa                <  1:213408/149‑1 (MQ=255)
                                                                                                                                      aaCACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGAc               >  3:389887/1‑149 (MQ=255)
                                                                                                                                         acTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAgt            <  4:501910/149‑1 (MQ=255)
                                                                                                                                         acTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAgt            >  4:194165/1‑149 (MQ=255)
                                                                                                                                           tGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAgtgt          <  1:25353/149‑1 (MQ=255)
                                                                                                                                            gCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGtt         <  2:582115/149‑1 (MQ=255)
                                                                                                                                            gCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGtt         <  2:136389/149‑1 (MQ=255)
                                                                                                                                            gCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGtt         <  4:132230/149‑1 (MQ=255)
                                                                                                                                                   aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGa                         >  3:181888/1‑126 (MQ=255)
                                                                                                                                                   aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGa                         <  4:181888/126‑1 (MQ=255)
                                                                                                                                                   aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGc  >  2:496670/1‑149 (MQ=255)
                                                                                                                                                   aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGc  >  4:295623/1‑149 (MQ=255)
                                                                                                                                                   aCGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGc  <  2:361621/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACAAGATCGACCCACGCGATTGTCGTCTGAAGGTCGGCAAAGAGATGTTTACATTGTTTGGGCCACAGTTTGTGCGCGAACTTCAACAGCGTGGTTTTGATATCTTTCTTGACCTGAAATTCCACGATATCCCCAACACTGCAGCGCACGCTGTCGCTGCTGCAGCTGACTTAGGCGTGTGGATGGTGAATGTTCATGCCTCTGGTGGGGCGCGTATGATGACCGCAGCGCGTGAGGCACTGGTTCCGTTTGGCAAAGATGCACCGCTTTTGATTGCTGTGACAGTGTTGACCAGC  >  NC_000913/1342021‑1342316

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