Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 784,012 T→A 7.9% Q271L (CAG→CTG)  zitB ← zinc efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913784,0120TA 7.9% 847.7 / 47.3 239Q271L (CAG→CTG) zitBzinc efflux system
Reads supporting (aligned to +/- strand):  ref base T (111/109);  new base A (11/8);  total (122/117)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 239 total aligned reads displayed.

TCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGA  >  NC_000913/783865‑784159
                                                                                                                                                   |                                                                                                                                                   
tCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                    >  2:569571/1‑149 (MQ=255)
tCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                    >  4:743864/1‑149 (MQ=255)
 cACGTGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                   <  3:730092/149‑1 (MQ=255)
 cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                   <  3:556424/149‑1 (MQ=255)
          tGTTCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAg                                                                                                                                          <  1:1307898/149‑1 (MQ=255)
          tGTCCCATTAATGGTGATGAAGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAg                                                                                                                                          <  3:1207078/149‑1 (MQ=255)
              ccATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCACCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTATGGTTGGATCTGATCCAACAAg                                                                                                                                          >  4:121258/1‑145 (MQ=255)
                aTTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcg                                                                                                                                    >  4:782143/1‑149 (MQ=255)
                     tGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatc                                                                                                                               <  4:663259/149‑1 (MQ=255)
                         gatgatGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCgtgt                                                                                                                           >  1:945221/1‑149 (MQ=255)
                             atgTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                       <  2:930306/149‑1 (MQ=255)
                              tgtgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatc                                                                                                                               >  4:460811/1‑140 (MQ=255)
                              tgtgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCgg                                                                                                                      >  3:411649/1‑149 (MQ=255)
                                tgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                   >  2:447438/1‑118 (MQ=255)
                                  aaTGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                   <  4:822662/116‑1 (MQ=255)
                                  aaTGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                   >  3:822662/1‑116 (MQ=255)
                                    tGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaac                                                                                                                                             >  2:516749/1‑120 (MQ=255)
                                     ggCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAcc                                                                                                               >  3:382462/1‑149 (MQ=255)
                                      gCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                              >  1:58108/1‑149 (MQ=255)
                                        cAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGc                                                                                                            >  1:23807/1‑149 (MQ=255)
                                         aGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGca                                                                                                           >  2:1164465/1‑149 (MQ=255)
                                             acCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGc                                                                                                                        <  1:443728/132‑1 (MQ=255)
                                               cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtg                                                                                                                                  <  4:680129/120‑1 (MQ=255)
                                               cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCgg                                                                                                                      >  3:881194/1‑132 (MQ=255)
                                               cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                     >  1:69306/1‑149 (MQ=255)
                                               cccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                     >  4:1095473/1‑149 (MQ=255)
                                                 cTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGc                                                                                                            >  2:215850/1‑140 (MQ=255)
                                                 cTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAAc                                                                                                   >  3:1005534/1‑149 (MQ=255)
                                                 cTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAAc                                                                                                   >  2:744050/1‑149 (MQ=255)
                                                  ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACg                                                                                                  >  2:1210954/1‑149 (MQ=255)
                                                    cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtc                                                                                                <  2:1380361/149‑1 (MQ=255)
                                                    cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtc                                                                                                <  4:87407/149‑1 (MQ=255)
                                                      ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                              <  2:613495/133‑1 (MQ=255)
                                                      ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCt                                                                                                              >  1:613495/1‑133 (MQ=255)
                                                      ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                                                                              >  2:694308/1‑149 (MQ=255)
                                                      ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                                                                              >  2:163719/1‑149 (MQ=255)
                                                       tcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcatc                                                                                             >  4:602860/1‑149 (MQ=255)
                                                         aGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCAc                                                                                           >  1:1003590/1‑149 (MQ=255)
                                                           gTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                                                                         <  2:58108/149‑1 (MQ=255)
                                                             ggCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                       >  4:265468/1‑117 (MQ=255)
                                                              tcattatGGCCCATGACAAGGTTGATATTCTATCTGAATGGTGGCGTGCTCAATCTGCTAGGCATCCAACAGGAAGCGTTGGATCTGATTCAACAAGGCATCGTGATCGTGTGTAGGGACCAGCTGCACATGCAACGACATCACCGGCt                                                                                      >  4:1279435/7‑149 (MQ=255)
                                                               cAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCtt                                                                                     <  3:1095473/149‑1 (MQ=255)
                                                                aGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATg                                                                                                       >  3:792039/1‑130 (MQ=255)
                                                                        ccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAg                                                                                                                                          >  3:745809/1‑87 (MQ=255)
                                                                        ccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                            >  2:644175/1‑149 (MQ=255)
                                                                          aTGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAcc                                                                                                               >  4:918435/1‑112 (MQ=255)
                                                                            gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATcgc                                                                                                                <  1:333037/109‑3 (MQ=255)
                                                                            gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGc                                                                                                            >  3:665403/1‑113 (MQ=255)
                                                                              cAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                                                                      <  3:171397/149‑1 (MQ=255)
                                                                               aaGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATc                                                                                                                  <  1:1393505/104‑1 (MQ=255)
                                                                                 ggTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtac                                                                   <  3:585218/149‑1 (MQ=255)
                                                                                         ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtc                                                                                                >  3:986699/1‑112 (MQ=255)
                                                                                          tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                          <  1:862501/149‑1 (MQ=255)
                                                                                          tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                          >  4:167945/1‑149 (MQ=255)
                                                                                           ccATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                                                                >  1:812541/1‑126 (MQ=255)
                                                                                             aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                          >  4:537433/1‑146 (MQ=255)
                                                                                             aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                          <  3:537433/146‑1 (MQ=255)
                                                                                             aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTg                                                       >  2:1153463/1‑149 (MQ=255)
                                                                                              tCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCaa                                                                                                    <  4:58036/103‑1 (MQ=255)
                                                                                              tCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtct                                                                                   >  4:299016/1‑120 (MQ=255)
                                                                                              tCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtct                                                                                   <  3:299016/120‑1 (MQ=255)
                                                                                              tCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGc                                                      <  1:744050/149‑1 (MQ=255)
                                                                                                 gAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGcaca                                                   >  3:1174376/1‑149 (MQ=255)
                                                                                                 gAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGcaca                                                   <  1:1210954/149‑1 (MQ=255)
                                                                                                   aTGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATGTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACAtt                                                 <  3:919633/149‑1 (MQ=255)
                                                                                                      gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                              <  4:1126605/149‑1 (MQ=255)
                                                                                                        ggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGa                                            <  2:334241/149‑1 (MQ=255)
                                                                                                              cTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCg                                      >  2:769086/1‑149 (MQ=255)
                                                                                                                cAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCggg                                    <  1:644175/149‑1 (MQ=255)
                                                                                                                 aaTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                        <  1:557455/112‑1 (MQ=255)
                                                                                                                 aaTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGa                                   <  4:1301094/149‑1 (MQ=255)
                                                                                                                 aaTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGa                                   <  1:535939/149‑1 (MQ=255)
                                                                                                                    cTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGAttt                                <  2:864789/149‑1 (MQ=255)
                                                                                                                     tGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGa                                   >  1:673533/1‑145 (MQ=255)
                                                                                                                     tGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGa                                   <  2:673533/145‑1 (MQ=255)
                                                                                                                      gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTcc                              <  2:251660/149‑1 (MQ=255)
                                                                                                                         aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCgg                           <  1:1164465/149‑1 (MQ=255)
                                                                                                                          gTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGcaca                                                   <  2:1396405/124‑1 (MQ=255)
                                                                                                                             aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                                                                <  1:160908/92‑1 (MQ=255)
                                                                                                                                cATCAGGTAGTGTTGGATCTGATCCAACAAGGGATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCGCGGCACATac                    >  4:932100/1‑149 (MQ=255)
                                                                                                                                cATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATac                    >  2:312544/1‑149 (MQ=255)
                                                                                                                                  tCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacga                  <  4:1424178/149‑1 (MQ=255)
                                                                                                                                  tCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacga                  >  3:1424178/1‑149 (MQ=255)
                                                                                                                                     tgccttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCAGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGc               >  4:871048/6‑149 (MQ=255)
                                                                                                                                     ggTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcca                                                                      >  1:917605/1‑94 (MQ=255)
                                                                                                                                      gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                              <  1:5028/117‑1 (MQ=255)
                                                                                                                                      gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                              >  2:5028/1‑117 (MQ=255)
                                                                                                                                      gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGAc                                           >  4:794588/1‑120 (MQ=255)
                                                                                                                                       ccttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCAGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCtt             <  3:871048/146‑1 (MQ=255)
                                                                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                                                                         <  3:1298651/70‑1 (MQ=255)
                                                                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc            >  1:176605/3‑149 (MQ=255)
                                                                                                                                        cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCACCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc            >  3:649184/3‑149 (MQ=255)
                                                                                                                                         ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca           >  4:510716/2‑149 (MQ=255)
                                                                                                                                         ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca           >  3:220922/2‑149 (MQ=255)
                                                                                                                                         ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca           >  1:283868/2‑149 (MQ=255)
                                                                                                                                         ttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca           <  4:220922/148‑1 (MQ=255)
                                                                                                                                            ttGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                              <  1:267023/111‑1 (MQ=255)
                                                                                                                                             tGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagttc       <  1:1153463/149‑1 (MQ=255)
                                                                                                                                               gATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGc                                                                                       <  2:848362/67‑1 (MQ=255)
                                                                                                                                                 tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCg   >  3:1372027/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGA  >  NC_000913/783865‑784159

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