Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,305,438 G→A 5.3% A387V (GCG→GTG)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,305,4380GA 5.3% 243.8 / 2.3 75A387V (GCG→GTG) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (2/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

ACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGT  >  NC_000913/4305290‑4305580
                                                                                                                                                    |                                                                                                                                              
aCCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCg                                                                                                                                                >  5:48857/1‑149 (MQ=255)
 ccGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGc                                                                                                                                               <  4:281213/149‑1 (MQ=255)
  cgtcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGc                                                                                                                                               >  3:185098/1‑148 (MQ=255)
   gtcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCa                                                                                                                                             <  3:214354/149‑1 (MQ=255)
    tcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAg                                                                                                                                            <  1:199493/149‑1 (MQ=255)
    tcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAg                                                                                                                                            <  8:166208/149‑1 (MQ=255)
    tcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAg                                                                                                                                            <  2:331843/149‑1 (MQ=255)
      gtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGt                                                                                                                                          <  4:185098/149‑1 (MQ=255)
       tAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCACCAGGTg                                                                                                                                         >  3:362971/1‑149 (MQ=255)
       tAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCACCAGGTg                                                                                                                                         >  3:357126/1‑149 (MQ=255)
          tAGCCAAGATAGAAGTTATACACCGCGCGGGAGGTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAg                                                                                                                                      <  2:125232/149‑1 (MQ=255)
            gCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATCGCCGCGGCTGGCCCAGGTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGGGTAGCCCTGGTTCGCCAGGTGCAgg                                                                                                                                     <  5:124358/148‑1 (MQ=255)
            gCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAggg                                                                                                                                    >  5:176655/1‑149 (MQ=255)
            gCCAAGATAGAAGTTATACACCGCGCCGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAgg                                                                                                                                     <  8:163981/148‑1 (MQ=255)
             ccAAGATAGAAGTTATTCACCGCGGGGGCGGTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGGGCAGGCGGCAGCTTATTCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAgggg                                                                                                                                   <  5:270075/149‑2 (MQ=255)
               aaGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCt                                                                                                                                 >  1:385007/1‑149 (MQ=255)
                    aGAAGTTTTACTCCGCGCGGGCGTTGTGGCTGACTGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTc                                                                                                                             <  6:48857/148‑1 (MQ=255)
                      aaGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGt                                                                                                                           >  7:144655/1‑148 (MQ=255)
                      aaGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGt                                                                                                                           >  1:58505/1‑148 (MQ=255)
                        gTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGaa                                                                                                                        >  4:290142/1‑149 (MQ=255)
                        gTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGaa                                                                                                                        <  8:144655/149‑1 (MQ=255)
                         ttATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAAt                                                                                                                       >  8:219567/1‑149 (MQ=255)
                            taCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTa                                                                                                                    >  3:117776/1‑149 (MQ=255)
                              cacCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTAc                                                                                                                   <  4:163258/148‑1 (MQ=255)
                              cacCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTAc                                                                                                                   <  2:385007/148‑1 (MQ=255)
                               acCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACtt                                                                                                                 <  6:176655/149‑1 (MQ=255)
                                           ttGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcgcg                                                                                                     >  3:68973/1‑149 (MQ=255)
                                           ttGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcgcg                                                                                                     <  7:44196/149‑1 (MQ=255)
                                                  tGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTa                                                                                              >  2:375200/1‑149 (MQ=255)
                                                   gACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                                             >  4:139824/1‑149 (MQ=255)
                                                          ccGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGACAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                                       >  2:345248/1‑148 (MQ=255)
                                                           cGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTa                                                                                              >  8:269397/1‑140 (MQ=255)
                                                           cGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTa                                                                                              <  7:269397/140‑1 (MQ=255)
                                                             tAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAAt                                                                                    >  6:184069/1‑148 (MQ=255)
                                                                  gCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTg                                                                               <  6:95828/148‑1 (MQ=255)
                                                                  gCCGCGGCTGGCCCAGTTATGGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGa                                                                              <  5:361698/149‑1 (MQ=255)
                                                                    cgcgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGata                                                                            <  4:117776/149‑1 (MQ=255)
                                                                      cgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATg                                                                          >  8:313529/1‑149 (MQ=255)
                                                                        gCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCt                                                                        <  3:184830/149‑1 (MQ=255)
                                                                         cTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCtt                                                                       <  1:375200/149‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCACTATAATCATTGATATGCCTAtt                                                                    <  6:283081/149‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                      <  3:46479/147‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                      <  7:216491/147‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTAt                                                                     <  7:258311/148‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCACCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                      <  4:357126/147‑1 (MQ=255)
                                                                            gCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCACCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                                      <  4:362971/147‑1 (MQ=255)
                                                                                 gTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGGGAATGTACTAGATGGAATCGCGTAATTTACCAATATAATCATTGATATGCTTATAGccc                                                                >  6:12341/1‑148 (MQ=255)
                                                                                     ttGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGAc                                                           >  4:172216/1‑149 (MQ=255)
                                                                                        gCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCgg                                                        >  7:23341/1‑149 (MQ=255)
                                                                                          aaGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGc                                                       >  7:27400/1‑148 (MQ=255)
                                                                                          aaGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGc                                                       >  7:21395/1‑148 (MQ=255)
                                                                                           aGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGc                                                       >  2:138724/1‑147 (MQ=255)
                                                                                           aGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGc                                                       >  3:363770/1‑147 (MQ=255)
                                                                                           aGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCa                                                     <  1:138724/149‑1 (MQ=255)
                                                                                            gTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTTTTGCCCCAGCCCGGCCAg                                                    >  6:323173/1‑149 (MQ=255)
                                                                                                aGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGgtgt                                                >  1:63689/1‑149 (MQ=255)
                                                                                                           tAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                                      >  3:133042/1‑148 (MQ=255)
                                                                                                            aaTCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                                      >  7:312487/1‑147 (MQ=255)
                                                                                                            aaTCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                                      <  1:98070/147‑1 (MQ=255)
                                                                                                                     gCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTcc                            <  3:139824/148‑1 (MQ=255)
                                                                                                                     gCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCcg                           >  3:225577/1‑149 (MQ=255)
                                                                                                                      ccGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCcgc                          >  7:84881/1‑149 (MQ=255)
                                                                                                                      ccGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCcgc                          <  2:58505/149‑1 (MQ=255)
                                                                                                                        gATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgt                        <  5:184069/149‑1 (MQ=255)
                                                                                                                           ttcattcaTGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTAACGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgtcgt                     >  2:164347/1‑149 (MQ=255)
                                                                                                                             cattcatGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTa                   >  5:377679/1‑149 (MQ=255)
                                                                                                                                  atGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCca              >  3:170914/1‑149 (MQ=255)
                                                                                                                                     gtgtAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCcaca            <  8:27400/148‑1 (MQ=255)
                                                                                                                                     gtgtAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCcaca            <  8:21395/148‑1 (MQ=255)
                                                                                                                                        tAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGc                                                                         <  7:210229/84‑1 (MQ=255)
                                                                                                                                        tAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGc                                                                         >  8:210229/1‑84 (MQ=255)
                                                                                                                                          gCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTa                   >  1:44160/1‑136 (MQ=255)
                                                                                                                                          gCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTa                   <  2:44160/136‑1 (MQ=255)
                                                                                                                                             cTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCACAGACCGGCCAGGTTCGCGGAATAAACAGCACTTCCGCGTCGTTACCCCAAATATCCAGct   >  2:73288/1‑148 (MQ=255)
                                                                                                                                              tGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGAGATGCTTATTGCCCCAGACCGGCCAGTTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGt  >  2:225426/1‑149 (MQ=255)
                                                                                                                                                gTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGACATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGt  <  1:133461/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGTTACCCCACATATCCAGCGT  >  NC_000913/4305290‑4305580

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: ‑