Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,876 C→A 6.2% R163S (AGG→AGT ykgH ← putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,8760CA 6.2% 341.4 / 10.2 96R163S (AGG→AGTykgHputative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (47/43);  new base A (3/3);  total (50/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 = 9.66e-01

CATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGCA  >  NC_000913/324729‑325023
                                                                                                                                                   |                                                                                                                                                   
cATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                    >  2:151291/1‑149 (MQ=255)
  tttttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTg                                                                                                                                                  >  1:842169/1‑149 (MQ=255)
   ttttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTgg                                                                                                                                                 >  3:1434952/1‑149 (MQ=255)
        tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAgg                                                                                                                                            <  1:151291/149‑1 (MQ=255)
        tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCatt                                                                                                                                            <  1:334550/149‑3 (MQ=255)
                         aCAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGtata                                                                                                                           >  3:149759/1‑149 (MQ=255)
                           aGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGtat                                                                                                                                         <  4:20221/133‑1 (MQ=255)
                           aGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGtat                                                                                                                                         >  3:20221/1‑133 (MQ=255)
                           aGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaaa                                                                                                                         >  4:873540/1‑149 (MQ=255)
                           aGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaaa                                                                                                                         <  4:1434952/149‑1 (MQ=255)
                            gAAATTACAATGATTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAg                                                                                                                        >  2:366334/1‑149 (MQ=255)
                                ttACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGt                                                                                                                                                >  3:283167/1‑121 (MQ=255)
                                ttACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGt                                                                                                                                                <  4:283167/121‑1 (MQ=255)
                                ttACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCa                                                                                                                    <  4:149759/149‑1 (MQ=255)
                                  aCAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTTTTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACt                                                                                                                  <  1:366334/149‑1 (MQ=255)
                                         ttAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAAtt                                                                                                           <  3:873540/149‑1 (MQ=255)
                                              aaaaaTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAg                                                                                                                                             >  1:15933/1‑110 (MQ=255)
                                              aaaaaTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAg                                                                                                                                             <  2:15933/110‑1 (MQ=255)
                                                  aTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGc                                                                                                  >  2:103953/1‑149 (MQ=255)
                                                     ccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                               <  2:904264/149‑1 (MQ=255)
                                                     ccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                               >  1:904264/1‑149 (MQ=255)
                                                        gagCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGt                                                                                            <  2:842169/149‑1 (MQ=255)
                                                                               tatttattTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctgcc                                                                     <  1:103953/149‑1 (MQ=255)
                                                                                  ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCa                                                                                                                                     >  4:264829/1‑82 (MQ=255)
                                                                                  ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCa                                                                                                                                     <  3:264829/82‑1 (MQ=255)
                                                                                  ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGttt                                                                                              <  4:1455299/121‑1 (MQ=255)
                                                                                  ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGttt                                                                                              >  3:1455299/1‑121 (MQ=255)
                                                                                  ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCAAATCGAGAACTTTCCTGCCAtt                                                                  >  1:1367998/1‑149 (MQ=255)
                                                                                     ttttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATg                                                               >  2:289215/1‑149 (MQ=255)
                                                                                         ttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgc                                                                          <  4:1443192/134‑1 (MQ=255)
                                                                                         ttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgc                                                                          >  3:1443192/1‑134 (MQ=255)
                                                                                            tttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                        <  2:883151/101‑1 (MQ=255)
                                                                                            tttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                        >  1:883151/1‑101 (MQ=255)
                                                                                               cTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     >  3:701602/1‑149 (MQ=255)
                                                                                                    gAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATa                                                >  2:757047/1‑149 (MQ=255)
                                                                                                       cTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                    >  1:602560/1‑110 (MQ=255)
                                                                                                       cTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCCGAAGGCGTTTGTATCTCATc                                                                                    <  2:602560/110‑1 (MQ=255)
                                                                                                          ataGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     <  4:1206439/138‑1 (MQ=255)
                                                                                                          ataGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     >  3:1206439/1‑138 (MQ=255)
                                                                                                          ataGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCAttt                                          >  1:55515/1‑149 (MQ=255)
                                                                                                            aGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                    >  2:201180/1‑105 (MQ=255)
                                                                                                            aGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                    <  1:201180/105‑1 (MQ=255)
                                                                                                             gTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctgc                                                                      >  4:1323950/1‑118 (MQ=255)
                                                                                                             gTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctgc                                                                      <  3:1323950/118‑1 (MQ=255)
                                                                                                              ttACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                    <  3:1219767/103‑1 (MQ=255)
                                                                                                              ttACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATc                                                                                    >  4:1219767/1‑103 (MQ=255)
                                                                                                              ttACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGata                                      >  4:771050/1‑149 (MQ=255)
                                                                                                                acacAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGaa                                                                                                     <  2:983319/84‑1 (MQ=255)
                                                                                                                acacAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGaa                                                                                                     >  1:983319/1‑84 (MQ=255)
                                                                                                                 cacaGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTc                                                         >  1:712724/1‑127 (MQ=255)
                                                                                                                 cacaGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTc                                                         <  2:712724/127‑1 (MQ=255)
                                                                                                                  acaGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATag                                  >  2:1052897/1‑149 (MQ=255)
                                                                                                                    aGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATagag                                <  1:757047/149‑1 (MQ=255)
                                                                                                                    aGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATagag                                >  1:979956/1‑149 (MQ=255)
                                                                                                                       cATCACCGTTCGGATACAATGACCAGTTCCTGGTCa                                                                                                                                              <  4:1052041/36‑1 (MQ=255)
                                                                                                                       cATCACCGTTCGGATACAATGACCAGTTCCTGGTCa                                                                                                                                              >  3:1052041/1‑36 (MQ=255)
                                                                                                                       cATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgc                                                                          <  1:1064524/104‑1 (MQ=255)
                                                                                                                       cATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgc                                                                          >  2:1064524/1‑104 (MQ=255)
                                                                                                                          cACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAttttt                                                                                                                                >  4:585591/1‑47 (MQ=255)
                                                                                                                          cACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAttttt                                                                                                                                <  3:585591/47‑1 (MQ=255)
                                                                                                                           aCCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCga                                                                                  >  2:1424480/1‑92 (MQ=255)
                                                                                                                           aCCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCga                                                                                  <  1:1424480/92‑1 (MQ=255)
                                                                                                                           aCCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCa                         >  1:1395766/1‑149 (MQ=255)
                                                                                                                              gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCa                                                                                                                    <  1:956518/55‑1 (MQ=255)
                                                                                                                              gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCa                                                                                                                    >  2:956518/1‑55 (MQ=255)
                                                                                                                              gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGa                                        >  4:489671/1‑131 (MQ=255)
                                                                                                                              gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGa                                        <  3:489671/131‑1 (MQ=255)
                                                                                                                                   gATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGat                 >  3:741600/1‑149 (MQ=255)
                                                                                                                                    aTACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGata                >  1:893113/1‑149 (MQ=255)
                                                                                                                                     tACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGatat               <  2:979956/149‑1 (MQ=255)
                                                                                                                                      tacctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAgg                               >  2:511055/5‑132 (MQ=255)
                                                                                                                                      tacctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAgg                               <  1:511055/128‑1 (MQ=255)
                                                                                                                                       cAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATAttt             <  2:1395766/149‑1 (MQ=255)
                                                                                                                                       cAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATAttt             <  2:893113/149‑1 (MQ=255)
                                                                                                                                        cctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTg                                                                                                          <  1:169959/53‑1 (MQ=37)
                                                                                                                                        cctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTg                                                                                                          >  2:169959/3‑55 (MQ=37)
                                                                                                                                        aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGTTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTc            >  2:492409/1‑149 (MQ=255)
                                                                                                                                        aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGt                                                                                                >  3:1067748/1‑65 (MQ=255)
                                                                                                                                        aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGt                                                                                                <  4:1067748/65‑1 (MQ=255)
                                                                                                                                        aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAc                                               >  3:43458/1‑114 (MQ=255)
                                                                                                                                        aaTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAc                                               <  4:43458/114‑1 (MQ=255)
                                                                                                                                         ctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc           >  2:315460/2‑149 (MQ=255)
                                                                                                                                         aTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc           <  3:771050/149‑1 (MQ=255)
                                                                                                                                         aTGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc           <  2:1367998/149‑1 (MQ=255)
                                                                                                                                          tGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAt                                                                                          >  1:221064/1‑69 (MQ=255)
                                                                                                                                          tGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAt                                                                                          <  2:221064/69‑1 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTa                                                                 <  1:615869/91‑1 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTa                                                                 >  2:615869/1‑91 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     >  3:225287/1‑103 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     <  4:225287/103‑1 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGatat               >  4:1060208/1‑141 (MQ=255)
                                                                                                                                             ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGatat               <  3:1060208/141‑1 (MQ=255)
                                                                                                                                              cAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTc      >  4:303072/1‑149 (MQ=255)
                                                                                                                                               aGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTca     <  1:289215/149‑1 (MQ=255)
                                                                                                                                                gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     >  1:968603/1‑100 (MQ=255)
                                                                                                                                                gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGAc                                                     <  2:968603/100‑1 (MQ=255)
                                                                                                                                                  tCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagca  <  1:1052897/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGCA  >  NC_000913/324729‑325023

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: